实验五数字信号代码.doc

产生xt信号代码 function xt=xtg %滤波器指标：通带截止频率fp=120hz,阻带截止频率fs=150hz,换算成数字频率为wp=2*pi*fp/fs=0.24*pi %通带最大衰减0.1dB，阻带截止频率ws=0.3*pi，阻带最小衰减60dB %实验五信号x(t)产生,并显示信号的幅频特性曲线 %xt=xtg(N) 产生一个长度为N,有加性高频噪声的单频调幅信号xt,采样频率Fs=1000Hz %载波频率fc=Fs/10=100Hz,调制正弦波频率f0=fc/10=10Hz. N=1000;Fs=1000;T=1/Fs;Tp=N*T;%采样时间为Tp=1s t=0:T:(N-1)*T; fc=Fs/10;f0=fc/10; %载波频率fc=Fs/10，单频调制信号频率为f0=Fc/10; mt=cos(2*pi*f0*t); %产生单频正弦波调制信号mt，频率为f0 ct=cos(2*pi*fc*t); %产生载波正弦波信号ct，频率为fc xt=mt.*ct; %相乘产生单频调制信号xt,抑制载波信号 nt=2*rand(1,N)-1; %产生随机噪声nt %=======设计高通滤波器hn,用于滤除噪声nt中的低频成分,生成高通噪声======= fp=150; fs=200;Rp=0.1;As=70; % 滤波器指标 fb=[fp,fs];m=[0,1]; % 计算remezord函数所需参数f,m,dev dev=[10^(-As/20),(10^(Rp/20)-1)/(10^(Rp/20)+1)];%允许的波纹幅度，求解通带波纹系数，和阻带波纹系数 [n,fo,mo,W]=remezord(fb,m,dev,Fs); % 确定remez函数所需参数，最低阶数m，误差加权值w，归一化边界频率f hn=remez(n,fo,mo,W); % 调用remez函数进行设计,用于滤除噪声nt中的低频成分 yt=filter(hn,1,10*nt); %滤除随机噪声中低频成分，生成高通噪声yt %================================================================ xt=xt+yt; %噪声加信号 fst=fft(xt,N);k=0:N-1;f=k/Tp; subplot(3,1,1);plot(t,xt);grid;xlabel(t/s);ylabel(x(t)); axis([0,Tp/5,min(xt),max(xt)]);title((a) 信号加噪声波形) subplot(3,1,2);plot(f,abs(fst)/max(abs(fst)));grid;title((b) 信号加噪声的频谱) axis([0,Fs/2,0,1.2]);xlabel(f/Hz);ylabel(幅度) 用窗函数法或等波纹法设计滤波器代码 clear all;close all Fs=1000;T=1/Fs; %采样频率 N=1000;xt=xtg ; %调用信号产生函数mstg产生由三路抑制载波调幅信号相加构成的复合信号st %低通滤波器设计与实现 fp=120;fs=150; wp=0.23*pi;ws=0.3*pi;rp=0.1;rs=60; B=ws-wp; %计算过度带带宽 M=ceil(11*pi/B)-1; %确保M为奇数 wc=(wp+ws)/2/pi;%计算理想高通滤波器通带截止频率（关于pi的归一化） hn=fir1(M,wc,blackman(M+1));%布莱科曼窗，返回6dB截止频率wc和M阶 ，hn为滤波器的单位脉冲响应 %以下为用窗函数法设计法的绘图部分（滤波器损耗函数，滤波器输出信号波形） Hw=abs(fft(hn,1024)); % 求设计的滤波器频率特性 ywt=fftfilt(hn,xt,N); %调用函数fftfilt对xt滤波，ywt为滤波后的信号 %fftfit ,hn是系统单位脉冲响应，xt是输入序列，N是输入序列的长度，ywt是hn和xt卷积的结果 %以下为用窗函数法设计法的绘图部分（滤波器损耗函数，滤波 figure(2); subplot(3,1,1); myplot(hn,1);%有限长hn的系数A=1 title((a)低通滤波器幅频特性); subplot(3,1,2); yt=y_w(t); tplot(ywt,T,yt); title(滤除噪声后的信号波形) %画出ywt信号的频谱。 N=1000 Fs=1000;T=1/Fs;Tp=N*T; %采样频率Fs=10kHz，Tp为采样时间 t=0:T:(