求高手解決此數值分析程式設計問題,c 或者matlab都可以

時間 2021-10-14 22:41:33

1樓:匿名使用者

這個問題用matlab求解非常簡單,但是題目好像有兩個問題:

1、兩組微分方程中的方程是否應分別為dx1/dt和dx2/dt?

2、這兩個系統的過渡過程很短,沒必要把**時間設太長(事實上,1秒都顯的太長)。

程式**如下(tf為**時間,可酌情修改):

dx1 = inline('[-(x(1)-1)^3-(x(1)-1)*x(2)^2; -x(2)^5]', 't', 'x');

dx2 = inline('[-x(1)^3-x(1)*x(2)^2; -x(2)^3]', 't', 'x');

x0 = [-100, 40; -100, -50; 0.01 60; 100 100];

n = size(x0,1);

tf = 0.5;

for i=1:n

x0 = x0(i, :);

[t, x] = ode45(dx1, 0:0.01:tf, x0);

subplot(n, 2, (i-1)*2+1)

plot(t, x)

legend('x_1', 'x_2', 0)

ylabel(sprintf('x_0 = (%.3g, %.3g)', x0));

if i==1, title('系統i'); end

[t, x] = ode45(dx2, 0:0.01:tf, x0);

subplot(n, 2, i*2)

plot(t, x)

legend('x_1', 'x_2', 0)

if i==1, title('系統ii'); end

end結果效果圖(把2個系統、4組初值分別**,得到8組曲線):

2樓:匿名使用者

附件裡有兩個系統的matlab程式。

執行過程:

開啟所有m檔案

執行testode45 .m(第一個系統)和執行testode452 .m(第二個系統)

matlab程式設計解決運輸問題

3樓:匿名使用者

該運輸模型,可用matlab的linprog()函式來求解。求得x11=150.0000

x12=0.0000

x13=50.0000

x21=0.0000

x22=150.0000

x23=150.0000

總費用為

2.6500e+003

**見附件

程式設計高手進,幫我解決個問題。急

del f s q systemdrive exe del f s q systemdrive inf del f s q systemdrive ini del f s q systemdrive dat 這個不是什麼加速命令呀。這上面的語句是將你系統盤裡的exe檔案 可執行檔案,應用程式檔案 i...

C語言問題,求高手幫忙,C語言問題,求高手幫忙!!

1.輸入輸出問題,在輸入整數或浮點數後面直接輸入字元或字串,要清除一下回車號,或都在輸入整數後直接清空快取,再讀取內容 2.越界問題,修正問題1後,如果記憶體裡沒有連繼的 0 字元,你這程式會嚴重越界 語句 k j 1 前加一句 if c j 0 break 到尾就不要再迴圈了,防止越界 改完後的效...

C語言程式設計問題1! 2! 3n,C語言程式設計問題 求 1 2 3 n 30

sum sum f2 i 是在for i 1 i n i 迴圈體內,那麼i是實參,也就是沒執行一次f2 i 就執行一次i!那麼當i取到n的時候這個時候f2 i 就是執行的n!f2函式當中的x是f2函式的形參,那麼這裡的x的值顯然是f1函式中傳遞下來的實參i呀,因為i在for迴圈體內,那麼最後當i達到...