版塊導(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í)熱門版塊排行榜

返回列表

li52

金蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 1399
散金: 90
紅花: 1
帖子: 54
在線: 19.7小時(shí)
蟲號(hào): 3714002
注冊(cè): 2015-03-06
性別: GG
專業(yè): 固體力學(xué)

[求助] Matlab求解二階常微分方程邊值問(wèn)題已有2人參與

我用bvp4c求解了一個(gè)簡(jiǎn)單的二階常微分方程邊值問(wèn)題
原方程 $\frac{\mathrmk77yzwn^{2}y}{\mathrmafbdewxx^{2}}+k\cdot \cos y=0$ ，其中k是一個(gè)常數(shù)
邊界條件 $y\left ( 0 \right )=0, \frac{\mathrmoz2ufgiy\left ( 1 \right )}{\mathrm82qb2oqx}=0$
我想得到不同的k時(shí)的結(jié)果。需要從別的結(jié)果中提取k的值，數(shù)量較大，一個(gè)一個(gè)修改不太現(xiàn)實(shí)。
據(jù)我嘗試，描述微分方程的function的參數(shù)只能是(x,y)，不能再添加一個(gè)k。
怎樣能整理成一個(gè)function，把k也當(dāng)作參數(shù)呢？

附：k=1的時(shí)候，MATLAB計(jì)算過(guò)程

先轉(zhuǎn)換成一階微分方程組
$\frac{\mathrmycbcxrxy_{1}}{\mathrmdm27zlix}=y_{2}$
$\frac{\mathrmexj27rpy_{2}}{\mathrmig7lctkx}=-\cos y_{1}$
$y_{1}\left ( 0 \right )=0$ ， $y_{2}\left ( 1 \right )=0$

求解過(guò)程：
1.magFun.m  描述微分方程
  function dy = magFun(x,y)
  k = 1;
  dy = zeros(2,1);
  dy(1) = y(2);
  dy(2) = -k*cos(y(1));
  end
2.magBC.m  描述邊界條件
  function res = magBC( ya,yb )
  res = [ ya(1); yb(2) ];
  end
3.magSol.m  求解
  function magSol
  solinit = bvpinit(linspace(0,1,5),[1 0]);
  sol = bvp4c(@magFun,@magBC,solinit);
  x = linspace(0,1);
  y = deval(sol,x);
  plot(x,y(1,: ))

[ Last edited by li52 on 2016-10-8 at 17:10 ]

回復(fù)此樓

» 收錄本帖的淘帖專輯推薦

程序

» 猜你喜歡

0702一志愿吉大B區(qū)求調(diào)劑已經(jīng)有5人回復(fù)
求調(diào)劑已經(jīng)有0人回復(fù)
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)? 已經(jīng)有165人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑有論文已經(jīng)有0人回復(fù)
請(qǐng)問(wèn)還有沒(méi)有用Latex寫文章的小伙伴們？已經(jīng)有0人回復(fù)
光學(xué)工程學(xué)碩調(diào)劑信息已經(jīng)有32人回復(fù)
歡迎加入課題組已經(jīng)有0人回復(fù)
散金幣，求好運(yùn)，祝面上順利！已經(jīng)有76人回復(fù)
華東師范大學(xué)芯片設(shè)計(jì)徐瓏真實(shí)水平咋樣已經(jīng)有1人回復(fù)
調(diào)劑已經(jīng)有0人回復(fù)
【新加坡】納米電子器件項(xiàng)目組有“聯(lián)合培養(yǎng)博士生”名額已經(jīng)有0人回復(fù)

1樓 2016-10-08 15:48:15

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

dingd

鐵桿木蟲 (職業(yè)作家)

計(jì)算強(qiáng)帖: 4
應(yīng)助: 1641 (講師)
金幣: 15037.3
散金: 101
紅花: 234
帖子: 3410
在線: 1223.7小時(shí)
蟲號(hào): 291104
注冊(cè): 2006-10-28

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

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
li52: 金幣+10 2016-10-09 09:37:29

可以試試用1stOpt計(jì)算。假如k范圍=【1,4】，變幅0.1

CODE:

LoopConstant k=[1:0.1:4];

Variable x=[0,1],y=0,y'=[,0];

Plot x[x],y,y';

ODEFunction y''=-k*cos(y);

