查看单个帖子
旧 2008-09-13, 13:14   #1
yiyu590
初级会员
 
注册日期: 2008-09-13
年龄: 44
帖子: 1
声望力: 0
yiyu590 正向着好的方向发展
默认 [求助]有关卷积和微分方程求解画图的问题

有两个函数如下:
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)。

谁能帮忙呢?
请用具体程序帮忙解答,谢谢!
yiyu590 当前离线   回复时引用此帖