查看: 286|回复: 0|关注: 0

[未答复] matlab求解非线性方程组

[复制链接]

新手

10 麦片

财富积分


050


1

主题

2

帖子

0

最佳答案
发表于 2019-3-13 19:11:13 | 显示全部楼层 |阅读模式
求大神帮忙,Cc  Cp1 Cp12为变量,分别为4.03e-6 4.04e-6 88e-12 ;H(1)、H(2)、H(3)为方程组,
无法通过解方程得到三个变量的值;试过fsolve,比solve好一点,但结果也不对,求大神帮忙~
k=57.735/12949;

C2=0.04885e-6/(k^2);
Rs=766*k^2;Ls=164.96*k^2;
Rm=7.141800e6;Lm=4.004459e4;
Rt1=4192.1*k^2;Lt1=5.3373*k^2;
Rt2=0.0033333;Lt2=0.0001061;
w=2*pi*100;
syms Cc Cp1 Cp12
Z1=1/(j*w*C2);Z2=Rs+j*w*Ls;
Z3=(Rt1+j*Lt1*w)+(Rt2+j*Lt2*w)+(Rt1+j*Lt1*w)*(Rt2+j*Lt2*w)/(Rm+j*Lm*w);
Z4=(Rt1+j*Lt1*w)+(Rm+j*Lm*w)+(Rt1+j*Lt1*w)*(Rm+j*Lm*w)/(Rt2+j*Lt2*w);
Z5=(Rm+j*Lm*w)+(Rt2+j*Lt2*w)+(Rm+j*Lm*w)*(Rt2+j*Lt2*w)/(Rt1+j*Lt1*w);
X1=1/(j*w*Cc);X2=1/(j*w*Cp1);X3=1/(j*w*Cp12);
A=[Z1+(X1*Z2)/(X1+Z2)];B=(X3*Z3)/(X3+Z3);C=(X2+Z4)/(X2*Z4);
H(1)=-(k/Z1)*[A+B+A*B*C]-(0.0107-j*6.9912e-4)
k=57.735/12949;
C2=0.04885e-6/(k^2);
Rs=766*k^2;Ls=164.96*k^2;
Rm=7.141800e6;Lm=4.004459e4;
Rt1=4192.1*k^2;Lt1=5.3373*k^2;
Rt2=0.0033333;Lt2=0.0001061;
w=2*pi*200;
syms Cc Cp1 Cp12
Z1=1/(j*w*C2);Z2=Rs+j*w*Ls;
Z3=(Rt1+j*Lt1*w)+(Rt2+j*Lt2*w)+(Rt1+j*Lt1*w)*(Rt2+j*Lt2*w)/(Rm+j*Lm*w);
Z4=(Rt1+j*Lt1*w)+(Rm+j*Lm*w)+(Rt1+j*Lt1*w)*(Rm+j*Lm*w)/(Rt2+j*Lt2*w);
Z5=(Rm+j*Lm*w)+(Rt2+j*Lt2*w)+(Rm+j*Lm*w)*(Rt2+j*Lt2*w)/(Rt1+j*Lt1*w);
X1=1/(j*w*Cc);X2=1/(j*w*Cp1);X3=1/(j*w*Cp12);
A=[Z1+(X1*Z2)/(X1+Z2)];B=(X3*Z3)/(X3+Z3);C=(X2+Z4)/(X2*Z4);
H(2)=-(k/Z1)*[A+B+A*B*C]-(0.0571-j*0.0014)

k=57.735/12949;
C2=0.04885e-6/(k^2);
Rs=766*k^2;Ls=164.96*k^2;
Rm=7.141800e6;Lm=4.004459e4;
Rt1=4192.1*k^2;Lt1=5.3373*k^2;
Rt2=0.0033333;Lt2=0.0001061;
w=2*pi*300;
syms Cc Cp1 Cp12
Z1=1/(j*w*C2);Z2=Rs+j*w*Ls;
Z3=(Rt1+j*Lt1*w)+(Rt2+j*Lt2*w)+(Rt1+j*Lt1*w)*(Rt2+j*Lt2*w)/(Rm+j*Lm*w);
Z4=(Rt1+j*Lt1*w)+(Rm+j*Lm*w)+(Rt1+j*Lt1*w)*(Rm+j*Lm*w)/(Rt2+j*Lt2*w);
Z5=(Rm+j*Lm*w)+(Rt2+j*Lt2*w)+(Rm+j*Lm*w)*(Rt2+j*Lt2*w)/(Rt1+j*Lt1*w);
X1=1/(j*w*Cc);X2=1/(j*w*Cp1);X3=1/(j*w*Cp12);
A=[Z1+(X1*Z2)/(X1+Z2)];B=(X3*Z3)/(X3+Z3);C=(X2+Z4)/(X2*Z4);
H(3)=-(k/Z1)*[A+B+A*B*C]-(0.1393-j*0.0021)
C=solve(H(1),H(2),H(3),Cc,Cp1,Cp12);
Cc=vpa(C.Cc)
Cp1=vpa(C.Cp1)
Cp12=vpa(C.Cp12)
您需要登录后才可以回帖 登录 | 注册

本版积分规则

关闭

站长推荐上一条 /4 下一条

快速回复 返回顶部 返回列表