1樓:匿名使用者
#include
#include
#include
void main()
for(j=0;j<=8;j++)
u[8][0]=u[8][0]+0.5*0.5*0.5*0.5*0.5*0.5*0.5*0.5*0.5*0.5;
for(j=0;j<=8;j++)
}for(j=1;j<=8;j++)
cout<<"\n";}}
2樓:淘淘vs跑跑
這裡假設t=0時,x0=1,y0=2,z0=3.
function hh
[t,x]=ode23t(@fun,0:1,[1 2 3])plot(t,x)
function yhat=fun(t,x)dxdt=-0.85*x(1)+0.25*x(2)+0.35*x(3)-0.8*x(1)^2*x(2);
dydt=0.8*x(1)^2*x(2)-1.2*x(2);
dzdt=0.9*x(1)+0.8*x(2)-0.3*x(3);
yhat=[dxdt;dydt;dzdt];
結果:t =
0 0.0145
0.0435
0.0726
0.1120
0.1730
0.2341
0.2951
0.3718
0.4485
0.5435
0.6386
0.7337
0.8287
0.9238
1.0000
x =1.0000 2.0000 3.0000
0.9879 1.9879 3.0228
0.9653 1.9629 3.0677
0.9459 1.9363 3.1110
0.9239 1.8985 3.1675
0.8979 1.8374 3.2500
0.8799 1.7748 3.3269
0.8680 1.7121 3.3984
0.8595 1.6348 3.4816
0.8565 1.5599 3.5578
0.8583 1.4717 3.6437
0.8645 1.3893 3.7212
0.8738 1.3127 3.7912
0.8852 1.2421 3.8546
0.8980 1.1772 3.9121
0.9090 1.1291 3.9544
3樓:
用dsolve求不出解析解,用數值解法吧,最好能提供初始條件,這裡假設t=0時,x0=1,y0=2,z0=3.
function hh
[t,x]=ode23t(@fun,0:1,[1 2 3])plot(t,x)
function yhat=fun(t,x)dxdt=-0.85*x(1)+0.25*x(2)+0.35*x(3)-0.8*x(1)^2*x(2);
dydt=0.8*x(1)^2*x(2)-1.2*x(2);
dzdt=0.9*x(1)+0.8*x(2)-0.3*x(3);
yhat=[dxdt;dydt;dzdt];
結果:t =
0 0.0145
0.0435
0.0726
0.1120
0.1730
0.2341
0.2951
0.3718
0.4485
0.5435
0.6386
0.7337
0.8287
0.9238
1.0000
x =1.0000 2.0000 3.0000
0.9879 1.9879 3.0228
0.9653 1.9629 3.0677
0.9459 1.9363 3.1110
0.9239 1.8985 3.1675
0.8979 1.8374 3.2500
0.8799 1.7748 3.3269
0.8680 1.7121 3.3984
0.8595 1.6348 3.4816
0.8565 1.5599 3.5578
0.8583 1.4717 3.6437
0.8645 1.3893 3.7212
0.8738 1.3127 3.7912
0.8852 1.2421 3.8546
0.8980 1.1772 3.9121
0.9090 1.1291 3.9544 一定對la
4樓:匿名使用者
la fang~~~~~
在matlan中執行時方法如下:
儲存檔案
function yhat=fun(t,x)dxdt=-0.85*x(1)+0.25*x(2)+0.35*x(3)-0.8*x(1)^2*x(2);
dydt=0.8*x(1)^2*x(2)-1.2*x(2);
dzdt=0.9*x(1)+0.8*x(2)-0.3*x(3);
yhat=[dxdt;dydt;dzdt];
然後輸入
t=0,x0=1,y0=2,z0=3.
[t,x]=ode23t(@fun,0:1,[1 2 3])plot(t,x)
答案就出來了
aa3301169
不要高估我們的智商~~~>_<
5樓:匿名使用者
游泳的天使22別吹拉
有本事答上來呀
6樓:游泳的天使
簡單總算我們數學專業的人一展身手了
7樓:咖啡是一種罪
真深奧~~~~~~~~~~
matlab求定積分解析解,matlab求一個定積分解析解
小學輔導資料小屋 這個求不出來的。warning explicit integral could not be found.e的表示式的積分找不到。只能通過數值計算的方法求出近似值,這個近似值也保持了足夠的精度 dx 0.01 x 0 dx 2 y exp exp x x s dx cumtrapz...
matlab不等式求解目標函式最優解
方程數大於未知量個數,屬於超定方程組,一般沒有精確解,只能求最小二乘解。matlab中可以利用左除或偽逆法求解。clear c1 10496.108 c2 10500.275 c3 10504.442 c4 10508.609 c5 10512.776 a 1,c1,c1 2,c1 3 1,c2,c...
如何用matlab求解0 1規劃問題
3 0 1型整數規劃 0 1型整數規劃是整數規劃中的特殊情形,它的變數j x 僅取值0 或1。這時j x 稱 為0 1變數,或稱二進位制變數。j x 僅取值0 或1 這個條件可由下述約束條件 0 1 j x 整數所代替,是和一般整數規劃的約束條件形式一致的。在實際問題中,如果引入0 1變 量,就可以...