matlab中fftseq怎么用,MATLAB运行中出现 Undefined command/function 'fftseq'. 求大神纠错...

MATLAB运行中出现 Undefined command/function 'fftseq'. 求大神纠错0

基于MATLAB的FM调制

t0=0.2;%设定初始值

tz=0.0001; %设定时间步长

fz=1/tz; %设定抽样频率

t=[-t0:tz:t0]; %产生时间向量

fc=1000; %载波信号频率

kf=100; %设定调频指数

df=0.05; %设定分辨率

m_fun=cos(400*pi*t);

int_m(1)=0; %对m_fun积分

for i=1:length(t)-1

int_m(i+1)=int_m(i)+m_fun(i)*tz;

end

x=cos(2*pi*fc*t+2*pi*kf*int_m); %调制信号

[M,m_fun,df1]=fftseq(m_fun,tz,df); %对原始信号快速傅里叶变换

M=M/fz;

f=[0:df1:df1*(length(m_fun)-1)]-fz/2;

[X,x,df1]=fftseq(x,tz,df); %对已调信号快速傅里叶变换

X=X/fz;

figure(1); %生成原始信号的时域波形

plot(t,m_fun(1:length(t)), 'linewidth',3);

axis([-0.01 0.01 -1.5 1.5]);

title('原始信号的时域波形')；

xlable('时间');

legend('m(t) ');

figure(2); %生成原始信号的频域图形

plot(f,abs(fftshift(M)), 'linewidth',3);

axis([-400400 -0.01 0.1]);

title('原始信号的频域波形')；

xlable('频率');

legend('M(f) ');

figure(3); %生成已调信号的时域波形

plot(t,x(1:length(t)), 'linewidth',3);

axis([-0.015 0.015 -1.5 1.5]);

title('已调信号的时域波形')；

xlable('时间');

legend('x(t) ');

figure(4); %生成已调信号的频域波形

plot(f,abs(fftshift(X)), 'linewidth',3);

axis([-15001500 0 0.1]);

title('已调信号的频域波形')；

xlable('频率');

legend('X(f) ');