资源描述：
《数字信号处理matlab编程作业》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、数字信号处理MATLAB编程作业姓名：白焱学号：2012001020006M2.2:代码：（以其中一组实例，其他组相同，其他类似情况不在说明）n1=-10:10;xn=5*cos(1.5*pi*n+0.75*pi)+4*cos(0.6*pi*n)-sin(0.5*pi*n);stem(xn);xlabel('n');ylabel('xe[n]');结果：M2.9代码：clear;clc;x=input('Typeinthereferencesequence=');y=input('Typeinthesecondsequence=');%Computethecorrelationsequence
2、n1=length(y)-1;n2=length(x)-1;r=conv(x,fliplr(y));k=(-n1):n2';stem(k,r);xlabel('Lagindex');ylabel('Amplitude');v=axis;axis([-n1n2v(3:end)]);结果：X[n]自相关序列y[n]自相关序列w[n]自相关序列X[n]与y[n]互相关序列X[n]与w[n]互相关序列M3.1代码：rr=input('valueofr=');xita=input('valueofxita=');fem(1)=0;fem(2)=-2*rr*cos(xita);fem(3)=rr*rr;n
3、um=[1];den=fem;%Computethefrequencyresponsew=0:pi/(k-1):pi;h=freqz(num,den,w);%Plotthefrequencyresponsesubplot(2,2,1)plot(w/pi,real(h));gridtitle('Realpart')xlabel('omega/pi');ylabel('Amplitude')subplot(2,2,2)plot(w/pi,imag(h));gridtitle('Imaginarypart')xlabel('omega/pi');ylabel('Amplitude')subp
4、lot(2,2,3)plot(w/pi,abs(h));gridtitle('MagnitudeSpectrum')xlabel('omega/pi');ylabel('Magnitude')subplot(2,2,4)plot(w/pi,angle(h));gridtitle('PhaseSpectrum')xlabel('omega/pi');ylabel('Phase,radians')结果：（r=0.8,θ=0.35）M4.1:代码：%Program4_2%ProgramtoDesignButterworthLowpassFilter%%Typeinthefilterorder
5、andpassbandedgefrequencyN=input('Typeinfilterorder=');Wn=input('3-dBcutoffangularfrequency=');%Determinethetransferfunction[num,den]=butter(N,Wn,'s');%Computeandplotthefrequencyresponseomega=[0:200:12000*pi];h=freqs(num,den,omega);plot(omega/(2*pi),20*log10(abs(h)));xlabel('Frequency,Hz');ylabel('Ga
6、in,dB');结果：M4.7：代码：[N,Wn]=cheb1ord(0.3157894,1,0.5,30,'s');[B,A]=cheby1(N,0.5,Wn,'s');[num,den]=lp2bs(B,A,2*pi*sqrt(700)*10^6,2*pi*15e6);figure(1)omega=[0:0.01:10];h=freqs(B,A,omega);gain=20*log10(abs(h));plot(omega,gain);grid;axis([04-705]);xlabel('Omega');ylabel('Gain,dB');title('AnalogLowpassFil
7、ter');figure(2)omega=[0:10000:160e6*pi];h=freqs(num,den,omega);gain=20*log10(abs(h));plot(omega/(2*pi),gain);grid;axis([080e6-705]);xlabel('FrequencyinHz');ylabel('Gain,dB');title('AnalogBandstopFilte