c261.jpg

c262.jpg

贊一下(2人)

回復(fù)此樓

2樓2016-10-08 19:50:27

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

li52

金蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 1399
散金: 90
紅花: 1
帖子: 54
在線: 19.7小時(shí)
蟲號(hào): 3714002
注冊(cè): 2015-03-06
性別: GG
專業(yè): 固體力學(xué)

引用回帖:

2樓: Originally posted by dingd at 2016-10-08 19:50:27
可以試試用1stOpt計(jì)算。假如k范圍=【1,4】，變幅0.1

LoopConstant k=;
Variable x=,y=0,y'=;
Plot x,y,y';
ODEFunction y''=-k*cos(y);

c261.jpg

c262.jpg
...

謝謝！
Google了下，1stOpt應(yīng)該也是一種數(shù)學(xué)軟件吧。但是之前的計(jì)算都是在MATLAB上進(jìn)行的，再重新學(xué)習(xí)新軟件，然后在軟件間轉(zhuǎn)換數(shù)據(jù)確實(shí)比較頭疼。

贊一下

回復(fù)此樓

3樓2016-10-08 21:21:59

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

yz457694

鐵桿木蟲 (正式寫手)

計(jì)算強(qiáng)帖: 1
應(yīng)助: 28 (小學(xué)生)
金幣: 7193
散金: 1000
紅花: 16
帖子: 721
在線: 168.7小時(shí)
蟲號(hào): 2820058
注冊(cè): 2013-11-22
性別: GG
專業(yè): 傳遞過(guò)程

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

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
li52: 金幣+40, ★★★★★最佳答案 2016-10-09 09:37:09

function dy = magFun(x,y)
  global k;
  dy = zeros(2,1);
  dy(1) = y(2);
  dy(2) = -k*cos(y(1));
end


function res = magBC( ya,yb )
  res = [ ya(1); yb(2) ];
end

  function magSol
  clear,clc
  global k; %全局變量k
  yt=[]; %存儲(chǔ)要繪圖的y數(shù)據(jù)
  x = linspace(0,1,10);
  for k=1:5
  solinit = bvpinit(x,[1 0]);
  sol = bvp4c(@magFun,@magBC,solinit);

  y = deval(sol,x);
  yt=[yt;y(1,

];
  end
  plot(x,yt)
  legend('1','2','3','4','5')

figure.jpg

贊一下

回復(fù)此樓

4樓2016-10-09 07:16:05

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

yz457694

鐵桿木蟲 (正式寫手)

計(jì)算強(qiáng)帖: 1
應(yīng)助: 28 (小學(xué)生)
金幣: 7193
散金: 1000
紅花: 16
帖子: 721
在線: 168.7小時(shí)
蟲號(hào): 2820058
注冊(cè): 2013-11-22
性別: GG
專業(yè): 傳遞過(guò)程

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

★ ★
dbb627: 金幣+2, 謝謝指導(dǎo) 2016-10-09 08:59:54

表情干擾了。。。大致上就是把k設(shè)為了全局變量，然后加個(gè)循環(huán)，就可以自己設(shè)置k的范圍來(lái)得到不同結(jié)果了

CODE:

function dy = magFun(x,y)

  global k;

  dy = zeros(2,1);

  dy(1) = y(2);

  dy(2) = -k*cos(y(1));

end

function res = magBC( ya,yb )

  res = [ ya(1); yb(2) ];

end

  function magSol

  clear,clc

  global k; %全局變量k

  yt=[];    %存儲(chǔ)要繪圖的y數(shù)據(jù)

  x = linspace(0,1,10);

  for k=1:5

  solinit = bvpinit(x,[1 0]);

  sol = bvp4c(@magFun,@magBC,solinit);

  y = deval(sol,x);

  yt=[yt;y(1,:)];

  end

  plot(x,yt)

  legend('1','2','3','4','5')

贊一下

回復(fù)此樓

» 本帖已獲得的紅花（最新10朵）

li52

5樓2016-10-09 07:19:13

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

li52

金蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 1399
散金: 90
紅花: 1
帖子: 54
在線: 19.7小時(shí)
蟲號(hào): 3714002
注冊(cè): 2015-03-06
性別: GG
專業(yè): 固體力學(xué)

送紅花一朵

引用回帖:

5樓: Originally posted by yz457694 at 2016-10-09 07:19:13
表情干擾了。。。大致上就是把k設(shè)為了全局變量，然后加個(gè)循環(huán)，就可以自己設(shè)置k的范圍來(lái)得到不同結(jié)果了

function dy = magFun(x,y)
  global k;
  dy = zeros(2,1);
  dy(1) = y(2);
  dy(2) = -k*cos(y(1)) ...

謝謝！
寫代碼的時(shí)候就糾結(jié)變量在各個(gè)function里，不能全局使用，原來(lái)可以定義成全局變量！基礎(chǔ)沒(méi)有打好啊
PS : 加上 ) 就成了表情了，幾個(gè)空格就好了。

