1、n=0;t=n*T;plot(t,xn);t/sylabel(yn);axis(0,t(end),min(xn),1.2*max(xn);myplot.m(实验一、四需要)%(1)myplot;计算时域离散系统损耗函数并绘制曲线图。 function myplot(B,A)%B为系统函数分子多项式系数向量%A为系统函数分母多项式系数向量H,W=freqz(B,A,1000)m=abs(H);plot(W/pi,20*log10(m/max(m);grid on;omega/pi幅度(dB)axis(0,1,-80,5);title(损耗函数曲线mstem.m(实验一、三需要)function
2、mstem(Xk)%mstem(Xk)绘制频域采样序列向量Xk的幅频特性图M=length(Xk);k=0:M-1;wk=2*k/M;%产生M点DFT对应的采样点频率(关于pi归一化值)stem(wk,abs(Xk),box on;%绘制M点DFT的幅频特性图w/pi幅度axis(0,2,0,1.2*max(abs(Xk);mpplot.m(实验一需要)%(2)mpplot;计算时域离散系统损耗函数和相频特性函数,并绘制曲线图。function mpplot(B,A,Rs)%mpplot(B,A,Rs)%时域离散系统损耗函数和相频特性绘图%Rs为滤波器阻带最小衰减,省略则幅频曲线最小值取80d
3、Bif nargin3 ymin=-80;else ymin=-Rs-20;end;%确定幅频曲线纵坐标最小值H,W=H,W=freqz(B,A,1000)subplot(2,2,1);axis(0,1,ymin,5);subplot(2,2,3);plot(W/pi,p/pi);相位/pi(b)相频特性曲线mfftplot.m(实验一需要)function mfftplot(xn,N)%mfftplot(xn,N)计算序列向量xn的N点fft并绘制其幅频特性曲线Xk=fft(xn,N);%计算信号xn的频谱的N点采样%=以下为绘图部分=N-1;wk=2*k/N;m=abs(Xk);mm=ma
4、x(m);plot(wk,m/mm);axis(0,2,0,1.2);幅度特性曲线mstg.m(实验四需要)function st=mstg%产生信号序列向量st,并显示st的时域波形和频谱%st=mstg 返回三路调幅信号相加形成的混合信号,长度N=1600N=1600; %N为信号st的长度Fs=10000;T=1/Fs;Tp=N*T; %采样频率Fs=10kHz,Tp为采样时间t=0:T:(N-1)*T;f=k/Tp;fc1=Fs/10; %第1路调幅信号的载波频率fc1=1000Hzfm1=fc1/10; %第1路调幅信号的调制信号频率fm1=100Hzfc2=Fs/20; %第2路调
5、幅信号的载波频率fc2=500Hzfm2=fc2/10; %第2路调幅信号的调制信号频率fm2=50Hzfc3=Fs/40; %第3路调幅信号的载波频率fc3=250Hzfm3=fc3/10; %第3路调幅信号的调制信号频率fm3=25Hzxt1=cos(2*pi*fm1*t).*cos(2*pi*fc1*t); %产生第1路调幅信号xt2=cos(2*pi*fm2*t).*cos(2*pi*fc2*t); %产生第2路调幅信号xt3=cos(2*pi*fm3*t).*cos(2*pi*fc3*t); %产生第3路调幅信号st=xt1+xt2+xt3; %三路调幅信号相加fxt=fft(st,
6、N); %计算信号st的频谱%=以下为绘图部分,绘制st的时域波形和幅频特性曲线=subplot(3,1,1);plot(t,st);grid;s(t)axis(0,Tp/8,min(st),max(st);(a) s(t)的波形subplot(3,1,2);stem(f,abs(fxt)/max(abs(fxt),(b) s(t)的频谱axis(0,Fs/5,0,1.2);f/Hzxtg.m(实验五需要)function xt=xtg(N)%实验五信号x(t)产生,并显示信号的幅频特性曲线%xt=xtg(N) 产生一个长度为N,有加性高频噪声的单频调幅信号xt,采样频率Fs=1000Hz%载
7、波频率fc=Fs/10=100Hz,调制正弦波频率f0=fc/10=10Hz.Fs=1000;fc=Fs/10;f0=fc/10; %载波频率fc=Fs/10,单频调制信号频率为f0=Fc/10;mt=cos(2*pi*f0*t); %产生单频正弦波调制信号mt,频率为f0ct=cos(2*pi*fc*t); %产生载波正弦波信号ct,频率为fcxt=mt.*ct; %相乘产生单频调制信号xtnt=2*rand(1,N)-1; %产生随机噪声nt%设计高通滤波器hn,用于滤除噪声nt中的低频成分,生成高通噪声%=fp=150; fs=200;Rp=0.1;As=70; % 滤波器指标fb=fp
8、,fs;m=0,1; % 计算remezord函数所需参数f,m,devdev=10(-As/20),(10(Rp/20)-1)/(10(Rp/20)+1);n,fo,mo,W=remezord(fb,m,dev,Fs); % 确定remez函数所需参数hn=remez(n,fo,mo,W); % 调用remez函数进行设计,用于滤除噪声nt中的低频成分yt=filter(hn,1,10*nt); %滤除随机噪声中低频成分,生成高通噪声ytxt=xt+yt; %噪声加信号fst=fft(xt,N);plot(t,xt);x(t)axis(0,Tp/5,min(xt),max(xt);(a) 信
9、号加噪声波形plot(f,abs(fst)/max(abs(fst);(b) 信号加噪声的频谱axis(0,Fs/2,0,1.2);10.1 系统响应及系统稳定性close all;clear all;clc;%容1:调用filter解差分方程,由系统对u(n)的响应判断稳定性A=1,-0.9;B=0.05,0.05; %系统差分方程系数向量B和Ax1n=1 1 1 1 1 1 1 1 zeros(1,50); %产生信号x1(n)=R8(n)x2n=ones(1,128); %产生信号x2(n)=u(n)hn=impz(B,A,58); %求系统单位脉冲响应h(n)y=h(n);tstem(
10、hn,y); %调用函数tstem绘图(a) 系统单位脉冲响应h(n)box ony1n=filter(B,A,x1n); %求系统对x1(n)的响应y1(n)subplot(2,2,2);y1(n)tstem(y1n,y);(b) 系统对R8(n)的响应y1(n)y2n=filter(B,A,x2n); %求系统对x2(n)的响应y2(n)subplot(2,2,4);y2(n)tstem(y2n,y);(c) 系统对u(n)的响应y2(n)%容2:调用conv函数计算卷积x1n=1 1 1 1 1 1 1 1 ;h1n=ones(1,10) zeros(1,10);h2n=1 2.5 2.5 1 zeros(1,10);y21n=conv(h1n,x1n);y22n=conv(h2n,x1n);figure(2)h1(n)tstem(h1n,y);(d) 系统单位脉冲响应h1(n)y21(n)tstem(y21n,y);(e) h1(n)与R8(n)的卷积y21(n)h2(n)tstem(h2n,y);(f) 系统单位脉冲响应h2(n)y22(n)tstem(y22n,y);(g) h2(n)与R8(n)的卷积y22(n)%容3:谐振器分析un=ones(1,256); %产生信号u(n)255;xsin=sin(