題主給出的積分函數,用matlab 求解,可以這樣做
1、確定x在-300,300中,取若幹個x(i)數據,如
x=-300:10:300;
2、使用for循環語句,將x(i)代入積分函數I(x)中,使用int函數計算其積分值,同時將值賦值給Ix數組變量中
for n=1:600/t+1
syms x
A_1=int(cos((pi*(x^2))/2),x,a1(1,n),a2(1,n));
A_2=int(sin((pi*(x^2))/2),x,a1(1,n),a2(1,n));
Ix(n)=1000*(lambda/(2*L))*((A_1)^2+(A_2)^2);
end
3、最後用plot函數,繪出其x—I(x)的圖形
運行程序,可以得到如下結果。