fm.rar_FM_FM matlab_matlab 棰戠巼璋冨埗_角度调制

clear clc echo on t0=0.15; %信号持续时间 tz=0.0005; %采样时间间隔 fc=200; %载波频率 kf=50; %调制系数 fz=1/tz; t=[0:tz:t0]; %定义时间序列 df=0.25; %定义信号序列 m=[ones(1,t0/(3*tz)),-2*ones(1,t0/(3*tz)),zeros(1,t0/(3*tz)+1)]; %100个1，100个-2，101个0 int_m(1)=0; %next int_m(1:301) is generated,int_m(1)=0;int_m(j)=sum[m(i),i=1,...,j-1]; echo off for i=1:length(t)-1 int_m(i+1)=int_m(i)+m(i)*tz; end echo on [M,m,df1]=fftseq(m,tz,df); M=M/fz; f=[0:df1:df1*(length(m)-1)]-fz/2; u=cos(2*pi*fc*t+2*pi*kf*int_m); %调制信号调制在载波上 [U,u,df1]=fftseq(u,tz,df); U=U/fz; pause; clf subplot(2,2,1) plot(t,m(1:length(t))) %作出信号波形 axis([0 0.15 -2.1 2.1]); xlabel('时间') title('未调信号'); pause subplot(2,2,2) plot(t,u(1:length(t))) %作出调频信号波形 axis([0 0.15 -2.1 2.1]); xlabel('时间') title('调频信号'); pause subplot(2,2,3) plot(f,abs(fftshift(M))) %作出原信号频谱 xlabel('频率') title('原信号频谱'); pause subplot(2,2,4) plot(f,abs(fftshift(U))) %作出调频信号频谱 xlabel('频率') title('调频信号频谱');