版塊導(dǎo)航: 正在加載中...

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號(hào)能夠正常使用，請(qǐng)盡快對(duì)帳號(hào)進(jìn)行手機(jī)號(hào)驗(yàn)證，感謝您的理解與支持！

24小時(shí)熱門版塊排行榜

返回列表

當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖

longwen8

木蟲 (文壇精英)

應(yīng)助: 0 (幼兒園)
金幣: 5385.8
散金: 13079
紅花: 48
沙發(fā): 83
帖子: 10804
在線: 1496.5小時(shí)
蟲號(hào): 4254310
注冊(cè): 2015-11-29
性別: GG
專業(yè): 環(huán)境工程

[求助] 微分方程組的四階龍格庫塔公式求解matlab版已有1人參與

微分方程組的四階龍格庫塔公式，求解matlab版的代碼。下面是我的代碼，運(yùn)行不出來。
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時(shí)的誤差曲線’);

回復(fù)此樓

» 猜你喜歡

推薦一款可以AI輔助寫作的Latex編輯器SmartLatexEditor，超級(jí)好用，推薦試試已經(jīng)有11人回復(fù)
2026-CJ開始申報(bào)了已經(jīng)有1人回復(fù)
物理學(xué)I論文潤色/翻譯怎么收費(fèi)? 已經(jīng)有229人回復(fù)
西安電子科學(xué)大學(xué)杭州研究院劉麗香教授招收智能多模態(tài)傳感器和微型儲(chǔ)能器件方向博士已經(jīng)有13人回復(fù)
重慶交通大學(xué)光子學(xué)微結(jié)構(gòu)與器件課題組招收2026年碩士研究生信息已經(jīng)有1人回復(fù)
一志愿鄭大材料學(xué)碩298分，求調(diào)劑已經(jīng)有6人回復(fù)
尋合作：應(yīng)力腐蝕多尺度模擬已經(jīng)有3人回復(fù)
考研交流已經(jīng)有0人回復(fù)
【新加坡】納米電子器件項(xiàng)目組有“聯(lián)合培養(yǎng)博士生”名額已經(jīng)有0人回復(fù)

1樓 2020-12-13 22:11:20

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

longwen8

木蟲 (文壇精英)

應(yīng)助: 0 (幼兒園)
金幣: 5385.8
散金: 13079
紅花: 48
沙發(fā): 83
帖子: 10804
在線: 1496.5小時(shí)
蟲號(hào): 4254310
注冊(cè): 2015-11-29
性別: GG
專業(yè): 環(huán)境工程

引用回帖:

2樓: Originally posted by hzlhm at 2020-12-18 18:12:48
一、幾個(gè)自定義函數(shù)位置錯(cuò)誤。應(yīng)為這樣放置
1、主程序，
function rungekutta()
clc,clear
x0=0;xn=1.0;y0=1;h=0.05;%h為步長
=rungekutta4(x0,xn,y0,h);
end
2、四階龍格庫塔函數(shù)程序
function =rungeku ...

什么結(jié)果？大神。

回復(fù)此樓

3樓2020-12-18 19:43:08

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

查看全部 8 個(gè)回答

hzlhm

至尊木蟲 (著名寫手)

應(yīng)助: 387 (碩士)
金幣: 17807
紅花: 53
帖子: 2879
在線: 606.3小時(shí)
蟲號(hào): 1517335
注冊(cè): 2011-11-30
性別: GG
專業(yè): 常微分方程與動(dòng)力系統(tǒng)

【答案】應(yīng)助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
獨(dú)孤神宇: 金幣+5, 鼓勵(lì)交流 2020-12-18 19:25:10
longwen8: 金幣+20, ★★★★★最佳答案 2020-12-18 19:42:37

一、幾個(gè)自定義函數(shù)位置錯(cuò)誤。應(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ù)書寫錯(cuò)誤。應(yīng)為
xlabel('x軸'),ylabel('誤差');
三、修改后運(yùn)行得到如下結(jié)果

贊一下(1人)

回復(fù)此樓

QQ：2120156492

2樓2020-12-18 18:12:48

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

hzlhm

至尊木蟲 (著名寫手)

應(yīng)助: 387 (碩士)
金幣: 17807
紅花: 53
帖子: 2879
在線: 606.3小時(shí)
蟲號(hào): 1517335
注冊(cè): 2011-11-30
性別: GG
專業(yè): 常微分方程與動(dòng)力系統(tǒng)

引用回帖:

3樓: Originally posted by longwen8 at 2020-12-18 19:43:08
什么結(jié)果？大神。...

