数字信号处理
「已注销」
这个作者很懒,什么都没留下…
展开
专栏收录文章
- 默认排序
- 最新发布
- 最早发布
- 最多阅读
- 最少阅读
-
数字信号处理基本指令(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 · 526 阅读 · 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 · 1109 阅读 · 0 评论 -
数字信号基本指令(4)实验四 离散信号FFT的实现实验
caiyangpinlv=200; %fft-1 figure(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 · 368 阅读 · 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 · 1779 阅读 · 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 · 253 阅读 · 0 评论 -
数字信号处理基本指令(1)
实验一 离散时间信号产生,图形显示及文件读存实验 %% 数字信号处理实验1 内容1 n1=-1; n2=8; n0=0; n3=5; if((n0<n1)||(n0>n2)||(n1>n2)) error('参数必须满足n1<=n0<=n2') end n=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 · 398 阅读 · 0 评论
分享