用matlab语言编通信模拟调频调制系统,包括积分微分等程序.doc

傅立叶变换子程序function M,m,df=fftseq(m,ts,df)fs=1/ts;if nargin=2 %nargin为输入参量的个数 n1=0;else n1=fs/df;endn2=length(m);n=2(max(nextpow2(n1),nextpow2(n2); % nextpow2(n)取n最接近的较大2次幂M=fft(m,n); % M为信号m的傅里叶变换，n为快速傅里叶变换的点数，%及为基n-FFT变换m=m,zeros(1,n-n2); %构建新的m信号df=fs/n; %重新定义频率分辨率主程序close allt0=8; %信号的持续时间,用来定义时间轴的取值范围ts=0.0001; %抽样间隔fs=1/ts; %抽样频率fc=10000; %载波频率t=-t0/2:ts:t0/2; %时间轴的取值kf=5; %压控振荡系数df=0.25; %所需的频率分辨率，用在求傅里叶变换时%它表示FFT的最小频率间隔m=2*pi*t; %调制信号%*jifen(1)=0; %求信号m(t)的积分for i=1:length(t)-1 jifen(i+1)=jifen(i)+m(i)*ts;end%*M,m,df1=fourier(m,ts,df); %对调制信号m(t)求傅里叶变换，需调用子函数fourierM=M/fs; %缩放，便于在频谱图上整体观察f=1:length(m)-fs/2; %时间轴上对应的频率向量u=1.414*cos(2*pi*fc*t+2*pi*kf*jifen); %调制后的信号U,u,df1=fourier(u,ts,df); %对调制后的信号u(t)求傅里叶变换U=U/fs; %缩放%*t1=0:ts:ts*(length(u)-1); %对已调信号进行解调z=hilbert(u); %希尔伯特变换yq =z.*exp(-j*2*pi*fc*t1);dem =diff(unwrap(angle(yq)*fs)/(2*pi*kf); %unwrap函数用于恢复相位并去掉2相位叠加的影响for i=1:(length(t)-1)/2 dem(i)=-dem(i);end%*figure(1)subplot(2,1,1)plot(t,m(1:length(t) %此时的m信号是重新构建的信号 %因为在对m求傅里叶变换时m=m,zeros(1,n-n2)xlabel(时间t)title(调制信号的时域波形)subplot(2,1,2)plot(t,u(1:length(t)axis(-1.5,1.5,-1.5,1.5);xlabel(时间t)title(已调信号的时域波形)figure(2)subplot(2,1,1)plot(f,abs(fftshift(M) %fftshift：将FFT中的DC分量移到频谱中心axis(58000,63000,0,2);xlabel(频率f)title(调制信号的频谱图)subplot(2,1,2)plot(f,abs(fftshift(U)axis(58000,63000,0,0.5);xlabel(频率f)title(已调信号的频谱图)figure(3)subplot(2,1,1)plot(t,m(1:length(t)xlabel(时间t)title(调制信号的时域波形)subplot(2,1,2)plo