 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 2024 | 回復(fù): 7 | ||
| 當前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖 | ||
longwen8木蟲 (文壇精英)
|
[求助]
微分方程組的四階龍格庫塔公式求解matlab版 已有1人參與
|
|
|
微分方程組的四階龍格庫塔公式，求解matlab版的代碼。下面是我的代碼，運行不出來。 function varargout=rungekutta(varargin) clc,clear x0=0;xn=1.0;y0=1;h=0.05;%h為步長 [y,x]=rungekutta4(x0,xn,y0,h); Function z=f(x,y); z=y-2*x/y; function [y,x]=rungekutta4(x0,xn,y0,h) x=x0:h:xn; n=(xn-x0)/h; y1=x; y1(1)=y0; for i=1:n %龍格庫塔法的算法 K1=f(x(i),y1(i)); K2=f(x(i)+h/2,y1(i)+h/2*K1); K3= f(x(i)+h/2,y1(i)+h/2*K2); K4= f(x(i)+h,y1(i)+h*K3); y1(i+1)=y1(i)+h/6*(K1+2*K2+2*K3+K4); yy1(i)=(1+2*x(i+1)) ^0.5; error(i)=y1(i+1)-yy1(i); yy(i+1)=yy1(i); e(i+1)=error(i); end y=y1; n=(xn-x0)/h; fprintf(‘i x(i) yi數(shù)值解 y(i)真值 誤差\n’); fprintf(‘-------------------------------------------------------------\n’); for i=1:n fprintf(‘%2d %12.4f %14.8f %14.8f %12.8f\n’,i, x(i+1), y(i+1), yy(i+1), error(i)); end plot(x,e,‘r.’);%畫出誤差曲線 xlable(‘x軸’),ylable(‘誤差’); title(‘步長為0.05時的誤差曲線’); |
回復(fù)此樓» 猜你喜歡
西安電子科學(xué)大學(xué)杭州研究院劉麗香教授招收智能多模態(tài)傳感器和微型儲能器件方向博士
已經(jīng)有13人回復(fù)
-
重慶交通大學(xué)光子學(xué)微結(jié)構(gòu)與器件課題組招收2026年碩士研究生信息
已經(jīng)有1人回復(fù)
-
物理學(xué)I論文潤色/翻譯怎么收費?
已經(jīng)有131人回復(fù)
-
一志愿鄭大材料學(xué)碩298分，求調(diào)劑
已經(jīng)有6人回復(fù)
-
尋合作：應(yīng)力腐蝕多尺度模擬
已經(jīng)有3人回復(fù)
-
考研交流
已經(jīng)有0人回復(fù)
-
【新加坡】納米電子器件項目組有“聯(lián)合培養(yǎng)博士生”名額
已經(jīng)有0人回復(fù)
-
【新加坡】納米電子器件項目組有“聯(lián)合培養(yǎng)博士生”名額
已經(jīng)有0人回復(fù)
7樓2020-12-21 12:41:45
【答案】應(yīng)助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
獨孤神宇: 金幣+5, 鼓勵交流 2020-12-18 19:25:10
longwen8: 金幣+20, ★★★★★最佳答案 2020-12-18 19:42:37
獨孤神宇: 金幣+5, 鼓勵交流 2020-12-18 19:25:10
longwen8: 金幣+20, ★★★★★最佳答案 2020-12-18 19:42:37
|
一、幾個自定義函數(shù)位置錯誤。應(yīng)為這樣放置 1、主程序 ， function rungekutta() clc,clear x0=0;xn=1.0;y0=1;h=0.05;%h為步長 [y,x]=rungekutta4(x0,xn,y0,h); end 2、 四階龍格庫塔函數(shù)程序 function [y,x]=rungekutta4(x0,xn,y0,h) 。。。 end 3、z函數(shù)程序 function z=f(x,y); z=y-2*x/y; end 二、xlable(‘x軸’),ylable(‘誤差’);這句代碼中函數(shù)書寫錯誤。應(yīng)為 xlabel('x軸'),ylabel('誤差'); 三、修改后運行得到如下結(jié)果 |
2樓2020-12-18 18:12:48
longwen8
木蟲 (文壇精英)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 5385.8
- 散金: 13079
- 紅花: 48
- 沙發(fā): 83
- 帖子: 10804
- 在線: 1496.5小時
- 蟲號: 4254310
- 注冊: 2015-11-29
- 性別: GG
- 專業(yè): 環(huán)境工程
3樓2020-12-18 19:43:08
|
運行結(jié)果： i x(i) yi數(shù)值解 y(i)真值 誤差 ------------------------------------------------------------- 1 0.0500 1.04880886 1.04880885 0.00000001 2 0.1000 1.09544514 1.09544512 0.00000002 3 0.1500 1.14017546 1.14017543 0.00000004 4 0.2000 1.18321600 1.18321596 0.00000005 5 0.2500 1.22474493 1.22474487 0.00000006 6 0.3000 1.26491113 1.26491106 0.00000007 7 0.3500 1.30384056 1.30384048 0.00000008 8 0.4000 1.34164088 1.34164079 0.00000009 9 0.4500 1.37840498 1.37840488 0.00000011 10 0.5000 1.41421368 1.41421356 0.00000012 11 0.5500 1.44913781 1.44913767 0.00000014 12 0.6000 1.48323985 1.48323970 0.00000015 13 0.6500 1.51657526 1.51657509 0.00000017 14 0.7000 1.54919353 1.54919334 0.00000019 15 0.7500 1.58113904 1.58113883 0.00000021 16 0.8000 1.61245178 1.61245155 0.00000023 17 0.8500 1.64316793 1.64316767 0.00000026 18 0.9000 1.67332034 1.67332005 0.00000028 19 0.9500 1.70293895 1.70293864 0.00000031 20 1.0000 1.73205115 1.73205081 0.00000034 圖像上傳不上來。 |
4樓2020-12-18 22:02:15
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 085600材料與化工（高分子）290分求調(diào)劑 +4 | wengyujian 2026-03-04 | 4/200 |
|
|---|---|---|---|---|
|
|
[考研] 085701環(huán)境工程 求調(diào)劑 +6 | xiiiia 2026-03-04 | 6/300 |
|
|
[考研] 292求調(diào)劑 +9 | yhk_819 2026-02-28 | 9/450 |
|
|
|
[考研]
085600 英一數(shù)二272求調(diào)劑
 5+6
