close all;home;clc;clear all;
%Inputs
As0=-1i*2.0658;
As1=2.4874;
As2=-1i*2.1950;
As3=2.4015;
As4=0;
Bs0=-0.1059;
Bs1=-1i*4.1687;
Bs2=4.4575;
Bs3=-1i*1.5183;
Bs4=2;
Cs0=-0.1476;
Cs1=-1i*3.0823;
Cs2=2.8836;
Cs3=0;
Cs4=0;
Ds0=-1i*2.0658;
Ds1=2.4874;
Ds2=-1i*2.1950;
Ds3=2.4015;
Ds4=0;
Pse0=2.0696;
Pse1=1i*2.7104;
Pse2=-0.8660;
Pse3=0;
Pse4=0;
%Polys
S1_As0=-1i*Cs0;
S1_As1=-1i*Cs1;
S1_As2=-1i*Cs2;
S1_As3=-1i*Cs3;
S1_As4=-1i*Cs4;
S1_Bs0=-1i*Ds0;
S1_Bs1=-1i*Ds1;
S1_Bs2=-1i*Ds2;
S1_Bs3=-1i*Ds3;
S1_Bs4=-1i*Ds4;
S1_Cs0=-1i*As0;
S1_Cs1=-1i*As1;
S1_Cs2=-1i*As2;
S1_Cs3=-1i*As3;
S1_Cs4=-1i*As4;
S1_Ds0=-1i*Bs0;
S1_Ds1=-1i*Bs1;
S1_Ds2=-1i*Bs2;
S1_Ds3=-1i*Bs3;
S1_Ds4=-1i*Bs4;
S1_As(:,1)=[S1_As0,S1_As1,S1_As2,S1_As3,S1_As4];
S1_Bs(:,1)=[S1_Bs0,S1_Bs1,S1_Bs2,S1_Bs3,S1_Bs4];
S1_Cs(:,1)=[S1_Cs0,S1_Cs1,S1_Cs2,S1_Cs3,S1_Cs4];
S1_Ds(:,1)=[S1_Ds0,S1_Ds1,S1_Ds2,S1_Ds3,S1_Ds4];
S1_Ps(:,1)=[Pse0,Pse1,Pse2,Pse3,Pse4];
C1 =S1_Ds(5,1)/S1_Bs(4,1);
S2_Cs0=S1_Cs0;
S2_Cs1=S1_Cs1-C1*S1_As0;
S2_Cs2=S1_Cs2-C1*S1_As1;
S2_Cs3=S1_Cs3-C1*S1_As2;
S2_Cs4=S1_Cs4-C1*S1_As3;
S2_Ds0=S1_Ds0;
S2_Ds1=S1_Ds1-C1*S1_Bs0;
S2_Ds2=S1_Ds2-C1*S1_Bs1;
S2_Ds3=S1_Ds3-C1*S1_Bs2;
S2_Ds4=S1_Ds4-C1*S1_Bs3;
S2_As(:,1)=[S1_As0,S1_As1,S1_As2,S1_As3,S1_As4];
S2_Bs(:,1)=[S1_Bs0,S1_Bs1,S1_Bs2,S1_Bs3,S1_Bs4];
S2_Cs(:,1)=[S2_Cs0,S2_Cs1,S2_Cs2,S2_Cs3,S2_Cs4];
S2_Ds(:,1)=[S2_Ds0,S2_Ds1,S2_Ds2,S2_Ds3,S2_Ds4];
S2_Ps(:,1)=[Pse0,Pse1,Pse2,Pse3,Pse4];
B1 =(S2_Ds(4,1)/S2_Bs(4,1))/1i;
S3_Cs0=S2_Cs0-1i*B1*S1_As0;
S3_Cs1=S2_Cs1-1i*B1*S1_As1;
S3_Cs2=S2_Cs2-1i*B1*S1_As2;
S3_Cs3=S2_Cs3-1i*B1*S1_As3;
S3_Cs4=S2_Cs4-1i*B1*S1_As4;
S3_Ds0=S2_Ds0-1i*B1*S1_Bs0;
S3_Ds1=S2_Ds1-1i*B1*S1_Bs1;
S3_Ds2=S2_Ds2-1i*B1*S1_Bs2;
S3_Ds3=S2_Ds3-1i*B1*S1_Bs3;
S3_Ds4=S2_Ds4-1i*B1*S1_Bs4;
S3_As(:,1)=[S1_As0,S1_As1,S1_As2,S1_As3,S1_As4];
S3_Bs(:,1)=[S1_Bs0,S1_Bs1,S1_Bs2,S1_Bs3,S1_Bs4];
S3_Cs(:,1)=[S3_Cs0,S3_Cs1,S3_Cs2,S3_Cs3,S3_Cs4]
S3_Ds(:,1)=[S3_Ds0,S3_Ds1,S3_Ds2,S3_Ds3,S3_Ds4]
S3_Ps(:,1)=[Pse0,Pse1,Pse2,Pse3,Pse4];
コメント