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)。
谁能帮忙呢？
请用具体程序帮忙解答，谢谢！
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)。
谁能帮忙呢？
请用具体程序帮忙解答，谢谢！