: )

贊一下

回復(fù)此樓

6樓2016-10-09 09:33:56

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

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 li52 的主題更新

返回列表

普通表情龍兔虎貓

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

[考研] 招08考數(shù)學(xué) +5	laoshidan 2026-03-20	12/600	2026-03-23 12:38 by 梨花珞晚風(fēng)
[考研] 291 求調(diào)劑 +4	化工2026屆畢業(yè)?/a> 2026-03-21	4/200	2026-03-23 09:59 by bingchuan
[考研] 一志愿華中農(nóng)業(yè)071010，總分320求調(diào)劑 +5	困困困困坤坤 2026-03-20	6/300	2026-03-22 17:41 by hxsm
[考研] 一志愿西北大學(xué) ，070300化學(xué)學(xué)碩，總分287，雙非一本，求調(diào)劑。 +3	晨昏線與星海 2026-03-20	3/150	2026-03-22 16:00 by ColorlessPI
[考研] 求調(diào)劑 +7	Auroracx 2026-03-22	7/350	2026-03-22 12:38 by 素顏傾城1988
[考研] 考研調(diào)劑 +4	來(lái)好運(yùn)來(lái)來(lái)來(lái) 2026-03-21	4/200	2026-03-22 12:15 by 星空星月
[考研] 求助 +5	夢(mèng)里的無(wú)言 2026-03-21	6/300	2026-03-21 17:51 by 學(xué)員8dgXkO
[考研] 299求調(diào)劑 +5	shxchem 2026-03-20	7/350	2026-03-21 17:09 by ColorlessPI
[考研] 求調(diào)劑 +3	白QF 2026-03-21	3/150	2026-03-21 13:12 by zhukairuo
[考研] 332求調(diào)劑 +3	鳳凰院丁真 2026-03-20	3/150	2026-03-21 10:27 by luoyongfeng
[考研] 南昌大學(xué)材料專碩311分求調(diào)劑 +6	77chaselx 2026-03-20	6/300	2026-03-21 07:24 by JourneyLucky
[考研] 一志愿天津大學(xué)化學(xué)工藝專業(yè)（081702）315分求調(diào)劑 +12	yangfz 2026-03-17	12/600	2026-03-21 03:30 by JourneyLucky
[考研] 274求調(diào)劑 +10	S.H1 2026-03-18	10/500	2026-03-20 23:51 by JourneyLucky
[考研] 一志愿南京理工大學(xué)085701資源與環(huán)境302分求調(diào)劑 +4	葵梓衛(wèi)隊(duì) 2026-03-18	6/300	2026-03-20 23:02 by JourneyLucky
[考研] 招收調(diào)劑碩士 +4	lidianxing 2026-03-19	12/600	2026-03-20 12:25 by lidianxing
[考研] 085410人工智能專碩317求調(diào)劑（0854都可以） +4	xbxudjdn 2026-03-18	4/200	2026-03-20 09:07 by 不168
[考研] 288求調(diào)劑，一志愿華南理工大學(xué)071005 +5	ioodiiij 2026-03-17	5/250	2026-03-19 18:22 by zcl123
[考研] 收復(fù)試調(diào)劑生 +4	雨后秋荷 2026-03-18	4/200	2026-03-18 14:16 by elevennnne
[考研] 材料，紡織，生物（0856、0710），化學(xué)招生啦 +3	Eember. 2026-03-17	9/450	2026-03-18 10:28 by Eember.
[考博] 26博士申請(qǐng) +3	1042136743 2026-03-17	3/150	2026-03-17 23:30 by 輕松不少隨