matlab解复杂超越方程组,matlab解超越方程组求助！ - 计算模拟 - 小木虫 - 学术 科研 互动社区...

再贴一下文件中代码，方面观看。

lc

clear

close all

y0 = [1e9;4e20;0];

tspan=[0 1e-7];

tic

[t,y] = ode23('rate_eq',tspan,y0);

toc

y(:,1) = max(y(:,1),1);

figure

subplot(3,1,1);

plot(t,y(:,1));

xlabel(' 时间(s)');

ylabel('光子数密度(m^{-3})');

subplot(3,1,2);

plot(t,y(:,2));

xlabel('时间(s)');

ylabel('反转粒子束密度(m^{-3})');

subplot(3,1,3);

plot(t,y(:,3));

xlabel('时间(s)');

ylabel('基态粒子数密度{-3})');

function dy = rate_eq(t,y)

sigma = 7.6e-23;

sigma_o = 1e-20;

sigma_s = 2.3e-21;

n1 = 1.972;

n2 = 3.48;

L = 0.1;

l = 0.006;

d = 0.00058;

R = 0.85;

c = 2.997963e8;

lc = n1*l+n2*d+(0.23-l-d);

tr = 2*lc/c;

n0 = 1.2e22;

beta = 2.6e-10;h = 6.623e-34;lamda = 1.06e-6;niu = c/lamda;

B = 6*beta*h*niu*c*(0.3)^2;%此项可以为0；

y(1) = max(y(1),1);

dy = [  y(1)*(2*sigma*y(2)*l-2*sigma_s*y(3)*d-2*sigma_o*(n0-y(3))*d-B*y(1)-(log(1/R)+L))/tr;

-4*sigma*y(1)*y(2)/(h*niu);

4*y(1)*(sigma_o*(n0-y(3))-sigma_s*y(3))/(h*niu)];

请虫友帮忙！