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

问题同上
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的绘图了,该如何处理

 采纳的回答

colorw
colorw 2022-11-17

0 个投票

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量级积分是有值的

更多回答(0 个)

类别

帮助中心File Exchange 中查找有关 Mathematics and Optimization 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!