大海的聲音

新蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 789
帖子: 168
在線: 172.1小時
蟲號: 2691385
注冊: 2013-09-30
性別: GG
專業(yè): 安全科學(xué)與工程

[求助] MATLAB求某一式子怎么求不出結(jié)果已有1人參與

如圖，
每個x對應(yīng)一個w，x有一列數(shù)，求其對應(yīng)的w；
x=[ ]; //數(shù)據(jù)輸入
c=[662.3012161 1.235047160711485 0.003402653539486 -0.421886852153956];
syms n;
u=symsum((-1)^n.*x.^(n*c(4)+2)/((c(4)*n+2)*c(3)^(n*c(4))*prod(n)),n,1,300);
w=c(1).*x.^2*(1-c(2))/2+c(1)*c(2).*x.^2/2+c(1)*c(2).*u-c(1).*x.^2.*(1-c(2)*(1-exp(-(x/c(3)).^c(4)))).^2/2;
plot(x,y,'*');
一直Busy，怎么沒有反應(yīng)，求各位指點迷津！

chuan.jpg

回復(fù)此樓

» 猜你喜歡

【求調(diào)劑】新能源材料本科，一志愿211，初試321 已經(jīng)有3人回復(fù)
材料考研調(diào)劑已經(jīng)有3人回復(fù)
材料調(diào)劑已經(jīng)有12人回復(fù)
英一數(shù)一408，總分284，二戰(zhàn)真誠求調(diào)劑已經(jīng)有14人回復(fù)
085410 一志愿211 22408分?jǐn)?shù)359求調(diào)劑已經(jīng)有4人回復(fù)
271求調(diào)劑已經(jīng)有19人回復(fù)
385分生物學(xué)（071000）求調(diào)劑已經(jīng)有3人回復(fù)
一志愿安徽大學(xué)計算機(jī)科學(xué)與技術(shù)學(xué)碩，331分求調(diào)劑已經(jīng)有3人回復(fù)
318求調(diào)劑，計算材料方向已經(jīng)有8人回復(fù)
291求調(diào)劑已經(jīng)有25人回復(fù)

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

求救，MATLAB程序顯示太大運行不了已經(jīng)有3人回復(fù)
怎樣用matlab輸出這個式子？已經(jīng)有11人回復(fù)
求高手給做一道題，用MATLAB或1stOpt，請把結(jié)果直接給出來，我的軟件就不行，謝謝已經(jīng)有5人回復(fù)
MATLAB無法計算出這個含【級數(shù)和超幾何函數(shù)】的式子已經(jīng)有6人回復(fù)
為什么matlab算出的數(shù)不正確，在線等已經(jīng)有5人回復(fù)
matlab “Inner matrix dimensions must agree.”一直不知道哪里出錯了求助已經(jīng)有6人回復(fù)
matlab求三元4次符號代數(shù)方程組，出現(xiàn)warning且答案不正確，該怎么辦？已經(jīng)有12人回復(fù)
有個matlab程序，但是運行不下來，求高手已經(jīng)有23人回復(fù)
使用拉格朗日優(yōu)化，但Matlab 運行不出結(jié)果來已經(jīng)有7人回復(fù)
Fortran 求解復(fù)數(shù)矩陣SVD時，用到cgesvd，結(jié)果與matlab不同，為什么？已經(jīng)有15人回復(fù)
對于符號積分matlab積分積不出來怎么解決+有代碼已經(jīng)有14人回復(fù)
用matlab運行代碼出不來圖像 �。� 已經(jīng)有3人回復(fù)
matlab解不等式組已經(jīng)有6人回復(fù)
編一個matlab程序，求某一向量在向量場里的坐標(biāo) 已經(jīng)有3人回復(fù)
求助，matlab公式繪圖，公式里面有兩個數(shù)組，編出小程序，但不出圖像已經(jīng)有4人回復(fù)
大家有沒有那種idea理論推導(dǎo)挺合理，但matlab仿真仿不出來的情況��？已經(jīng)有6人回復(fù)
matlab編程大家看這兩個式子該如何編寫呢？已經(jīng)有3人回復(fù)
matlab 積分計算算不出結(jié)果，求幫助！��！已經(jīng)有10人回復(fù)
如何用matlab求一個式子的整數(shù)解已經(jīng)有4人回復(fù)
matlab線性矩陣不等式求解做優(yōu)化問題已經(jīng)有10人回復(fù)
matlab 計算過程中工作空間變量保存及讀取求解決內(nèi)存不足的方法已經(jīng)有5人回復(fù)
MATLAB初學(xué)者，這樣的式子怎么就是畫不了圖，總是出錯!求指導(dǎo)!!!!!!!!!!!!! 已經(jīng)有8人回復(fù)
請幫忙用Matlab算一下附件中式子的結(jié)果，有急用，謝謝！已經(jīng)有6人回復(fù)
一個復(fù)雜式子的MATLAB解算已經(jīng)有8人回復(fù)
【求助】利用MATLAB求解的外圍輪廓包絡(luò)線為什么不與圖形外圍輪廓重合？怎樣才可以精確已經(jīng)有8人回復(fù)

