求助一个信号处理中的微积分问题
mathematica吧
全部回复
仅看楼主
level 5
陈醋鲈鱼1 楼主
要求如下图所示
就是定义了两个信号和两个多变量函数,要求计算两信号对应的多变量函数在一定区间上积分所得的关于t的单变量函数,并画出相关图象
我初步写了两段代码但是都报错或者有问题,麻烦大佬帮我看一下
f0=150;
m=400;
T=1;
f1=550;
k=(f0-f1)/(f0*f1*T);
fs=2500;
Ts=1/fs;
xlfm[t_]:=Exp[I*2Pi*(f0*t+m*t^2*0.5)];
xhfm[t_]:=Exp[I*2Pi*((1/k)Log[1+k*f0*t])];
Xnblfm=Integrate[xlfm[t]*Conjugate[xlfm[t-tao]]*Exp[2Pi*I*fd*t],{tao,-Infinity,Infinity},{fd,-100,100}];
Xwblfm=Sqrt[Abs[a]]*Integrate[xlfm[t]*Conjugate[xlfm[a(t-tao)]],{tao,-Infinity,Infinity},{a,0.8,1.2}];
Xnhfm=Integrate[xhfm[t]*Conjugate[xhfm[t-tao]]*Exp[2Pi*I*fd*t],{tao,-Infinity,Infinity},{fd,-100,100}];
Xwbhfm=Sqrt[Abs[a]]*Integrate[xhfm[t]*Conjugate[xhfm[a(t-tao)]],{tao,-Infinity,Infinity},{a,0.8,1.2}];
后来改了一下是这样的
f0=150;
m=400;
T=1;
f1=550;
k=(f0-f1)/(f0*f1*T);
fs=2500;
Ts=1/fs;
xlfm=Exp[I*2Pi*(f0*t+m*t^2*0.5)];
xlfm1=Exp[I*2Pi*(f0*(t-tao)+m*(t-tao)^2*0.5)];
xlfma=Exp[I*2Pi*(f0*a(t-tao)+m*a (t-tao)^2*0.5)];
xhfm=Exp[I*2Pi*((1/k)Log[1+k*f0*t])];
xhfm1=Exp[I*2Pi*((1/k)Log[1+k*f0*(t-tao)])];
xhfma=Exp[I*2Pi*((1/k)Log[1+k*f0*a(t-tao)])];
Xnblfm=Integrate[xlfm*Conjugate[xlfm1]*Exp[2Pi*I*fd*t],{tao,-0.1,0.1},{fd,-100,100}];
Xwblfm=Sqrt[Abs[a]]*Integrate[xlfm*Conjugate[xlfma],{tao,-0.1,0.1},{a,0.8,1.2}];
Xnhfm=Integrate[xhfm*Conjugate[xhfm1]*Exp[2Pi*I*fd*t],{tao,-0.1,0.1},{fd,-100,100}];
Xwbhfm=Sqrt[Abs[a]]*Integrate[xhfm*Conjugate[xhfma],{tao,-0.1,0.1},{a,0.8,1.2}];
错误信息见下图
2018年04月25日 02点04分 1
level 8
专业性太强没看懂是个啥
但是你如果积分不行
试试NIntegrate???
后面你的Abs[a]里面的a不是变量??和Double Integral的a没有关系??
2018年04月27日 23点04分 3
1