| 查看: 1619 | 回复: 0 | ||
[求助]
超越方程求解的实现
|
|
程序大致如下，主要是想通过求解超越方程，将解得的结果用于循环。但是在超越方程部分编程出现问题，想请大家帮忙完成注释部分的超越方程的编写。谢谢！【公式稍微有点长】 clear;clc L=1e-3; C=40e-6; E=48; T=0.0002; R=5; T1=0.0001;T2=0.0001;e=0; % R=13;%注意该值的正负性I=6; for m=1:300 K=0.1+0.001*(m-1); i(1)=0; v(1)=0; a1=(((- 4*C*R^2 + L)/L)^(1/2) - 1)/(2*C*R); a2=-(((- 4*C*R^2 + L)/L)^(1/2) + 1)/(2*C*R); for n=1:1000 %%%%%%%%%%%%%%%%%%%%%%%% %在该处求解超越方程，解得的结果T1，从而可以得到T2=T-T1然后再进行下面的循环计算. %超越方程为1.5-24*K-K*X1-3*T1;其中的X1 如下： %X1=(E*(L*((a2*exp(T1*a1))/(a1 - a2) - (a1*exp(T1*a2))/(a1 - a2) +(exp(T1*a1)/(a1 - a2) - exp(T1*a2)/(a1 - a2))/(C*R) + 1) - (L*(exp(T1*a1)/(a1 - a2) - exp(T1*a2)/(a1 - a2)))/(C*R)))/L - v(n)*((a2*exp(T1*a1))/(a1 - a2) - (a1*exp(T1*a2))/(a1 - a2) + (exp(T1*a1)/(a1 - a2) - exp(T1*a2)/(a1 - a2))/(C*R)) + (i(n)*(exp(T1*a1)/(a1 - a2) - exp(T1*a2)/(a1 - a2)))/C; %%%%%%%%%%%%%%%%%%%%%%%% i(n+1)= i(n)*(((a1*exp(T1*a2))/(a1 - a2) - (a2*exp(T1*a1))/(a1 - a2))*((a1*exp(T2*a2))/(a1 - a2) - (a2*exp(T2*a1))/(a1 - a2)) - ((exp(T1*a1)/(a1 - a2) - exp(T1*a2)/(a1 - a2))*(exp(T2*a1)/(a1 - a2) - exp(T2*a2)/(a1 - a2)))/(C*L)) - v(n)*((((a1*exp(T2*a2))/(a1 - a2) - (a2*exp(T2*a1))/(a1 - a2))*(exp(T1*a1)/(a1 - a2) - exp(T1*a2)/(a1 - a2)))/L - ((exp(T2*a1)/(a1 - a2) - exp(T2*a2)/(a1 - a2))*((a2*exp(T1*a1))/(a1 - a2) - (a1*exp(T1*a2))/(a1 - a2) + (exp(T1*a1)/(a1 - a2) - exp(T1*a2)/(a1 - a2))/(C*R)))/L) + (E*(L*((((a1*exp(T2*a2))/(a1 - a2) - (a2*exp(T2*a1))/(a1 - a2))*(exp(T1*a1)/(a1 - a2) - exp(T1*a2)/(a1 - a2)))/L - ((exp(T2*a1)/(a1 - a2) - exp(T2*a2)/(a1 - a2))*((a2*exp(T1*a1))/(a1 - a2) - (a1*exp(T1*a2))/(a1 - a2) + (exp(T1*a1)/(a1 - a2) - exp(T1*a2)/(a1 - a2))/(C*R) + 1))/L) + (L*(((a1*exp(T2*a2))/(a1 - a2) - (a2*exp(T2*a1))/(a1 - a2))*((a2*exp(T1*a1))/(a1 - a2) - (a1*exp(T1*a2))/(a1 - a2) + 1) + ((exp(T1*a1)/(a1 - a2) - exp(T1*a2)/(a1 - a2))*(exp(T2*a1)/(a1 - a2) - exp(T2*a2)/(a1 - a2)))/(C*L)))/R))/L; v(n+1)=v(n)*(((a2*exp(T1*a1))/(a1 - a2) - (a1*exp(T1*a2))/(a1 - a2) + (exp(T1*a1)/(a1 - a2) - exp(T1*a2)/(a1 - a2))/(C*R))*((a2*exp(T2*a1))/(a1 - a2) - (a1*exp(T2*a2))/(a1 - a2) + (exp(T2*a1)/(a1 - a2) - exp(T2*a2)/(a1 - a2))/(C*R)) - ((exp(T1*a1)/(a1 - a2) - exp(T1*a2)/(a1 - a2))*(exp(T2*a1)/(a1 - a2) - exp(T2*a2)/(a1 - a2)))/(C*L)) + i(n)*((((a1*exp(T1*a2))/(a1 - a2) - (a2*exp(T1*a1))/(a1 - a2))*(exp(T2*a1)/(a1 - a2) - exp(T2*a2)/(a1 - a2)))/C - ((exp(T1*a1)/(a1 - a2) - exp(T1*a2)/(a1 - a2))*((a2*exp(T2*a1))/(a1 - a2) - (a1*exp(T2*a2))/(a1 - a2) + (exp(T2*a1)/(a1 - a2) - exp(T2*a2)/(a1 - a2))/(C*R)))/C) - (E*(L*(((a2*exp(T2*a1))/(a1 - a2) - (a1*exp(T2*a2))/(a1 - a2) + (exp(T2*a1)/(a1 - a2) - exp(T2*a2)/(a1 - a2))/(C*R))*((a2*exp(T1*a1))/(a1 - a2) - (a1*exp(T1*a2))/(a1 - a2) + (exp(T1*a1)/(a1 - a2) - exp(T1*a2)/(a1 - a2))/(C*R) + 1) - ((exp(T1*a1)/(a1 - a2) - exp(T1*a2)/(a1 - a2))*(exp(T2*a1)/(a1 - a2) - exp(T2*a2)/(a1 - a2)))/(C*L)) - (L*(((exp(T2*a1)/(a1 - a2) - exp(T2*a2)/(a1 - a2))*((a2*exp(T1*a1))/(a1 - a2) - (a1*exp(T1*a2))/(a1 - a2) + 1))/C + ((exp(T1*a1)/(a1 - a2) - exp(T1*a2)/(a1 - a2))*((a2*exp(T2*a1))/(a1 - a2) - (a1*exp(T2*a2))/(a1 - a2) + (exp(T2*a1)/(a1 - a2) - exp(T2*a2)/(a1 - a2))/(C*R)))/C))/R))/L; end end |
回复此楼» 猜你喜欢
材料与化工（0856）304求 B区 调剂
已经有3人回复
-
268求调剂
已经有8人回复
-
265求调剂
已经有14人回复
-
350 本科985求调剂，求老登收留
已经有3人回复
-
【考研调剂】化学专业 281分，一志愿四川大学，诚心求调剂
已经有9人回复
-
初试 317
已经有3人回复
-
279分求调剂 一志愿211
已经有15人回复
-
286求调剂
已经有4人回复
-
279求调剂
已经有4人回复
-
求调剂
已经有3人回复
