通信原理课程设计(10电信)0612.ppt

%显示模拟调制的波形方法AM，文件S_AM.m %Signal dt=0.001; %时间采样间隔 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]=F2T (t,s_am); %调制信号频谱 PSD=(abs(Xf).^2)/T; %调制信号功率谱密度 figure(1) subplot(211); plot(t,s_am);hold on; %画出AM信号波形 plot(t,A+mt,r--); %标示AM的包络 title(AM调制信号及其包络); xlabel(t); subplot(212); %画出功率谱图形 plot(f,PSD); axis([-2*fc 2*fc 0 1.5*max(PSD)]); title(AM信号功率谱); xlabel(f); 同理，在命令窗口中键入x的二进制代码和函数名，就可以得到所对应的单极性归零码输出，如输入以下指令，将出现图示结果。 x=[1 0 1 1 0 0 1 0]; srz(x) 双极性归零码 双极性归零码的MATLAB实现同单极性也基本一样，只需将srz.m中的判断得到0信息后的语句 for j=1:t0 y((i-1)*t0+j)=0; 改为for j=1:t0/2 y((2*i-2)*t0/2+j)=-1; y((2*i-1)*t0/2+j)=0; 即可，其波形图如图示。 * * 通信原理课程设计 2012年12月 通信MATLAB仿真基础 信号 syms a b c solve(a*x^2+b*x+c=0,x) ans = 1/2/a*(-b+(b^2-4*a*c)^(1/2)) 1/2/a*(-b-(b^2-4*a*c)^(1/2)) %%% waveform and spectrum of sin signal close all K = 10; f1 = 1; f2 = 2; N = 2^k; dt = 0.01; %ms df = 1/(N*dt); % KHz T = N*dt; % 截短时间 Bs = N*df/2; % 系统带宽 f = [-Bs:df:Bs-df]; %频域横坐标 t = [-T/2:dt:T/2-dt]; %时域横坐标 s=2*sin(2*pi*f1*t)+3*sin(2*pi*f2*t); [f,S] = T2F(t,s); % S是s的傅氏变换 [t,a] = F2T(f,S); % a是S的傅氏反变换 a=real(a); as=abs(S); f0=max(f1,f2); figure subplot(2,1,1) %输出的频谱 plot(f,as,b); grid axis([-2*f0,+2*f0,min(as),max(as)]) xlabel(f (KHz)); ylabel(|S(f)| (V/KHz)) subplot(2,1,2) plot(t,a,black) %输出信号波形画图 grid axis([-2/f0,+2/f0,-5,5]) xlabel(t(ms)); ylabel(a(t)(V)) 傅里叶变换 function [f,sf] = T2F(t,st) % calculate a signals Fourier