1樓 2015-01-16 15:47:43

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

kabaisun

專家顧問 (正式寫手)

應(yīng)助: 20 (小學(xué)生)
金幣: 1311.9
散金: 681
紅花: 3
帖子: 315
在線: 113.8小時
蟲號: 1990141
注冊: 2012-09-10
性別: GG
專業(yè): 交通運輸經(jīng)濟(jì)學(xué)
管轄: 計算模擬區(qū)

引用回帖:

3樓: Originally posted by 大海的聲音 at 2015-01-16 19:35:37
謝謝，我這個x=由小到大無規(guī)律排列的一列數(shù)，有14000多個數(shù)，看來只能把n縮小點了，雖然這樣是有誤差的。...

x數(shù)據(jù)太多可以試試把u的計算式換成矩陣乘法，把syms符號計算去掉。

贊一下

回復(fù)此樓

嚴(yán)實誠勇

4樓2015-01-16 20:28:37

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

查看全部 10 個回答

kabaisun

專家顧問 (正式寫手)

應(yīng)助: 20 (小學(xué)生)
金幣: 1311.9
散金: 681
紅花: 3
帖子: 315
在線: 113.8小時
蟲號: 1990141
注冊: 2012-09-10
性別: GG
專業(yè): 交通運輸經(jīng)濟(jì)學(xué)
管轄: 計算模擬區(qū)

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

★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
大海的聲音: 金幣+4, ★★★很有幫助 2015-01-16 20:28:03

你這個式子求u時很費時，所以一直出不來結(jié)果，你要看怎么把計算u的時間縮短。如果x個數(shù)給的比較少的話，像取十個x就比較快；提升速度還可以把n調(diào)小些。

以下是修改的代碼，

CODE:

c=[662.3012161 1.235047160711485 0.003402653539486 -0.421886852153956];

x=1:10;

syms n ;

u=single(symsum((-1)^n.*x.^(n*c(4)+2)/((c(4)*n+2)*c(3)^(n*c(4))*prod(n)),n,1,100));

w=c(1).*x.^2*(1-c(2))/2+c(1)*c(2).*x.^2/2+c(1)*c(2).*u-c(1).*x.^2.*(1-c(2)*(1-exp(-(x/c(3)).^c(4)))).^2/2;

plot(x,w,'*');

贊一下

回復(fù)此樓

嚴(yán)實誠勇

2樓2015-01-16 19:28:31

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

大海的聲音

新蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 789
帖子: 168
在線: 172.1小時
蟲號: 2691385
注冊: 2013-09-30
性別: GG
專業(yè): 安全科學(xué)與工程

引用回帖:

2樓: Originally posted by kabaisun at 2015-01-16 19:28:31
你這個式子求u時很費時，所以一直出不來結(jié)果，你要看怎么把計算u的時間縮短。如果x個數(shù)給的比較少的話，像取十個x就比較快；提升速度還可以把n調(diào)小些。

以下是修改的代碼，

c=;
x=1:10;
syms n ;
u=single ...

