现代数字信号处理及应用仿真题答案.docxVIP

下载本文档

205
0
约5.72千字
约 17页
2019-07-03 发布于浙江
举报
版权申诉

现代数字信号处理及应用仿真题答案.docx

1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
7、VIP文档为合作方或网友上传，每下载1次，网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档

仿真作业姓名：李亮学号：S130101083 4.17程序 clc; clear; for i=1:500 sigma_v1=0.27; b(1)=-0.8458; b(2)=0.9458; a(1)=-(b(1)+b(2)); a(2)=b(1)*b(2); datlen=500; rand(state,sum(100*clock)); s=sqrt(sigma_v1)*randn(datlen,1); x=filter(1,[1,a],s); %% sigma_v2=0.1; u=x+sqrt(sigma_v2)*randn(datlen,1); d=filter(1,[1,-b(1)],s); %% w0=[1;0]; w=w0; M=length(w0); N=length(u); mu=0.005; for n=M:N ui=u(n:-1:n-M+1); y(n)=w*ui; e(n)=d(n)-y(n); w=w+mu.*conj(e(n)).*ui; w1(n)=w(1); w2(n)=w(2); ee(:,i)=mean(e.^2,2); end end ep=mean(ee); plot（ep）; xlabel(迭代次数);ylabel(MSE);title(学习曲线); plot(w1); hold; plot(w2); 仿真结果：步长0.015仿真结果步长0.025仿真结果步长0.005仿真结果 4.18 程序 data_len = 512; %样本序列的长度 trials = 100; %随机试验的次数 A=zeros(data_len,2);EA=zeros(data_len,1); B=zeros(data_len,2);EB=zeros(data_len,1); for m = 1: trials a1 = -0.975; a2 = 0.95; sigma_v_2 =0.0731; v = sqrt(sigma_v_2) * randn(data_len, 1, trials);%产生v(n) u0 = [0 0]; num = 1; den = [1 a1 a2]; Zi = filtic(num, den, u0); %滤波器的初始条件 u = filter(num, den, v, Zi); %产生样本序列u(n) %（2）用LMS滤波器来估计w1和w2 mu1 = 0.05; mu2 = 0.005; w1 = zeros(2, data_len); w2 = zeros(2, data_len); e1 = zeros(data_len, 1); e2 = zeros(data_len, 1); d1 = zeros(data_len, 1); d2 = zeros(data_len, 1); %LMS迭代过程 for n =3 :data_len - 1 w1( :, n+1) = w1( :, n) + mu1 * u(n-1 : -1: n-2, : , m) * conj(e1(n)); w2( :, n+1) = w2( :, n) + mu2 * u(n-1 : -1: n-2, : , m) * conj(e2(n)); d1(n+1) = w1( : , n+1) * u(n: -1: n-1, :, m); d2(n+1) = w2( : , n+1) * u(n: -1: n-1, :, m); e1(n+1) = u(n+1, : ,m) - d1(n+1); e2(n+1) = u(n+1, : ,m) - d2(n+1); end A = A + conj(w1); EA = EA +e1.^2; B = B + conj(w2); EB = EB + e2.^2; end %剩余均方误差和失调参数 wopt=zeros(2,trials); Jmin=zeros(1,trials); sum_eig=zeros(trials,1); for m=1:trials; rm=xcorr(u(:,:,m),biased); R=[rm(512),rm(513);rm(511),rm(512)]; p=[rm(511);rm(510)]; wopt(:,m)=R\p; [v,d]=eig(R); Jmin(m)=rm(512)-p*wopt(:,m); sum_eig(m)=d(1,1)+d(2,2); end sJmin=sum(Jmin)/trials; e1_100t