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

3 次查看（过去 30 天）
gowat
gowat 2022-11-17
回答： colorw 2022-11-17
问题同上
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
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 个）

类别

在 Help Center 和 File Exchange 中查找有关 Logical 的更多信息

标签

Community Treasure Hunt

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

Start Hunting!