变压器负载运行仿真时遇到的问题

wxlhhm

已知， ，下面是我编的程序
SN=10e3;U1N=380;U2N=220;r1=0.14;r2=0.035;x1=0.22;x2=0.055;rm=30;xm=310;ZL=4+j*3;
I1N=SN/U1N;
I2N=SN/U2N;
k=U1N/U2N;
Z1=r1+j*x1;
rr2=k^2*r2;
xx2=k^2*x2;
ZZ2=rr2+j*xx2;
ZZL=k^2*ZL;
Zm=rm+j*xm;
Zd=Z1+1/(1/Zm+1/(ZZ2+ZZL));
U11=U1N;
I11=U11/Zd;
E11=(U11-I11*Z1);
I221=E11/(ZZ2+ZZL);
I21=k*I221;
U221=I221*ZZL;
U21=U221/k;
cospsi1=cos(angle(Zd));
cospsi2=cos(angle(ZL));
P1=abs(U11)*abs(I11)*cospsi1;
P2=abs(U21)*abs(I21)*cospsi2;
eta=P2/P1;
Im1=E11/Zm;
pFe=abs(Im1)^2*rm;
pCu1=abs(I11)^2*r1;
pCu2=abs(I21)^2*r2;
disp('原边电流='),disp(abs(I11))
disp('副边电流='),disp(abs(I21))
disp('副边电压='),disp(abs(U21))
disp('原边功率因数='),disp(cospsi1)
disp('原边功率='),disp(P1)
disp('副边功率因数='),disp(cospsi2)
disp('副边功率='),disp(P2)
disp('效率='),disp(eta)
disp('励磁电流='),disp(abs(Im1))
disp('铁损耗='),disp(pFe)
disp=('原边铜损耗='),disp(pCu1)
disp=('副边铜损耗='),disp(pCu2)
第一次将{disp=('原边铜损耗='),disp(pCu1)
disp=('副边铜损耗='),disp(pCu2)}误打成了
disp=('原边铜损耗='),disp(pcu1)
disp=('副边铜损耗='),disp(pcu2),即Cu打成了cu,运行结果为：
原边电流=
   25.5752
副边电流=
   42.7447
副边电压=
  213.7237
原边功率因数=
    0.7725
原边功率=
  7.5072e+003
副边功率因数=
    0.8000
副边功率=
  7.3085e+003
效率=
    0.9735
励磁电流=
    1.1998
铁损耗=
   43.1836

disp =
原边铜损耗=
??? Undefined function or variable "pcu1".
Error in ==> bianyaqifuzaiyunxing at 38
disp=('原边铜损耗='),disp(pcu1)
可是改正后，在调试却调试不出来了，一直出来下面的情况：
Error in ==> bianyaqifuzaiyunxing at 28
disp('原边电流='),disp(abs(I11))
??? Index exceeds matrix dimensions.
哪位高手帮忙解决下啊，谢谢了！

yjhou

你的E11计算的结果是E11 =

  3.7366e+002 -2.0724e+000i，提示说是下标必须是实数或者逻辑数。

wxlhhm

你好，谢谢版主，我把I11=U11/Zd改成I11=abs(U11/Zd)，调试下，还是出现了disp =
原边铜损耗=
??? Undefined function or variable "pcu1".
Error in ==> bianyaqifuzaiyunxing at 38
disp=('原边铜损耗='),disp(pcu1)
可是改正后，在调试却调试不出来了，一直出来下面的情况：
Error in ==> bianyaqifuzaiyunxing at 28
disp('原边电流='),disp(abs(I11))
??? Index exceeds matrix dimensions.
这样的提示啊？能不能帮我改下程序啊！谢谢啊！我是一个初学者！麻烦了！