謝謝，我這個x=[0.000006 0.000008 0.00001 0.000012 0.000014 0.00003 .....]由小到大無規(guī)律排列的一列數(shù)，有14000多個數(shù)，看來只能把n縮小點了，雖然這樣是有誤差的。

贊一下

回復(fù)此樓

3樓2015-01-16 19:35:37

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

大海的聲音

新蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 789
帖子: 168
在線: 172.1小時
蟲號: 2691385
注冊: 2013-09-30
性別: GG
專業(yè): 安全科學(xué)與工程

引用回帖:

4樓: Originally posted by kabaisun at 2015-01-16 20:28:37
x數(shù)據(jù)太多可以試試把u的計算式換成矩陣乘法，把syms符號計算去掉。...

這個不好換吧？我這里面u表示的是無窮級數(shù)求和。

贊一下

回復(fù)此樓

5樓2015-01-16 20:31:51

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

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

[考研] 英一數(shù)一408，總分284，二戰(zhàn)真誠求調(diào)劑 +12	12.27 2026-03-30	14/700	2026-04-02 00:18 by 欣喜777
[考研] 301求調(diào)劑 +12	A_JiXing 2026-04-01	12/600	2026-04-01 22:47 by peike
[考研] 材料專碩322分 +10	哈哈哈吼吼吼哈 2026-04-01	10/500	2026-04-01 22:19 by Dyhoer
[考研] 303分 0807學(xué)碩求調(diào)劑 +3	TYC3632 2026-04-01	3/150	2026-04-01 19:24 by lwk2004
[考研] 265求調(diào)劑 +11	yelck 2026-04-01	12/600	2026-04-01 19:12 by 549790059
[考研] 一志愿085600中科院寧波所276分求調(diào)劑 +22	材料學(xué)257求調(diào)劑 2026-03-28	23/1150	2026-04-01 19:03 by 客爾美德
[考研] 求調(diào)劑 +4	圖鑒212 2026-03-30	5/250	2026-04-01 15:32 by 圖鑒212
[考研] 086000生物與醫(yī)藥初試274求調(diào)劑 +5	小叮當(dāng)來了 2026-03-30	6/300	2026-04-01 14:25 by mxj1203
[考研] 291求調(diào)劑 +3	迷蒙木木 2026-04-01	4/200	2026-04-01 11:07 by 逆水乘風(fēng)
[考研] 一志愿武理材料工程302調(diào)劑環(huán)化或化工 +15	Doleres 2026-03-31	16/800	2026-04-01 09:49 by lfj11
[考研] 材料工程專碩求調(diào)劑 +10	hyl3153942 2026-03-29	10/500	2026-03-31 16:31 by hypershenger
[考研] 學(xué)碩274求調(diào)劑 +17	Li李魚 2026-03-26	17/850	2026-03-31 15:19 by 客爾美德
[考研] 0703化學(xué) +20	妮妮ninicgb 2026-03-27	20/1000	2026-03-31 13:33 by 無際的草原
[考研] 323分食品與營養(yǎng)調(diào)劑 +3	嘿ooo 2026-03-31	3/150	2026-03-31 09:38 by longlotian
[考研] 332求調(diào)劑 +6	@MZB382400 2026-03-28	6/300	2026-03-30 16:57 by 無際的草原
[考研] 348求調(diào)劑 +6	小懶蟲不懶了 2026-03-28	6/300	2026-03-30 10:29 by Evan_Liu
[考研] 11408軟件工程求調(diào)劑 +3	Qiu學(xué)ing 2026-03-28	3/150	2026-03-28 21:50 by zhq0425
[考研] 求調(diào)劑 +7	爭取九點睡 2026-03-28	8/400	2026-03-28 21:07 by 爭取九點睡
[考研] 本科新能源科學(xué)與工程，一志愿華理能動285求調(diào)劑 +7	AZMK 2026-03-28	11/550	2026-03-28 21:01 by xxxsssccc
[考研] 315分求調(diào)劑 +7	26考研上岸版26 2026-03-26	7/350	2026-03-28 04:05 by fmesaito