函数不可积的情况下如何对其原函数进行绘图。

问题同上

syms psi t;

r0=38.9971;

r1=4.5580;

R=(r0+r1)*cos(pi/14+psi)+sqrt(100-((r0+r1)*sin(pi/14+psi)).^2);

S=8.1324*10^-8/(pi*R*R);

c=3*10^8;

t1=(t-(2*R)/c)/(10*10^-9);

a=pi/3.6*10^-4;

F=10^6*t1*exp(-t1^2-2*(tan(psi))^2/a^2);

G=F*S;

int (G,psi,-a,a)

用int并不能求出原函数，更加不能进行后面对的G的绘图了，该如何处理

Answer 1

[colorw]

syms psi t;
r0=38.9971;
r1=4.5580;
R=(r0+r1)*cos(pi/14+psi)+sqrt(100-((r0+r1)*sin(pi/14+psi)).^2);
S=8.1324*10^-8/(pi*R*R);
c=3*10^8;
t1=(t-(2*R)/c)/(10*10^-9);
a=pi/3.6*10^-4;
F=10^6*t1*exp(-t1^2-2*(tan(psi))^2/a^2);

G=F*S;
G=matlabFunction(G);
J=@(t)integral(@(psi)G(psi,t),-a,a);
t=linspace(-1e-6,1e-6,1e3);
y=arrayfun(J,t);
plot(t,y)

t在10^-6量级积分是有值的