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
就是定义了两个信号和两个多变量函数,要求计算两信号对应的多变量函数在一定区间上积分所得的关于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}];
错误信息见下图





