第八讲__MATLAB在通信网中的应用.pptVIP

第八讲Matlab在通信网中的应用

信号的时域分析单边衰减指数信号t=0:0.01:10;A=1;A=-0.4;ft=A*exp(a*t);plot(t,ft)xlabel(‘t’)ylabel(‘f(t)’)

正弦信号t=0:0.001:8;A=1;w0=2*pi;phi=pi/6;ft=A*sin(w0*t+phi);plot(t,ft)xlabel(‘t’)ylabel(‘f(t)’)信号的时域分析

抽样函数信号t=-3*pi:pi/100:3*pi;ft=sinc(t/pi);plot(t,ft)xlabel(‘t’)ylabel(‘Sa(t)’)信号的时域分析

矩形脉冲信号t=0:0.001:4T=1;ft=rectpuls(t-2*T,T);plot(t,ft);axis([0,4,-0.5,1.5])xlabel(‘t’)ylabel(‘rectpuls(t)’)信号的时域分析

产生单位脉冲序列ks=-4;ke=4;n=2;k=[ks:ke];x=[(k-n)==0];stem(k,x);xlabel(‘k’);信号的时域分析

正弦序列k=0:39;fk=sin(pi/6*k);stem(k,fk)xlabel(‘k’)ylabel(‘sinusoidalsignal’)信号的时域分析

系统的时域分析连续时间系统的冲激响应ts=0;te=5;dt=0.01;sys=tf([10],[12100]);t=ts:dt:te;y=impulse(sys,t);plot(t,y);xlabel(‘Time(sec)’)ylabel(‘h(t)’)

离散系统的单位脉冲响应k=0:10;a=[132];b=[1];h=impz(b,a,k);subplot(2,1,1)stem(k,h)title(‘approximationofh(k)’);hk=-(-1).^k+2*(-2).^k;subplot(2,1,2)stem(k,hk)title(‘h(k)intheory’);系统的时域分析

信号的频域分析求周期序列的DFSN＝32；M=4;%定义周期方波序列的参数f=[ones(1,M+1)zeros(1,N-2*M-1)ones(1,M)];%产生序列F＝fft(f);%计算DFS系数m=0:N-1;subplot(2,1,1);stem(m,real(F));title(‘realofF[m]’);xlabel(‘m’);subplot(2,1,2);stem(m,imag(F));title(‘imaginaryofF[m]’);xlabel(‘m’);

连续周期信号的傅立叶级数的计算N=8;%计算n=-N到－1的傅立叶系数n1=-N:-1;c1=-4*j*sin(n1*pi/2)/pi^2./n1.^2;%计算n=0时的傅立叶系数c0=0;%计算n=1到N的傅立叶系数n2=1:N;c2=-4*j*sin(n2*pi/2)/pi^2./n2.^2;cn=[c1c0c2];n=-N:N;subplot(2,1,1);stem(n,abs(cn));ylabel(‘AmplitudeofCn’);subplot(2,1,2);stem(n,angle(cn));ylabel(‘PhaseofCn’);xlabel(‘\omega/\omega0’);信号的频域分析

系统的频域分析输入为余弦的连续系统响应RC=0.04;t=linspace(-2,2,1024);w1=5;w2=100;H1=j*w1/(j*w1+1/RC);H2=j*w2/(j*w2+1/RC);f=cos(5*t);y=abs(H1)*cos(w1*t+angle(H1))+abs(H2)*cos(w2*t+angle(H2));subplot(2,1,1);plot(t,f);ylabel(‘f(t)’);xlabel(‘t(sec)’)subplot(2,1,2);ylabel(‘y(t)’);xlabel(‘t(sec)’);

产生调制信号〔PM〕Fm=10;Fc=100;Fs=1000;N=1000;k=0:N-1;t=k/Fs;x=sin(2.0*pi*Fm*t);Xf=abs(fft(x,N));y2=modulate(x,Fc,Fs,’pm’);subplot(2,1,1);plot(t(1:200),y2(1:200));xlabel(‘Time(s)’);axis([0,0.2,-1,1]);title(‘ModulatedSignal(PM)’);Yf2=abs(fft(y2,N));subplot(2,1,2);stem(Yf2(1:200));xlabel(‘Frequency(