AM的模拟调制程序要出现AM调制信号及其包络和AM信号功率谱

jiayue

程序如下：
%显示模拟调制的波形方法AM，文件S_AM.m
%Signal
dt=0.0001;
fmax=1;
fc=10;
T=5;
N=T/dt;
t=[0:N-1]*dt;
mt=sqrt(2)*cos(2*pi*fmax*t);
%AM modulation
A=2;
s_am=(A+mt).*cos(2*pi*fc*t);
%Power Spectrum Density
[f,Xf]=FFT_SHIFT(t,s_am);
PSD=(abs(Xf).^2)/T;
figure(1)
subplot(211);
plot(t,s_am);hold on;
plot(t,A+mt,'r--');
title('AM调制信号及其包络');
xlabel('t');
subplot(212);
plot(f,PSD);
axis([-2*fc 2*fc 0 1.5*max(PSD)]);
title('AM信号功率谱');
xlabel('f');
出现的错误是：
??? Undefined function or method 'FFT_SHIFT' for input arguments of type
'double'.
Error in ==> AM at 14
[f,Xf]=FFT_SHIFT(t,s_am);
哪位好心人帮我解决下，问题出在哪要怎么修改，谢谢了！

jiayue

程序如下：
%显示模拟调制的波形方法AM，文件S_AM.m
%Signal
dt=0.0001;
fmax=1;
fc=10;
T=5;
N=T/dt;
t=[0:N-1]*dt;
mt=sqrt(2)*cos(2*pi*fmax*t);
%AM modulation
A=2;
s_am=(A+mt).*cos(2*pi*fc*t);
%Power Spectrum Density
double [f,Xf]=FFT_SHIFT(t,s_am);
double PSD=(abs(Xf).^2)/T;
figure(1)
subplot(211);
plot(t,s_am);hold on;
plot(t,A+mt,'r--');
title('AM调制信号及其包络');
xlabel('t');
subplot(212);
plot(f,PSD);
axis([-2*fc 2*fc 0 1.5*max(PSD)]);
title('AM信号功率谱');
xlabel('f');
出现的错误：
??? Undefined function or variable 'f'.
Error in ==> AM at 23
plot(f,PSD);
这是第二种情况，我修改了下，加了两个double ，但又出现了错误，哪位好心人可以帮我修改下，谢谢！

kongfeng0824

这是关于气象的吗？你可以详细说明一下你的研究目的吗？

jiayue

kongfeng0824 发表于 2013-4-21 05:37

                               
登录/注册后可看大图

这是关于气象的吗？你可以详细说明一下你的研究目的吗？

不是关于气象的，是关于matlab的，基于matlab的模拟调制系统仿真

kongfeng0824

jiayue 发表于 2013-4-22 13:33

                               
登录/注册后可看大图

不是关于气象的，是关于matlab的，基于matlab的模拟调制系统仿真

这个还真的不知道。不好意思。