運(yùn)行結(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
圖像上傳不上來。

贊一下

回復(fù)此樓

QQ：2120156492

4樓2020-12-18 22:02:15

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

longwen8

木蟲 (文壇精英)

應(yīng)助: 0 (幼兒園)
金幣: 5385.8
散金: 13079
紅花: 48
沙發(fā): 83
帖子: 10804
在線: 1496.5小時(shí)
蟲號(hào): 4254310
注冊(cè): 2015-11-29
性別: GG
專業(yè): 環(huán)境工程

引用回帖:

4樓: Originally posted by hzlhm at 2020-12-18 22:02:15
運(yùn)行結(jié)果：
i    x(i)    yi數(shù)值解       y(i)真值    誤差
-------------------------------------------------------------
1    0.0500    1.04880886    1.04880885 0.00000001
2    ...

這么多結(jié)果，選哪一個(gè)。還是這種運(yùn)算，就沒有盡頭啊

發(fā)自小木蟲Android客戶端

贊一下

回復(fù)此樓

5樓2020-12-19 00:00:51

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

查看全部 8 個(gè)回答

最具人氣熱帖推薦 [查看全部]		作者	回/看	最后發(fā)表

[考研] 化學(xué)工程求調(diào)劑 +5	化工人999 2026-03-04	5/250	2026-03-04 16:02 by sslc1985
[考研] 290求調(diào)劑 +3	Bananaiy 2026-03-04	3/150	2026-03-04 15:52 by sslc1985
[考研] 材料化工調(diào)劑 +15	今夏不夏 2026-03-01	18/900	2026-03-04 15:44 by 每天只擺一小會(huì)
[考研] 281求調(diào)劑 +3	我是小小蔥蔥 2026-03-03	5/250	2026-03-04 14:23 by kakakapanpan
[考研] 347求調(diào)劑 +5	啊歐歐歐 2026-03-03	5/250	2026-03-04 14:10 by kakakapanpan
[考研] 一志愿西工大材料學(xué)，英一數(shù)二總分321分，求調(diào)劑。 +3	zz05zz 2026-03-02	5/250	2026-03-04 09:16 by lbsjt
[考研] 理學(xué)，工學(xué)，農(nóng)學(xué)調(diào)劑，少走彎路，這里歡迎您！ +8	likeihood 2026-03-02	11/550	2026-03-03 19:39 by 一葉凡塵?
[考研] 276求調(diào)劑 +8	路lyh123 2026-02-28	10/500	2026-03-03 18:25 by xin吖
[考研] 江蘇省農(nóng)科院招調(diào)劑1名 +5	Qwertyuop 2026-03-01	5/250	2026-03-03 17:50 by 梅子鴨鴨
[考研] 化工270求調(diào)劑 +10	什么名字qwq 2026-03-02	10/500	2026-03-03 13:54 by NUAAZXWS
[考研] 268求調(diào)劑 +6	好運(yùn)連綿不絕 2026-03-02	6/300	2026-03-03 13:03 by 秋收
[考研] 324求調(diào)劑 +4	wxz2 2026-03-03	5/250	2026-03-03 09:25 by 楊楊楊紫
[考研] 307求調(diào)劑 +6	wyyyqx 2026-03-01	6/300	2026-03-03 09:24 by 2235787770
[考研] 調(diào)劑材料學(xué)碩 +4	詞凝Y 2026-03-02	4/200	2026-03-03 07:53 by njzyff
[考研] 材料復(fù)試調(diào)劑 +5	學(xué)材料的點(diǎn) 2026-03-01	6/300	2026-03-02 20:01 by hypershenger
[考研] 材料085601調(diào)劑 +5	多多子. 2026-03-02	5/250	2026-03-02 19:15 by zhukairuo
[考研] 295求調(diào)劑 +8	19171856320 2026-02-28	8/400	2026-03-02 11:19 by yuchj
[考研] 調(diào)劑 +3	13853210211 2026-03-02	4/200	2026-03-02 10:16 by 13853210211
[考研] 317一志愿華南理工電氣工程求調(diào)劑 +6	Soliloquy_Q 2026-02-28	11/550	2026-03-01 11:14 by 歌liekkas
[考研] 307求調(diào)劑 +4	73372112 2026-02-28	6/300	2026-03-01 00:04 by ll247

亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产 高清 中文字幕,99re热久久亚洲综合精品成人,熟妇 一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频

24小時(shí)熱門版塊排行榜

longwen8

[求助] 微分方程組的四階龍格庫塔公式求解matlab版 已有1人參與

» 猜你喜歡

longwen8

hzlhm

【答案】應(yīng)助回帖

hzlhm

longwen8

亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产高清中文字幕,99re热久久亚洲综合精品成人,熟妇一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频

[求助] 微分方程組的四階龍格庫塔公式求解matlab版已有1人參與