| 9 | 1/1 | 返回列表 |
| 查看: 1566 | 回復(fù): 8 | ||
foreverolia新蟲 (初入文壇)
|
[求助]
matlab求解非線性方程組報(bào)錯,請各位大神指點(diǎn)! 已有1人參與
|
|
clear >> syms n beta; >> f1=sym(\'(f01-4*(2-n)*n*f01^(n-1)*sigma1^2*beta*yita1)/(1+2*n*(n-1)*f01^(n-2)*sigma1^2*beta*yita1)=fc1\') f1 = (f01 + 4*f01^(n - 1)*n*sigma1^2*yita1*(n - 2)*beta)/(2*f01^(n - 2)*n*sigma1^2*yita1*(n - 1)*beta + 1) = fc1 >> f2=sym(\'(f02-4*(2-n)*n*f02^(n-1)*sigma2^2*beta*yita2)/(1+2*n*(n-1)*f02^(n-2)*sigma2^2*beta*yita2)=fc2\') f2 = (f02 + 4*f02^(n - 1)*n*sigma2^2*yita2*(n - 2)*beta)/(2*f02^(n - 2)*n*sigma2^2*yita2*(n - 1)*beta + 1) = fc2 >> fo1=4.92; >> fc1=3.76; >> fo2=2.45; >> fc2=2.31; >> sigma1=1.50; >> sigma2=0.635; >> [n,beta]=solve(f1,f2) Warning: Explicit solution could not be found. > In solve at 81 n = [ empty sym ] beta = [] f1和f2是所求的方程組,算法總是報(bào)錯,不知道該如何求解。。。 |
新蟲 (初入文壇)
木蟲 (正式寫手)

新蟲 (初入文壇)
木蟲 (正式寫手)
|
為圖方便,給你用1stOpt編程算了下。 代碼: 采用L-M+UGO結(jié)果: 可能結(jié)果1: n: 1 betta: 0.0879576864290957 yita1: 1.46535105823627 yita2: 0.986840117382106 可能結(jié)果2: n: 1 betta: 0.235202634137847 yita1: 0.547990839309009 yita2: 0.369044224009309 上述結(jié)果都不一定正確,因?yàn)槲也恢滥愕淖兞康娜≈捣秶,而且是否是連續(xù)或者整數(shù)變量,我假設(shè)n為非負(fù)整數(shù)變量,其余均為實(shí)數(shù)域連續(xù)變量了,你根據(jù)你的問題修改下吧,也可以采取其它算法試一下。 |

新蟲 (初入文壇)
|
高手,想跟您請教一下,plot函數(shù)畫圖的時候,為什么x=1:1:10和x=1:0.1:10得到的函數(shù)圖像走勢不同?應(yīng)該是一樣的呀。。。程序代碼如下:clear all; >> n=1:1:10; >> y=(4.92-3.72)/(2.45-2.34)-((1.5/0.635)^2)*((4.92/2.45).^(n-2)).*(((n-1)*3.72+(4-2*n)*4.92)./((n-1).*2.34+(4-2*n)*2.45)); >> plot(n,y) 以及 >> clear all; >> n=1:0.1:10; >> y=(4.92-3.72)/(2.45-2.34)-((1.5/0.635)^2)*((4.92/2.45).^(n-2)).*(((n-1)*3.72+(4-2*n)*4.92)./((n-1).*2.34+(4-2*n)*2.45)); >> plot(n,y) 兩種情況的函數(shù)圖像不同,您能給指點(diǎn)一下嗎? |
木蟲 (正式寫手)

新蟲 (初入文壇)
木蟲 (正式寫手)

| 9 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 北科281學(xué)碩材料求調(diào)劑 +16 | tcxiaoxx 2026-03-20 | 18/900 |
|
|---|---|---|---|---|
|
[考研] 一志愿哈工大,085400,320,求調(diào)劑 +4 | gdlf9999 2026-03-24 | 4/200 |
|
|
[考研] 332求調(diào)劑 +6 | 032500 2026-03-25 | 6/300 |
|
|
[考研] 299求調(diào)劑 +7 | 某某某某位 2026-03-21 | 8/400 |
|
|
[考研] 材料與化工304求B區(qū)調(diào)劑 +3 | 邱gl 2026-03-25 | 3/150 |
|
|
[考研] 生物學(xué)學(xué)碩求調(diào)劑 +7 | 小羊睡著了? 2026-03-23 | 10/500 |
|
|
[考研] 一志愿武理085500機(jī)械專業(yè)總分300求調(diào)劑 +3 | an10101 2026-03-24 | 7/350 |
|
|
[考研] 材料專碩331求調(diào)劑 +4 | 鮮當(dāng)牛 2026-03-24 | 4/200 |
|
|
[考研] 276求調(diào)劑。有半年電池和半年高分子實(shí)習(xí)經(jīng)歷 +9 | 材料學(xué)257求調(diào)劑 2026-03-23 | 10/500 |
|
|
[考研] 384求調(diào)劑 +3 | 子系博 2026-03-22 | 6/300 |
|
|
[考研] 工科0856求調(diào)劑 +5 | 沐析汀汀 2026-03-21 | 5/250 |
|
|
[考研] 接收2026碩士調(diào)劑(學(xué)碩+專碩) +4 | allen-yin 2026-03-23 | 6/300 |
|
|
[考研] 306求調(diào)劑 +5 | 來好運(yùn)來來來 2026-03-22 | 5/250 |
|
|
[考研] 材料學(xué)碩301分求調(diào)劑 +7 | Liyouyumairs 2026-03-21 | 7/350 |
|
|
[考研] 初試 317 +7 | 半拉月丙 2026-03-20 | 7/350 |
|
|
[考研]
|
Grand777 2026-03-21 | 3/150 |
|
|
[考研] 336求調(diào)劑 +5 | rmc8866 2026-03-21 | 5/250 |
|
|
[考研] 085601調(diào)劑 358分 +3 | zzzzggh 2026-03-20 | 4/200 |
|
|
[考研] 一志愿 南京航空航天大學(xué)大學(xué) ,080500材料科學(xué)與工程學(xué)碩 +5 | @taotao 2026-03-20 | 5/250 |
|
|
[考研] 一志愿西安交通大學(xué) 學(xué)碩 354求調(diào)劑211或者雙一流 +3 | 我想要讀研究生 2026-03-20 | 3/150 |
|