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

返回列表

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

三星2010

木蟲(chóng) (小有名氣)

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

[求助] 請(qǐng)matlab高手看一下這個(gè)問(wèn)題怎樣編程才好，附程序代碼，謝謝！

主要工作就是求解隱函數(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滿(mǎn)足一下關(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等分，對(duì)每一個(gè)a去求對(duì)應(yīng)的h并將這兩個(gè)參數(shù)輸出.
我自己編寫(xiě)的程序如下：
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');
程序運(yùn)行后總提示h未定義，不知道是哪里出了問(wèn)題，有沒(méi)有更好的辦法求解啊，我都被這個(gè)問(wèn)題折磨暈了

，請(qǐng)大家出手相助！

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

回復(fù)此樓

» 猜你喜歡

336求調(diào)劑已經(jīng)有3人回復(fù)
306求調(diào)劑已經(jīng)有9人回復(fù)
工科0856求調(diào)劑已經(jīng)有4人回復(fù)
工科材料085601 279求調(diào)劑已經(jīng)有10人回復(fù)
081700 調(diào)劑 267分已經(jīng)有4人回復(fù)
276求調(diào)劑。有半年電池和半年高分子實(shí)習(xí)經(jīng)歷已經(jīng)有9人回復(fù)
263求調(diào)劑已經(jīng)有9人回復(fù)
求老師收我已經(jīng)有3人回復(fù)
招08考數(shù)學(xué) 已經(jīng)有12人回復(fù)
接收2026碩士調(diào)劑(學(xué)碩+專(zhuān)碩) 已經(jīng)有4人回復(fù)

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

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

靜心，合力

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

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

chembetsey

木蟲(chóng) (小有名氣)

應(yīng)助: 125 (高中生)
金幣: 3820.1
散金: 5
紅花: 9
帖子: 262
在線: 281.8小時(shí)
蟲(chóng)號(hào): 1781694
注冊(cè): 2012-04-27
專(zhuān)業(yè): 理論和計(jì)算化學(xué)

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

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

贊一下

回復(fù)此樓

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

三星2010

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

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

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

chembetsey

木蟲(chóng) (小有名氣)

應(yīng)助: 125 (高中生)
金幣: 3820.1
散金: 5
紅花: 9
帖子: 262
在線: 281.8小時(shí)
蟲(chóng)號(hào): 1781694
注冊(cè): 2012-04-27
專(zhuān)業(yè): 理論和計(jì)算化學(xué)

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

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

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

贊一下

回復(fù)此樓

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

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

dbb627

榮譽(yù)版主 (著名寫(xiě)手)

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

【答案】應(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

木蟲(chóng) (小有名氣)

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

送鮮花一朵

引用回帖:

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
你好，按你的建議畫(huà)圖試了一下，并不是期望的結(jié)果，因?yàn)閍和h都有對(duì)應(yīng)的物理量含義并且是一一對(duì)應(yīng)的，而圖中出現(xiàn)的是多值對(duì)應(yīng)，很遺憾，這里仍表示感謝！

贊一下

回復(fù)此樓

靜心，合力

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

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

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

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

[考研] 306求調(diào)劑 +9	chuanzhu川燭 2026-03-18	9/450	2026-03-23 13:17 by luoyongfeng
[考研] 263求調(diào)劑 +6	yqdszhdap－ 2026-03-22	9/450	2026-03-23 12:57 by yqdszhdap－
[考研] 招08考數(shù)學(xué) +5	laoshidan 2026-03-20	12/600	2026-03-23 12:38 by 梨花珞晚風(fēng)
[考研] 307求調(diào)劑 +3	余意卿 2026-03-21	3/150	2026-03-23 10:32 by Iveryant
[考研] 一志愿070300浙大化學(xué)358分，求調(diào)劑！ +4	酥酥魚(yú).. 2026-03-21	4/200	2026-03-23 08:12 by Iveryant
[考研] 289材料與化工（085600）B區(qū)求調(diào)劑 +3	這么名字咋樣 2026-03-22	4/200	2026-03-22 17:56 by 云民大李老師
[考研] 一志愿中南化學(xué)（0703）總分337求調(diào)劑 +9	niko- 2026-03-19	10/500	2026-03-22 16:08 by ColorlessPI
[考研] 260求調(diào)劑 +3	朱芷琳 2026-03-20	4/200	2026-03-22 15:12 by 朱芷琳
[考研] 求調(diào)劑 +5	Zhangbod 2026-03-21	7/350	2026-03-22 13:13 by Zhangbod
[考研] 一志愿華中科技大學(xué)071000，求調(diào)劑 +4	沿岸有貝殼6 2026-03-21	4/200	2026-03-22 07:21 by ilovexiaobin
[考研] 材料求調(diào)劑 +5	@taotao 2026-03-21	5/250	2026-03-21 20:55 by lbsjt
[考研] 考研調(diào)劑 +3	呼呼？~+123456 2026-03-21	3/150	2026-03-21 20:04 by 無(wú)際的草原
[考研] 一志愿南大，0703化學(xué)，分?jǐn)?shù)336，求調(diào)劑 +3	收到VS 2026-03-21	3/150	2026-03-21 18:42 by 學(xué)員8dgXkO
[考研] 南昌大學(xué)材料專(zhuān)碩311分求調(diào)劑 +6	77chaselx 2026-03-20	6/300	2026-03-21 07:24 by JourneyLucky
[考研] 346求調(diào)劑[0856] +4	WayneLim327 2026-03-16	7/350	2026-03-21 04:02 by JourneyLucky
[考研] 求調(diào)劑 +3	Ma_xt 2026-03-17	3/150	2026-03-21 02:05 by JourneyLucky
[考研] 295求調(diào)劑 +4	一志愿京區(qū)211 2026-03-18	6/300	2026-03-20 23:41 by JourneyLucky
[考研] 一志愿中海洋材料工程專(zhuān)碩330分求調(diào)劑 +8	小材化本科 2026-03-18	8/400	2026-03-20 23:16 by JourneyLucky
[論文投稿] 申請(qǐng)回稿延期一個(gè)月，編輯同意了。但系統(tǒng)上的時(shí)間沒(méi)變，給編輯又寫(xiě)郵件了，沒(méi)回復(fù) 10+3	wangf9518 2026-03-17	4/200	2026-03-19 23:55 by babero
[碩博家園] 湖北工業(yè)大學(xué) 生命科學(xué)與健康學(xué)院-課題組招收2026級(jí)食品/生物方向碩士 +3	1喜春8 2026-03-17	5/250	2026-03-17 17:18 by ber川cool子