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

返回列表

上一頁(yè)

下一頁(yè)

蘋果核兒

金蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 381.5
紅花: 1
帖子: 171
在線: 44.2小時(shí)
蟲號(hào): 1398885
注冊(cè): 2011-09-12
專業(yè): 制造系統(tǒng)與自動(dòng)化

引用回帖:

10樓: Originally posted by 清晨的兔子 at 2012-06-04 08:56:55
不需要吧，那就在for i=length(T)后面加一句T=nt(i)，x=nx(i)...

那你這for程序的循環(huán)變量是哪個(gè)向量的下標(biāo)呢？

贊一下

回復(fù)此樓

把握當(dāng)下的時(shí)光，不要計(jì)較于自己一時(shí)的得與失�。�！

11樓2012-06-04 09:07:07

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

清晨的兔子

銀蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 272.1
散金: 60
紅花: 2
帖子: 135
在線: 39.6小時(shí)
蟲號(hào): 1283857
注冊(cè): 2011-05-03
性別: MM
專業(yè): 化工熱力學(xué)和基礎(chǔ)數(shù)據(jù)

引用回帖:

11樓: Originally posted by 蘋果核兒 at 2012-06-04 09:07:07
那你這for程序的循環(huán)變量是哪個(gè)向量的下標(biāo)呢？...

t=[0.07 0.13 0.28 0.44],那么t(1)=o.o7，t(2)=o.13，t(3)=o.28，t(4)=o.44
for i=1:length(T)，就相當(dāng)于for i=1:4，i=1的時(shí)候用對(duì)應(yīng)著的函數(shù)m,n,o,p等求出結(jié)果并進(jìn)行比較選擇，然后i=2時(shí)再進(jìn)行同樣的下一次的循環(huán)，其實(shí)就是求了4組數(shù)。

贊一下

回復(fù)此樓

做只勤勞的小螞蟻！

12樓2012-06-04 09:31:15

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

lgycjpcqu

金蟲 (正式寫手)

應(yīng)助: 35 (小學(xué)生)
金幣: 566.8
散金: 411
帖子: 371
在線: 176.9小時(shí)
蟲號(hào): 1327328
注冊(cè): 2011-06-20
性別: GG
專業(yè): 電化學(xué)

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

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

x=[0.07 0.13 0.28 0.44];
T=[283 296 312 320];
m=x.^2+x.*2;
n=2*x.^2+T.*3;
o=4*x.^2+5*T.^2;
p=m+n.*o;
q=3*(T.^2).^(-1);
qp=q-p;
y=exp(qp);
k=y-x;
L=length(x);
xcal=zeros(1,L);
for ii=1:L
if abs(k(ii))<1*10^(-5)
xcal(ii)=y(ii);
else
xcal(ii)=x(ii);
end
end
xcal

贊一下

回復(fù)此樓

13樓2012-06-04 09:53:54

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

csgt0

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

彩色掛圖

專家經(jīng)驗(yàn): +2
應(yīng)助: 367 (碩士)
貴賓: 0.457
金幣: 7807
散金: 8897
紅花: 20
沙發(fā): 7
帖子: 2809
在線: 515.8小時(shí)
蟲號(hào): 449046
注冊(cè): 2007-11-02
專業(yè): 金屬非晶態(tài)、準(zhǔn)晶和納米晶
管轄: 計(jì)算模擬

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

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

小妹子啊，你這還用解嗎？一眼就能看出q《《《《p，y等于0，你要abs(y-x)<=0.00001,只要取x<=0.00001就可以了。

贊一下

回復(fù)此樓

showmethemoney

14樓2012-06-04 09:59:04

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

清晨的兔子

銀蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 272.1
散金: 60
紅花: 2
帖子: 135
在線: 39.6小時(shí)
蟲號(hào): 1283857
注冊(cè): 2011-05-03
性別: MM
專業(yè): 化工熱力學(xué)和基礎(chǔ)數(shù)據(jù)

引用回帖:

14樓: Originally posted by csgt0 at 2012-06-04 09:59:04
小妹子啊，你這還用解嗎？一眼就能看出q《《《《p，y等于0，你要abs(y-x)<=0.00001,只要取x<=0.00001就可以了。

(⊙o⊙)…，那個(gè)啥啥，函數(shù)文件時(shí)隨便寫的啦，我只是想要程序啦

贊一下

回復(fù)此樓

做只勤勞的小螞蟻！

15樓2012-06-04 10:15:09

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

清晨的兔子

銀蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 272.1
散金: 60
紅花: 2
帖子: 135
在線: 39.6小時(shí)
蟲號(hào): 1283857
注冊(cè): 2011-05-03
性別: MM
專業(yè): 化工熱力學(xué)和基礎(chǔ)數(shù)據(jù)

