MATLAB上机实践报告

相关主题

1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。

clear all;

load('Hsoa2ib2.mat')

Hsoa2ib2=(a4+j*a5);f=a0;clear a0a1a2a3a4a5a6a7a8; load('Hsoa2ib1.mat')

Hsoa2ib1=(a4+j*a5);f=a0;clear a0a1a2a3a4a5a6a7a8; load('Hsoa1ib2.mat')

Hsoa1ib2=(a4+j*a5);f=a0;clear a0a1a2a3a4a5a6a7a8; load('Hsoa1ib1.mat')

Hsoa1ib1=(a4+j*a5);f=a0;clear a0a1a2a3a4a5a6a7a8; load('Hscb2cb2.mat')

Hscb2cb2=(a4+j*a5);f=a0;clear a0a1a2a3a4a5a6a7a8; load('Hscb2cb1.mat')

Hscb2cb1=(a4+j*a5);f=a0;clear a0a1a2a3a4a5a6a7a8; load('Hscb1cb2.mat')

Hscb1cb2=(a4+j*a5);f=a0;clear a0a1a2a3a4a5a6a7a8; load('Hscb1cb1.mat')

Hscb1cb1=(a4+j*a5);f=a0;clear a0a1a2a3a4a5a6a7a8; load('Hsca2cb2.mat')

Hsca2cb2=(a4+j*a5);f=a0;clear a0a1a2a3a4a5a6a7a8; load('Hsca2cb1.mat')

Hsca2cb1=(a4+j*a5);f=a0;clear a0a1a2a3a4a5a6a7a8; load('Hsca2ca2.mat')

Hsca2ca2=(a4+j*a5);f=a0;clear a0a1a2a3a4a5a6a7a8; load('Hsca2ca1.mat')

Hsca2ca1=(a4+j*a5);f=a0;clear a0a1a2a3a4a5a6a7a8; load('Hsca1cb2.mat')

Hsca1cb2=(a4+j*a5);f=a0;clear a0a1a2a3a4a5a6a7a8; load('Hsca1cb1.mat')

Hsca1cb1=(a4+j*a5);f=a0;clear a0a1a2a3a4a5a6a7a8; load('Hsca1ca2.mat')

Hsca1ca2=(a4+j*a5);f=a0;clear a0a1a2a3a4a5a6a7a8; load('Hsca1ca1.mat')

Hsca1ca1=(a4+j*a5);f=a0;clear a0a1a2a3a4a5a6a7a8; load('Hoa2ca2.mat')

Hoa2ca2=(a4+j*a5);f=a0;clear a0a1a2a3a4a5a6a7a8; load('Hoa2ca1.mat')

Hoa2ca1=(a4+j*a5);f=a0;clear a0a1a2a3a4a5a6a7a8; load('Hoa1ca2.mat')

Hoa1ca2=(a4+j*a5);f=a0;clear a0a1a2a3a4a5a6a7a8; load('Hoa1ca1.mat')

Hoa1ca1=(a4+j*a5);f=a0;clear a0a1a2a3a4a5a6a7a8;

load('Hcb2ib2.mat')

Hcb2ib2=(a4+j*a5);f=a0;clear a0a1a2a3a4a5a6a7a8;

load('Hcb2ib1.mat')

Hcb2ib1=(a4+j*a5);f=a0;clear a0a1a2a3a4a5a6a7a8;

load('Hcb2cb2.mat')

Hcb2cb2=(a4+j*a5);f=a0;clear a0a1a2a3a4a5a6a7a8;

load('Hcb2cb1.mat')

Hcb2cb1=(a4+j*a5);f=a0;clear a0a1a2a3a4a5a6a7a8;

load('Hcb1ib2.mat')

Hcb1ib2=(a4+j*a5);f=a0;clear a0a1a2a3a4a5a6a7a8;

load('Hcb1ib1.mat')

Hcb1ib1=(a4+1i*a5);f=a0;clear a0a1a2a3a4a5a6a7a8;

load('Hcb1cb2.mat')

Hcb1cb2=(a4+j*a5);f=a0;clear a0a1a2a3a4a5a6a7a8;

load('Hcb1cb1.mat')

Hcb1cb1=(a4+j*a5);f=a0;clear a0a1a2a3a4a5a6a7a8;

load('Hca2ca2.mat')

Hca2ca2=(a4+j*a5);f=a0;clear a0a1a2a3a4a5a6a7a8;

load('Hca2ca1.mat')

Hca2ca1=(a4+j*a5);f=a0;clear a0a1a2a3a4a5a6a7a8;

load('Hca1ca2.mat')

Hca1ca2=(a4+j*a5);f=a0;clear a0a1a2a3a4a5a6a7a8;

load('Hca1ca1.mat')

Hca1ca1=(a4+j*a5);f=a0;clear a0a1a2a3a4a5a6a7a8;

for k=2:401;w=2*pi*f

[Hcaca]=[Hca1ca1(k),Hca1ca2(k);Hca2ca1(k),Hca2ca2(k)]/(-w(k)^2); [Hcbcb]=[Hcb1cb1(k),Hcb1cb2(k);Hcb2cb1(k),Hcb2cb2(k)]/(-w(k)^2); [Hcbib]=[Hcb1ib1(k),Hcb1ib2(k);Hcb2ib1(k),Hcb2ib2(k)]/(-w(k)^2); [Hoaca]=[Hoa1ca1(k),Hoa1ca2(k);Hoa2ca1(k),Hoa2ca2(k)]/(-w(k)^2);

[Hscaca]=[Hsca1ca1(k),Hsca1ca2(k);Hsca2ca1(k),Hsca2ca2(k)]/(-w(k)^2);

[Hscacb]=[Hsca1cb1(k),Hsca1cb2(k);Hsca2cb1(k),Hsca2cb2(k)]/(-w(k)^2);

[Hscbcb]=[Hscb1cb1(k),Hscb1cb2(k);Hscb2cb1(k),Hscb2cb2(k)]/(-w(k)^2);

[Hsoaib]=[Hsoa1ib1(k),Hsoa1ib2(k);Hsoa2ib1(k),Hsoa2ib2(k)]/(-w(k)^2); KsD=inv(Hscaca*transpose(inv(Hscacb))*Hscbcb-Hscacb);

Ks11D(k)=KsD(1,1);Ks22D(k)=KsD(2,2);