要編的是這個方程,nT=88,C0,G0,Ck,Gk都是88*88矩陣,u(s)經(jīng)傅里葉變換為88*2049矩陣,B為單位陣,s為2049*1的矩陣。求X(s)。我編的程序如下,有錯誤,該怎么修改?求修改后程序。謝謝
f=Freq;%2049*1
s=j*w;%2049*1
w=2*pi*f;%2049*1
c=3e8;
epslonr=8.85e-12*36;
mur=4*pi*1e-7;
n=88;%未知變量個數(shù)
nT=88;%指數(shù)項(xiàng)個數(shù)
C0=diag(diag(L));%自電感88*88
G0=diag(diag(Gt));%自電容88*88
%C1=randn(n);
%G1=randn(n);
% X(s)=[];%未知狀態(tài)變量矢量
B=eye(88,88);%匹配矩陣,單位陣?
%L=zeros(n,k);
% us=zeros(k,1);%網(wǎng)絡(luò)端口電壓
td1=dis./c*sqrt(epslonr*mur);%延時時間
C=[];%電感
G=[];%電容
g=Gt-G0;%電容
c=L-C0;%電感
sum=0;
us=zeros(88,4097)
us(1, =Lig_FFT
for k=1:88%指數(shù)項(xiàng)數(shù)量
C=[C;c(1:k,1:k)];
G=[G;g(1:k,1:k)];
C((k-1)*n+1:k*n, =c(1:k,1:k);
G((k-1)*n+1:k*n, =g(1:k,1:k);
sum=sum+(s*C+G)*exp(-s*td(1:k,1:k));
end
X(s)=inv(s*C0+G0+sum)*B*us;
% i(s)=L'*X(s);
錯誤提示如下
??? Error using ==> mtimes
Inner matrix dimensions must agree.
Error in ==> lashi at 349
sum=sum+(s*C+G)*exp(-s*td(1:k,1:k));
,內(nèi)部矩陣維數(shù)不一致該怎么修改?
![求大神幫我看看matlab我編的程序錯在哪里,怎么修改?]()
圖片1.png |