引用回帖:

13樓: Originally posted by lgycjpcqu at 2012-06-04 09:53:54
x=;
T=;
m=x.^2+x.*2;
n=2*x.^2+T.*3;
o=4*x.^2+5*T.^2;
p=m+n.*o;
q=3*(T.^2).^(-1);
qp=q-p;
y=exp(qp);
k=y-x;
L=length(x);
xcal=zeros(1,L);
for ii=1:L
if abs(k(ii))<1*10^(-5)
...

帶進(jìn)去試過了，這個(gè)計(jì)算出來(lái)的結(jié)果跟x=[0.07 0.13 0.28 0.44]是一樣的啊，程序不對(duì)

贊一下

回復(fù)此樓

做只勤勞的小螞蟻！

16樓2012-06-04 10:17:27

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

csgt0

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

彩色掛圖

專家經(jīng)驗(yàn): +2
應(yīng)助: 367 (碩士)
貴賓: 0.457
金幣: 7807
散金: 8897
紅花: 20
沙發(fā): 7
帖子: 2809
在線: 515.8小時(shí)
蟲號(hào): 449046
注冊(cè): 2007-11-02
專業(yè): 金屬非晶態(tài)、準(zhǔn)晶和納米晶
管轄: 計(jì)算模擬

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

引用回帖:

15樓: Originally posted by 清晨的兔子 at 2012-06-04 10:15:09
(⊙o⊙)…，那個(gè)啥啥，函數(shù)文件時(shí)隨便寫的啦，我只是想要程序啦...

數(shù)據(jù)會(huì)影響程序的，估計(jì)你要的是這個(gè)

clc
clear
x0=[0.07 0.13 0.28 0.44];
T=[283 296 312 320];
for i=1:length(x0)
k(i)=1;
y(i)=x0(i);
while k(i)>0.00001
      x(i)=y(i);
      m=x(i)^2+2*x(i);
      n=2*x(i)^2+3*T(i);
      o=4*x(i)^2+5*T(i)^2;
      p=m+n.*o;
      q=3./(T(i)^2);
      y(i)=exp(q-p);
      k(i)=abs(y(i)-x(i));
end
end
x
k
y

但是用這個(gè)也行
%需要輸入的數(shù)組為x數(shù)組和T數(shù)組
function main
clc
clear
x0=[0.07 0.13 0.28 0.44];
options = optimset('Display','iter','TolFun',1e-5);
x=fsolve(@calk,x0,options)
end

%計(jì)算用到的函數(shù)如下，這是我隨便寫的函數(shù)，主要寫明每個(gè)函數(shù)里面的未知輸入量是哪個(gè)
function k=calk(x)
T=[283 296 312 320];
m=x.^2+2*x;
n=2*x.^2+3*T;
o=4*x.^2+5*T.^2;
p=m+n.*o;
q=3./(T.^2);
%我要得到的結(jié)果如下
y=exp(q-p);
%最終目的是將計(jì)算得到的y(i)與代入的初值x(i)比較
k=abs(y-x)-0.00001;
end
%如果k的絕對(duì)值abs(k)<=0.00001，則最終的輸出結(jié)果為xcal=y(i)，否則將得出的y(i)賦值給最初代入計(jì)算的x，接著進(jìn)行計(jì)算，直到abs(k)<=0.00001為止。

贊一下(2人)

回復(fù)此樓

showmethemoney

17樓2012-06-04 10:27:39

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

disini

新蟲 (初入文壇)

應(yīng)助: 1 (幼兒園)
金幣: 25.5
帖子: 2
在線: 28分鐘
蟲號(hào): 1678865
注冊(cè): 2012-03-09

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

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
清晨的兔子: 金幣+30, ★★★很有幫助 2012-06-04 14:58:52

x=[0.07 0.13 0.28 0.44]
T=[283 296 312 320]
m=x.^2+2*x
n=2*x.^2+3*T
o=4*x.^2+5*T.^2
p=m+n.*o
q=3*(T.^(-2))
for i=1:length(T);
  y(i)=exp(q(i)-p(i));
  k(i)=y(i)-x(i);
   if abs(k)<=0.00001
   xcal=y(i);
   end
  end

好了。
最后的y(i)全部是0，因?yàn)槟氵x的數(shù)值不合適，xcal沒有值，也是這個(gè)原因，不符合if設(shè)定的條件。

親你編程序要嚴(yán)謹(jǐn)啊，要認(rèn)真點(diǎn)寫才行，語(yǔ)法，邏輯都不能有一點(diǎn)錯(cuò)誤
比如說那個(gè).*,.^
這都不能漏。。就算舉例子也得仔細(xì)些。
這又不是寫小說

贊一下

回復(fù)此樓

18樓2012-06-04 12:55:08

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

gov_sb

新蟲 (初入文壇)

應(yīng)助: 2 (幼兒園)
金幣: 0.8
紅花: 1
帖子: 5
在線: 35分鐘
蟲號(hào): 1846519
注冊(cè): 2012-06-04

x=[1 2];
T=[10 20];
i = 1;
while i<=length(x)
  m=2*x(i)+3*T(i);
  n=2*x(i).^2+3*T(i).^2;
  y(i)=n./10-m;
  k=y(i)-x(i);
  if k<-2
i = i+1
  else
x(i) = y(i)
  end
end

樓主要的是這個(gè)吧，16樓的雖然也能行，
但是時(shí)間復(fù)雜度和可讀性有待改進(jìn)

贊一下(3人)

回復(fù)此樓

19樓2012-06-04 19:38:33

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

gov_sb

新蟲 (初入文壇)

應(yīng)助: 2 (幼兒園)
金幣: 0.8
紅花: 1
帖子: 5
在線: 35分鐘
蟲號(hào): 1846519
注冊(cè): 2012-06-04

忘記寫輸出了，嘎嘎
.........
if k<-2
i = i+1
在這里輸出k就是了
else
.......

贊一下(1人)

回復(fù)此樓

20樓2012-06-04 19:40:40

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

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

返回列表

上一頁(yè)

下一頁(yè)

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

[考研] 材料專碩306英一數(shù)二 +10	z1z2z3879 2026-03-16	13/650	2026-03-18 14:20 by 007_lilei
[考研] 070300化學(xué)319求調(diào)劑 +6	錦鯉0909 2026-03-17	6/300	2026-03-18 13:22 by Iveryant
[考研] 266求調(diào)劑 +4	陽(yáng)陽(yáng)哇塞 2026-03-14	8/400	2026-03-18 12:43 by Linda Hu
[考研] 303求調(diào)劑 +4	睿08 2026-03-17	6/300	2026-03-18 11:01 by Iveryant
[考博] 26博士申請(qǐng) +3	1042136743 2026-03-17	3/150	2026-03-17 23:30 by 輕松不少隨
[考研] 085600材料與化工求調(diào)劑 +5	緒幸與子 2026-03-17	5/250	2026-03-17 16:40 by laoshidan
[考研] 26考研求調(diào)劑 +6	丶宏Sir 2026-03-13	6/300	2026-03-17 16:13 by 醉在風(fēng)里
[考研] 材料與化工專碩調(diào)劑 +5	heming3743 2026-03-16	5/250	2026-03-17 14:03 by 勇敢太監(jiān)王公公
[考研] 302求調(diào)劑 +4	小賈同學(xué)123 2026-03-15	8/400	2026-03-17 10:33 by 小賈同學(xué)123
[考研] 東南大學(xué)364求調(diào)劑 +5	JasonYuiui 2026-03-15	5/250	2026-03-16 21:28 by 木瓜膏
[考研] 277材料科學(xué)與工程080500求調(diào)劑 +3	自由煎餅果子 2026-03-16	3/150	2026-03-16 14:10 by 運(yùn)氣yunqi
[考研] 本科南京大學(xué)一志愿川大藥學(xué)327 +3	麥田耕者 2026-03-14	3/150	2026-03-14 20:04 by 外星文明
[考研] 265求調(diào)劑 +4	威化餅07 2026-03-12	4/200	2026-03-14 17:23 by userper
[基金申請(qǐng)] 現(xiàn)在如何回避去年的某一個(gè)專家，不知道名字 +3	zk200107 2026-03-12	6/300	2026-03-14 17:13 by zk200107
[考研] 301求調(diào)劑 +6	Liyouyumairs 2026-03-11	6/300	2026-03-13 20:11 by JourneyLucky
[考研] 求b區(qū)學(xué)校調(diào)劑 +3	周56 2026-03-11	3/150	2026-03-13 16:20 by JourneyLucky
[考研] 277求調(diào)劑 +4	anchor17 2026-03-12	4/200	2026-03-13 11:15 by 白夜悠長(zhǎng)
[考研] 0817化學(xué)工程與技術(shù)考研312分調(diào)劑 +3	T123 tt 2026-03-12	3/150	2026-03-13 10:49 by houyaoxu
[考研] 290求調(diào)劑 +3	ADT 2026-03-13	3/150	2026-03-13 10:19 by peike
[考研] 283求調(diào)劑，材料、化工皆可 +8	蘇打水7777 2026-03-11	10/500	2026-03-13 09:06 by Linda Hu