数字信号处理
「已注销」
这个作者很懒,什么都没留下…
展开
-
数字信号处理基本指令(6)实验六 IIR滤波器的设计及应用实验
实验题目以下程序为主程序,本实验还需要同时写mstg函数,需同时有mstg函数建立在同一个文件夹下。st=mstg;Fs=10000;fp=280;fs=450;wp=2*fp/Fs;ws=2*fs/Fs;rp=0.1;rs=60;[N,wp0]=ellipord(wp,ws,rp,rs);[B,A]=ellip(N,rp,rs,wp0);y1t=filter(B,A,st);figure(2);[H1,w]=freqz(B,A,1000);m=abs(H1);plot(w/pi原创 2022-02-06 12:35:24 · 420 阅读 · 0 评论 -
数字信号处理基本指令(5)实验五 离散信号IFFT的实现实验
t=-1:0.001:1;x=0.7*sin(2*pi*20*t)+3*sin(80*3*pi*t);figure(1);plot(t,x)xlabel('时间/s');ylabel('函数值');title('原始函数 班级 学号 姓名');fs=200;N_data=64;N=128;n=0:N_data-1;t=n/fs;y=fft(x,N);mag=abs(y);f=(0:N-1)*fs/N;figure(2);plot(f(1:N/2),mag(1:N/2)*2/.原创 2022-02-06 12:25:30 · 1035 阅读 · 0 评论 -
数字信号基本指令(4)实验四 离散信号FFT的实现实验
caiyangpinlv=200;%fft-1figure(1);N_data=64;N=64;n=0:N_data-1;t=n/caiyangpinlv;x=2*sin(50*2*pi*t)+10*sin(20*2*pi*t);y=fft(x,N);mag=abs(y);f=(0:N-1)*caiyangpinlv/N;plot(f(1:N/2),mag(1:N/2)*2/N);xlabel('F/Hz');ylabel('振幅|X(z)|');title('班级 学号 姓名.原创 2022-02-06 12:18:07 · 319 阅读 · 0 评论 -
数字信号处理基本指令(3)实验三 吉布斯现象实验
(1)生成矩形波信号t=-6:0.01:6;x=0.5*square(pi*t);plot(t,x);title('班级 学号 姓名 矩形波信号');xlabel('时间t/s');ylabel('幅值');axis([-6 6 -0.6 0.8]);上图为矩形波信号(2)求矩形波信号的傅里叶展开式;(3)用plot函数绘制傅里叶级数展开式,t=-2:0.01:2;T=[1 3 5 7 9 15 25 51 101 151];for k=1:length(T)a=0;for.原创 2022-02-06 12:02:45 · 1636 阅读 · 0 评论 -
数字信号处理基本指令(2)
实验二 离散信号的卷积和实验np1=8;N1=np1;ns=0;nf1=8;n1=ns:nf1;np=0;%单位阶跃信号x1=stepseq(np,ns,nf1)-stepseq(np1,ns,nf1);hn=((0.8).^n1).*x1;%激励信号np2=4;N2=np2;nf2=4;n2=ns:nf2;x2=stepseq(np,ns,nf2)-stepseq(np2,ns,nf2);figure(1);subplot(211);stem(n1,hn);xla原创 2022-02-06 11:44:55 · 217 阅读 · 0 评论 -
数字信号处理基本指令(1)
实验一 离散时间信号产生,图形显示及文件读存实验%% 数字信号处理实验1 内容1n1=-1;n2=8;n0=0;n3=5;if((n0<n1)||(n0>n2)||(n1>n2)) error('参数必须满足n1<=n0<=n2')endn=n1:n2;x=[(n-n0)==0];y=[(n-n3)==0];figure(1);subplot(211);stem(n,x);xlabel('N');ylabel('theta(n)');t原创 2022-02-06 11:34:34 · 349 阅读 · 0 评论