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

登錄注冊

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

24小時熱門版塊排行榜

返回列表

三星2010

木蟲 (小有名氣)

應(yīng)助: 2 (幼兒園)
金幣: 1394.5
散金: 70
紅花: 1
帖子: 197
在線: 95.4小時
蟲號: 1210684
注冊: 2011-02-23
性別: GG
專業(yè): 力學(xué)中的基本問題和方法

[求助] 請matlab高手看一下這個問題怎樣編程才好，附程序代碼，謝謝！

主要工作就是求解隱函數(shù)方程，參數(shù)及方程如下：
taof=7.2;
Gf=0.576;
delta1=0.034;
deltaf=0.16;
b1=25;
E1=250000;
t1=5.12;
b2=400;
E2=10000;
t2=20;
L=1200;
eql=600;
beta=b2*E2*t2/(b1*E1*t1);
lamda=((taof^2/2/Gf*(1/E1/t1+b1/b2/E2/t2)))^0.5;
lamda1=(lamda^2*2*Gf/(delta1*taof))^0.5;
lamda2=(lamda^2*2*Gf/(deltaf-delta1)/taof)^0.5;
參數(shù)a和h滿足一下關(guān)系l
amda2/lamda1*sinh(lamda1*(L-a-h))- cot(lamda2*h)*cosh(lamda1*(L-a-h)-1/beta/sin(lamda2*h)=0
我能肯定的是a從0單調(diào)增大至892.806169，h的范圍大致是[135,150]我想要的結(jié)果是將a的變化區(qū)間若干等分，比如1000等分，對每一個a去求對應(yīng)的h并將這兩個參數(shù)輸出.
我自己編寫的程序如下：
for i=1:1:892.9/0.001
a(i)=i*0.001-0.001;
error=1;
  for j=135/0.001:1:150/0.001
h1(j)=j*0.001;
error1(j)=lamda2/lamda1*sinh(lamda1*(L-a(i)-h1(j)))- ...
            cot(lamda2*h1(j))*cosh(lamda1*(L-a(i)-h1(j)))- ...
            1/beta/sin(lamda2*h1(j));
   if error>abs(error1(j))
error=abs(error1(j));
   h(i)=h1(j);
   end
end
if L break;
end
ltpu(i)=lamda2/lamda1*tanh(lamda1*(L-h(i)-a(i)))*cos(lamda2*h(i))+sin(lamda2*h(i));
ltpb(i)=lamda2/((lamda1)^2*delta1*b1*E1*t1)*(1+cos(lamda2*h(i))/(beta*cosh(lamda1*(L-a(i)-h(i)))));
ltp(i)=ltpu(i)/ltpb(i);
lts(i)=deltaf+ltp(i)*a(i)/(b1*E1*t1);
end
plot(lts,ltp,'b');
程序運行后總提示h未定義，不知道是哪里出了問題，有沒有更好的辦法求解啊，我都被這個問題折磨暈了

，請大家出手相助！

[ Last edited by 三星2010 on 2012-12-26 at 16:36 ]

回復(fù)此樓

» 猜你喜歡

能源材料化學(xué)課題組招收碩士研究生8-10名已經(jīng)有6人回復(fù)
296求調(diào)劑已經(jīng)有5人回復(fù)
0703化學(xué)調(diào)劑，六級已過，有科研經(jīng)歷已經(jīng)有8人回復(fù)
268求調(diào)劑已經(jīng)有8人回復(fù)
材料，紡織，生物（0856、0710），化學(xué)招生啦已經(jīng)有7人回復(fù)
301求調(diào)劑已經(jīng)有5人回復(fù)
281求調(diào)劑（0805）已經(jīng)有7人回復(fù)
材料與化工求調(diào)劑已經(jīng)有6人回復(fù)
299求調(diào)劑已經(jīng)有4人回復(fù)
277調(diào)劑已經(jīng)有6人回復(fù)

» 本主題相關(guān)價值貼推薦，對您同樣有幫助:

matlab編程中遇到的問題已經(jīng)有4人回復(fù)
我想請問一個matlab編程的問題，sin（v/2）=3*v/14/pi 已經(jīng)有10人回復(fù)
請問如何用Matlab編程做這個圖？謝謝大家！已經(jīng)有5人回復(fù)
求高人指點 MATLAB 編程求解洛倫茲的最大李雅普諾夫指數(shù)程序！�。�！已經(jīng)有5人回復(fù)
100BB急求matlab編程中的問題解決，如何修改這個程序呢？已經(jīng)有4人回復(fù)
【求助】vb與matlab混合編程時圖形顯示窗口中，怎樣使圖形比窗口小一些？已經(jīng)有3人回復(fù)
【求助】用matlab編程，需要將數(shù)據(jù)自檢分類（高手進來指點下）已經(jīng)有5人回復(fù)
【求助】請問matlab編程語言如何打包成為獨立的應(yīng)用程序�。恐x謝已經(jīng)有8人回復(fù)

靜心，合力

1樓 2012-12-26 16:34:01

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

chembetsey

木蟲 (小有名氣)

應(yīng)助: 125 (高中生)
金幣: 3820.1
散金: 5
紅花: 9
帖子: 262
在線: 281.8小時
蟲號: 1781694
注冊: 2012-04-27
專業(yè): 理論和計算化學(xué)

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

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
三星2010: 金幣+10, ★有幫助, 我不能確定h到底是多少維的，可以這樣預(yù)定義嗎？ 2012-12-27 11:30:00

h是一個一維矩陣，你需要先定義一下，如：h=zeros[1000]

贊一下

回復(fù)此樓

2樓2012-12-26 23:57:16

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

chembetsey

木蟲 (小有名氣)

應(yīng)助: 125 (高中生)
金幣: 3820.1
散金: 5
紅花: 9
帖子: 262
在線: 281.8小時
蟲號: 1781694
注冊: 2012-04-27
專業(yè): 理論和計算化學(xué)

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

你可以定義一個足夠大的就是了，雖然有點浪費內(nèi)存。
如果你知道維數(shù)大小，就定義成指定的大小。

贊一下

回復(fù)此樓

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

三星2010

3樓2012-12-28 00:49:29

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

dbb627

榮譽版主 (著名寫手)

專家經(jīng)驗: +4
程序強帖: 2
應(yīng)助: 289 (大學(xué)生)
貴賓: 0.589
金幣: 24640.4
散金: 551
紅花: 61
沙發(fā): 1
帖子: 1246
在線: 1794.8小時
蟲號: 149791
注冊: 2005-12-29
性別: GG
專業(yè): 污染控制化學(xué)
管轄: 計算模擬

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

感謝參與，應(yīng)助指數(shù) +1

>> taof=7.2;
Gf=0.576;
delta1=0.034;
deltaf=0.16;
b1=25;
E1=250000;
t1=5.12;
b2=400;
E2=10000;
t2=20;
L=1200;
eql=600;
beta=b2*E2*t2/(b1*E1*t1);
lamda=((taof^2/2/Gf*(1/E1/t1+b1/b2/E2/t2)))^0.5;
lamda1=(lamda^2*2*Gf/(delta1*taof))^0.5;
lamda2=(lamda^2*2*Gf/(deltaf-delta1)/taof)^0.5;
syms a h
f=lamda2/lamda1*sinh(lamda1*(L-a-h))- cot(lamda2*h)*cosh(lamda1*(L-a-h))-1/beta/sin(lamda2*h);
ezplot(f,[0 900])

贊一下

回復(fù)此樓

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

三星2010

The more you learn, the more you know, the more you know, and the more you forget. The more you forget, the less you know. So why bother to learn.

4樓2012-12-28 09:05:17

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

三星2010

木蟲 (小有名氣)

應(yīng)助: 2 (幼兒園)
金幣: 1394.5
散金: 70
紅花: 1
帖子: 197
在線: 95.4小時
蟲號: 1210684
注冊: 2011-02-23
性別: GG
專業(yè): 力學(xué)中的基本問題和方法

送鮮花一朵

引用回帖:

4樓: Originally posted by dbb627 at 2012-12-28 09:05:17
>> taof=7.2;
Gf=0.576;
delta1=0.034;
deltaf=0.16;
b1=25;
E1=250000;
t1=5.12;
b2=400;
E2=10000;
t2=20;
L=1200;
eql=600;
beta=b2*E2*t2/(b1*E1*t1);
lamda=((taof^2/2/Gf*(1/E1/t1+b1/b2/ ...

dbb627
你好，按你的建議畫圖試了一下，并不是期望的結(jié)果，因為a和h都有對應(yīng)的物理量含義并且是一一對應(yīng)的，而圖中出現(xiàn)的是多值對應(yīng)，很遺憾，這里仍表示感謝！

贊一下

回復(fù)此樓

靜心，合力

5樓2012-12-28 13:29:32

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

三星2010

木蟲 (小有名氣)

應(yīng)助: 2 (幼兒園)
金幣: 1394.5
散金: 70
紅花: 1
帖子: 197
在線: 95.4小時
蟲號: 1210684
注冊: 2011-02-23
性別: GG
專業(yè): 力學(xué)中的基本問題和方法

送鮮花一朵

引用回帖:

3樓: Originally posted by chembetsey at 2012-12-28 00:49:29
你可以定義一個足夠大的就是了，雖然有點浪費內(nèi)存。
如果你知道維數(shù)大小，就定義成指定的大小。

謝謝你再次回復(fù)！我的疑問是假設(shè)求解前預(yù)定義h=zeros[100]
但實際計算后只用到了前50維，那后面h是自動調(diào)整為50維還是后面的50維元素是都置0（h維數(shù)不變）？，前者的話就好了，如果是后者的話就不是期望的情形。

贊一下

回復(fù)此樓

靜心，合力

6樓2012-12-28 13:39:03

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

dbb627

榮譽版主 (著名寫手)

專家經(jīng)驗: +4
程序強帖: 2
應(yīng)助: 289 (大學(xué)生)
貴賓: 0.589
金幣: 24640.4
散金: 551
紅花: 61
沙發(fā): 1
帖子: 1246
在線: 1794.8小時
蟲號: 149791
注冊: 2005-12-29
性別: GG
專業(yè): 污染控制化學(xué)
管轄: 計算模擬

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

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
三星2010: 金幣+30, ★★★很有幫助, 百度中也有個賬號是dbb627，莫非就是閣下，真乃高人也，你的代碼我看了一陣才明白，不過還是明白了，非常感謝！ 2012-12-28 19:48:07

引用回帖:

5樓: Originally posted by 三星2010 at 2012-12-28 13:29:32
dbb627
你好，按你的建議畫圖試了一下，并不是期望的結(jié)果，因為a和h都有對應(yīng)的物理量含義并且是一一對應(yīng)的，而圖中出現(xiàn)的是多值對應(yīng)，很遺憾，這里仍表示感謝！...

taof=7.2;
Gf=0.576;
delta1=0.034;
deltaf=0.16;
b1=25;
E1=250000;
t1=5.12;
b2=400;
E2=10000;
t2=20;
L=1200;
eql=600;
beta=b2*E2*t2/(b1*E1*t1);
lamda=((taof^2/2/Gf*(1/E1/t1+b1/b2/E2/t2)))^0.5;
lamda1=(lamda^2*2*Gf/(delta1*taof))^0.5;
lamda2=(lamda^2*2*Gf/(deltaf-delta1)/taof)^0.5;
syms a h
f=@(a,h)lamda2/lamda1*sinh(lamda1*(L-a-h))- cot(lamda2*h)*cosh(lamda1*(L-a-h))-1/beta/sin(lamda2*h);
x=0:5:890;
h=zeros(numel(x),1);
for i=1:numel(x)
opt=optimset('Display','off');
a=x(i);
h(i)=fsolve(@(h)f(a,h),135,opt);
end
plot(x,h)

贊一下(1人)

回復(fù)此樓

The more you learn, the more you know, the more you know, and the more you forget. The more you forget, the less you know. So why bother to learn.

7樓2012-12-28 14:39:01

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

三星2010

木蟲 (小有名氣)

應(yīng)助: 2 (幼兒園)
金幣: 1394.5
散金: 70
紅花: 1
帖子: 197
在線: 95.4小時
蟲號: 1210684
注冊: 2011-02-23
性別: GG
專業(yè): 力學(xué)中的基本問題和方法

引用回帖:

7樓: Originally posted by dbb627 at 2012-12-28 14:39:01
taof=7.2;
Gf=0.576;
delta1=0.034;
deltaf=0.16;
b1=25;
E1=250000;
t1=5.12;
b2=400;
E2=10000;
t2=20;
L=1200;
eql=600;
beta=b2*E2*t2/(b1*E1*t1);
lamda=((taof^2/2/Gf*(1/E1/t1+b1/b2/E2/t2))) ...

這個運行的很好！結(jié)果也可以接受，只是我覺得少點什么，比如若將來碰到的問題含有三個參數(shù)，該方法能否比較容易推廣至三維情形？最主要的是我現(xiàn)在還沒有弄明白我編的程序問題出在何處。。。

贊一下

回復(fù)此樓

靜心，合力

8樓2012-12-28 20:10:29

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

615780236

禁蟲 (知名作家)

★
感謝參與，應(yīng)助指數(shù) +1
jjdg: 金幣+1, 辛苦了 2012-12-30 01:55:12
jjdg: 編輯內(nèi)容 2012-12-30 01:55

本帖內(nèi)容被屏蔽

9樓2012-12-28 20:47:05

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

三星2010

木蟲 (小有名氣)

應(yīng)助: 2 (幼兒園)
金幣: 1394.5
散金: 70
紅花: 1
帖子: 197
在線: 95.4小時
蟲號: 1210684
注冊: 2011-02-23
性別: GG
專業(yè): 力學(xué)中的基本問題和方法

引用回帖:

9樓: Originally posted by 615780236 at 2012-12-28 20:47:05
taof=7.2;
Gf=0.576;
delta1=0.034;
deltaf=0.16;
b1=25;
E1=250000;
t1=5.12;
b2=400;
E2=10000;
t2=20;
L=1200;
eql=600;
beta=b2*E2*t2/(b1*E1*t1);
lamda=((taof^2/2/Gf*(1/E1/t1+b1/b2/E2/t2)))^ ...

你和dbb627一個人嗎，怎么寫一樣的代碼

贊一下

回復(fù)此樓

靜心，合力

10樓2012-12-29 09:43:57

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

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

返回列表

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

[考研] 293求調(diào)劑 +7	zjl的號 2026-03-16	12/600	2026-03-17 18:22 by 重科小霸王
[考研] 0703化學(xué)336分求調(diào)劑 +4	zbzihdhd 2026-03-15	5/250	2026-03-17 17:33 by ruiyingmiao
[考研] 085601求調(diào)劑 +4	Du.11 2026-03-16	4/200	2026-03-17 17:08 by ruiyingmiao
[考研] 材料與化工專碩調(diào)劑 +5	heming3743 2026-03-16	5/250	2026-03-17 14:03 by 勇敢太監(jiān)王公公
[考研] 機械專碩325，尋找調(diào)劑院校 +3	y9999 2026-03-15	5/250	2026-03-16 19:58 by y9999
[考研] 333求調(diào)劑 +3	文思客 2026-03-16	7/350	2026-03-16 18:21 by 文思客
[考研] 283求調(diào)劑 +10	小樓。 2026-03-12	14/700	2026-03-16 16:08 by 13811244083
[考研] 0856專碩279求調(diào)劑 +5	加油加油！? 2026-03-15	5/250	2026-03-15 11:58 by 2020015
[考研] 085601材料工程315分求調(diào)劑 +3	yang_0104 2026-03-15	3/150	2026-03-15 10:58 by peike
[考研] 中科大材料專碩319求調(diào)劑 +3	孟鑫材料 2026-03-13	3/150	2026-03-14 18:10 by houyaoxu
[考研] 材料371求調(diào)劑 +9	鱷魚? 2026-03-11	11/550	2026-03-13 22:53 by JourneyLucky
[考研] 材料與化工求調(diào)劑一志愿 985 總分 295 +8	dream…… 2026-03-12	8/400	2026-03-13 22:17 by 星空星月
[考研] 求調(diào)劑（材料與化工327） +4	愛吃香菜啦 2026-03-11	4/200	2026-03-13 22:11 by JourneyLucky
[考研] 四川大學(xué)085601材料工程專碩初試294求調(diào)劑 +4	祝我們好在冬天 2026-03-11	4/200	2026-03-13 21:39 by peike
[考研] 0703化學(xué)求調(diào)劑 +7	綠豆芹菜湯 2026-03-12	7/350	2026-03-13 17:25 by njzyff
[考研] 求b區(qū)學(xué)校調(diào)劑 +3	周56 2026-03-11	3/150	2026-03-13 16:20 by JourneyLucky
[考研] 考研調(diào)劑 +4	芬達46 2026-03-12	4/200	2026-03-13 16:04 by ruiyingmiao
[考研] 一志愿211化學(xué)學(xué)碩310分求調(diào)劑 +8	努力奮斗112 2026-03-12	9/450	2026-03-13 15:41 by JourneyLucky
[考研] 274求調(diào)劑 +3	S.H1 2026-03-12	3/150	2026-03-13 15:15 by JourneyLucky
[考研] 0856化學(xué)工程280分求調(diào)劑 +4	shenzxsn 2026-03-11	4/200	2026-03-13 11:55 by ymwdoctor

24小時熱門版塊排行榜

[求助] 請matlab高手看一下這個問題怎樣編程才好，附程序代碼，謝謝！

» 猜你喜歡

» 本主題相關(guān)價值貼推薦，對您同樣有幫助:

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

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

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

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

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

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

[求助] 請matlab高手看一下這個問題怎樣編程才好，附程序代碼，謝謝！

» 本主題相關(guān)價值貼推薦，對您同樣有幫助: