数字信号处理实验一.docx
- 文档编号:3470769
- 上传时间:2022-11-23
- 格式:DOCX
- 页数:15
- 大小:123.84KB
数字信号处理实验一.docx
《数字信号处理实验一.docx》由会员分享,可在线阅读,更多相关《数字信号处理实验一.docx(15页珍藏版)》请在冰豆网上搜索。
数字信号处理实验一
1.R=50;
m=0:
R-1;
x1=0;
s=2*m.*(0.9.^m);
d=rand(R,1)-0.5;
forn=1:
R;
d=rand(R,1)-0.5;
x=s+d';
x1=x1+x;
end
x1=x1/R;
subplot(2,2,1);stem(m,s);
xlabel('Timeindexn');ylabel('Amplitude');title('original');
subplot(2,2,2);stem(m,d);
xlabel('Timeindexn');ylabel('Amplitude');title('Noise');
subplot(2,2,3);stem(m,x);
xlabel('Timeindexn');ylabel('Amplitude');title('Noisecorrupted');
subplot(2,2,4);stem(m,x1);
xlabel('Timeindexn');ylabel('Amplitude');title('Ensemblemaverage');
2.
(1)a=input('Typeinrealexponent=');
b=input('Typeinimaginaryexponent=');
c=a+b*i;
K=input('Typeinthegainconstant=');
N=input('Typeinlengthofsequence=');
n=1:
N;
x=K*exp(c*n);
stem(n,real(x));
xlabel('Timeindexn');ylabel('Amplitude');
title('Realpart');
disp('PRESSRETURNforimaginarypart');
pause
stem(n,imag(x));
xlabel('Timeindexn');ylabel('Amplitude')
title('Imaginarypart');
(2)a=input('Typeinthebasenumber=');
b=input('Typeinthecoefficient=');
K=input('Typeinthegainconstant=');
N=input('Typeinlengthofsequence=');
n=1:
N;
x=K.*b.*(a.^n);
stem(n,x);
xlabel('Timeindexn');ylabel('Amplitude');
title('Realexponential');
A.a=1.2
B.a=0.9
(3)w=input('Typeinthew=');
K=input('Typeinthegainconstant=');
N=input('Typeinlengthofsequence=');
n=1:
N;
x=K*1.5*cos(w*n);
stem(n,x);
xlabel('Timeindexn');ylabel('Amplitude');
title('Sinusoidalsequences');
(a)w=0
(b)w=0.1*pi
(c)w=0.2*pi
(d)w=0.8*pi
(e)w=0.9*pi
(f)w=pi
(g)w=1.1*pi
(h)w=1.2*pi
3.1
functionprogram2_3
R=50;
d=rand(R,1)-0.5;
m=0:
1:
R-1;
s=2*m.*(0.9.^m);
x=s+d';
plot(m,d,'r-',m,s,'b--',m,x,'g:
')
xlabel('Timeindexn');ylabel('Amplitude');
legend('d[n]','s[n]','x[n]');
3.2
functionprogram2_32
R=50;
d=rand(R,1)-0.5;
m=0:
1:
R-1;
s=2*m.*(0.9.^m);
x=s+d';
M=input('Numberofinputsamples=');
b=ones(M,1)/M;
y=filter(b,1,x);
plot(m,s,'r-',m,y,'b--')
xlabel('Timeindexn');ylabel('Amplitude');
legend('s[n]','y[n]');
3.3
functionprogram2_33
R=50;
d=rand(R,1)-0.5;
m=0:
1:
R-1;
s=2*m.*(0.9.^m);
x=s+d';
M=input('Numberofinputsamples=');
fori=M:
R;
sum=0;
forl=1:
M;
sum=sum+x(i-l+1);
end
y(i)=sum/M;
end
plot(m,s,'r-',m,y,'b--')
xlabel('Timeindexn');ylabel('Amplitude');
legend('s[n]','y[n]');
3.4
A.Numberofinputsamples=5
B.Numberofinputsamples=7
C.Numberofinputsamples=9
4.1
functionprogram2_41
R=50;a=rand(1,R)-0.4;
b=round(a);
m=0:
R-1;
s=2*m.*(0.9.^m);
x=s+b;
M=input('MedianFilterlength=');
y=medfilt1(x,M);
subplot(2,1,1)
stem(m,x);axis([050-18]);
xlabel('n');ylabel('Amplitude');
title('TmpulseNoiseCorruptedSignal');
subplot(2,1,2)
stem(m,y);
xlabel('n');ylabel('Amplitude');
title('OutputofMedianFilter');
MedianFilterlength=3
4.3
R=50;
a=rand(1,R)-0.4;
b=round(a);
m=0:
R-1;
s=2*m.*(0.9.^m);
x=s+b;
M=input('Numberofinputsamples=');
x2=[zeros(1,(M-1)/2),x,zeros(1,(M-1)/2)];
fori=1:
R;
c=sort(x2(i:
i+M-1));
ifmod(M,2)==1;y(i)=c((M+1)/2);
elsey(i)=(c(M/2)+c(M/2+1))/2;
end
end
subplot(2,1,1)
stem(m,x);axis([050-18]);
xlabel('n');ylabel('Amplitude');
title('TmpulseNoiseCorruptedSignal');
subplot(2,1,2)
stem(m,y);
xlabel('n');ylabel('Amplitude');
title('OutputofMedianFilter');
Numberofinputsamples=3
5.
functionprogram2_5
a=input('Typeinthefirstsequence=');
b=input('Typeinthesecondsequence=');
c=conv(a,b);
M=length(c)-1;
n=0:
1:
M;
disp('outputsequence=');disp(c);
stem(n,c)
xlabel('Timeindexn');ylabel('Amplotude');
A.
Typeinthefirstsequence=[-201-13]
Typeinthesecondsequence=[120-1]
outputsequence=
-2-413151-3
B.Typeinthefirstsequence=[-201-13]
Typeinthesecondsequence=[10000]
outputsequence=
-201-130000
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 数字信号 处理 实验