MATLAB爱好者论坛-LabFans.com

MATLAB爱好者论坛-LabFans.com (https://www.labfans.com/bbs/index.php)
-   MATLAB论坛 (https://www.labfans.com/bbs/forumdisplay.php?f=6)
-   -   [求助]有关卷积和微分方程求解画图的问题 (https://www.labfans.com/bbs/showthread.php?t=4442)

yiyu590 2008-09-13 13:14

[求助]有关卷积和微分方程求解画图的问题
 
有两个函数如下:
o=exp(-pi*(t/0.3)^2);
l=exp(-t/0.6)*heaviside(t);
怎样求这两个函数的卷积呢?


我试着求了一下,程序和结果如下:
syms t r;
o=exp(-pi*(t/0.3)^2);
l=exp(-t/0.6)*heaviside(t);
g=int(exp(-pi*((t-r)/0.3)^2)*exp(-r/0.6)*heaviside(r),r,-inf,inf);
得到的结果为
g =
4810910854423451/295147905179352825856+3/40000*erf(1/12*(40*pi*t-3)/pi^(1/2))*exp(-1/48*(80*pi*t-3)/pi)+3/40000*exp(-1/48*(80*pi*t-3)/pi)

结果中含有一个误差函数。

请问,是不是我的解法错了,怎样才能使结果中没有误差函数呢?
其实,我知道可以用conv2求卷积的,但是结果中还是有误差函数吧。

第二个问题,我要讲得到的g带入非线性微分方程中求出解的曲线。方程如下:

Dy+(1+85*g/(85*2.3*10^(-15))*y=8.6/(85*2.3*10^(-15)。

谁能帮忙呢?
请用具体程序帮忙解答,谢谢!


所有时间均为北京时间。现在的时间是 12:38

Powered by vBulletin
版权所有 ©2000 - 2025,Jelsoft Enterprises Ltd.