MATLAB运行中出现 Undefined command/function 'fftseq'. 求大神纠错
基于MATLAB的FM调制t0=0.2;%设定初始值tz=0.0001;%设定时间步长fz=1/tz;%设定抽样频率t=[-t0:tz:t0];%产生时间向量fc=100...
基于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) '); 展开
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) '); 展开
2015-04-18 · 知道合伙人数码行家
关注
展开全部
这个谁 会去看这个呢
如有疑问请追问,
如满意记得采纳,
如有其他问题也可点我名字向我求助,
答题不易,
请采纳最快回答的正确答案!! 谢谢小主!!
如有疑问请追问,
如满意记得采纳,
如有其他问题也可点我名字向我求助,
答题不易,
请采纳最快回答的正确答案!! 谢谢小主!!
推荐律师服务:
若未解决您的问题,请您详细描述您的问题,通过百度律临进行免费专业咨询