|
vida_a 2026-03-01 | 47/2350 |
|
|
|
[考研] 322,求調(diào)劑 +3 | 菜菜愛玩 2026-03-04 | 3/150 |
|
|
|
[考研] 求調(diào)劑院校 +6 | 云朵452 2026-03-02 | 11/550 |
|
|
|
[考研] 085600求調(diào)劑 +4 | LRZZZZZZ 2026-03-02 | 6/300 |
|
|
|
[論文投稿]
EST拒稿重投
5+3
|
15102603076 2026-03-02 | 3/150 |
|
|
|
[考研] 環(huán)境調(diào)劑 +8 | chenhanheng 2026-03-02 | 8/400 |
|
|
|
[考研] 283求調(diào)劑 +9 | 鹿沫笙 2026-03-02 | 9/450 |
|
|
|
[考研] 281電子信息求調(diào)劑 +3 | jhtfeybgj 2026-03-02 | 6/300 |
|
|
|
[考研] 沒上岸的看過來 +3 | tangxiaotian 2026-03-01 | 3/150 |
|
|
|
[考研] 0856求調(diào)劑285 +11 | 呂仔龍 2026-02-28 | 11/550 |
|
|
|
[考研] 306分材料調(diào)劑 +5 | chuanzhu川燭 2026-03-01 | 6/300 |
|
|
|
[考研] 材料085601調(diào)劑 +5 | 多多子. 2026-03-02 | 5/250 |
|
|
|
[考研] 材料調(diào)劑 +3 | 恒順自然 2026-03-02 | 3/150 |
|
|
|
[考研] 一志愿華南理工大學(xué)材料與化工326分，求調(diào)劑 +3 | wujinrui1 2026-02-28 | 3/150 |
|
|
|
[考研] 291分工科求調(diào)劑 +9 | science餓餓 2026-03-01 | 10/500 |
|
|
|
[考研] 0856材料求調(diào)劑 +4 | 麻辣魷魚 2026-02-28 | 4/200 |
|
|
|
[論文投稿]
Optics letters投稿被拒求助
30+3
|
luckyry 2026-02-26 | 4/200 |
|