程序大致如下,主要是想通过求解超越方程,将解得的结果用于循环。但是在超越方程部分编程出现问题,想请大家帮忙完成注释部分的超越方程的编写。谢谢!【公式稍微有点长】
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)