chuzhenhui

金蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 600.7
散金: 35
帖子: 129
在線: 10.2小時(shí)
蟲號(hào): 763256
注冊(cè): 2009-05-05
性別: GG
專業(yè): 食品&包裝

[交流] 【求助】請(qǐng)教一個(gè)關(guān)于解方程的問題已有3人參與

問題描述：

在做一個(gè)gui程序，需要根據(jù)輸入的數(shù)值，進(jìn)行相應(yīng)的計(jì)算。比如這樣一個(gè)方程：

y=(w*c*k*x)/((1-k*x)*(1-k*x+c*k*x))

其中w，c，k，y為常數(shù)，x是變量。程序需要根據(jù)w，c，k，y的輸入，代入方程中計(jì)算x的值，后面程序再用x進(jìn)行其他的運(yùn)算。

假如輸入的常數(shù)分別是:w=0.0378;c=9.992;k=0.9247;y=0.026，我這么編：

w=0.0378;c=9.992;k=0.9247;y=0.026;

x=solve('(w*c*k*x)/((1-k*x)*(1-k*x+c*k*x))-y');

-------------------------
matlab出現(xiàn)錯(cuò)誤如下：

??? Error using ==> solve>getEqns at 182
' w1*c1*k1*awi)/((1-k1*awi)*(1-k1*awi+c1*k1*awi))=in ' is not a valid expression or equation.

Error in ==> solve at 67
[eqns,vars] = getEqns(varargin{:});

Error in ==> solvetest at 3
solve('w1*c1*k1*awi)/((1-k1*awi)*(1-k1*awi+c1*k1*awi))=in');
-----------------------------------------------------------------------------------
直接用輸入的數(shù)字代入方程，用solve是可以算出來(lái)的，x值有兩個(gè)解，我要取正解。但是像以上的程序，先對(duì)符號(hào)賦值，再用符號(hào)代入方程中用solve就不能用了。

我想解決的就是如果要算以上方程中的x,我應(yīng)該用什么方法？這樣用solve和fsolve都不行。望大俠賜教！謝謝~~~

回復(fù)此樓

» 猜你喜歡

請(qǐng)問共同通訊和共同一作的認(rèn)可度問題已經(jīng)有4人回復(fù)
286分調(diào)劑已經(jīng)有22人回復(fù)
材料調(diào)劑已經(jīng)有4人回復(fù)
085600 295分求調(diào)劑已經(jīng)有19人回復(fù)
材料調(diào)劑已經(jīng)有4人回復(fù)
302求調(diào)劑一志愿華中師范大學(xué) 已經(jīng)有3人回復(fù)
環(huán)境285分，過(guò)六級(jí)，求調(diào)劑已經(jīng)有5人回復(fù)
08工科求調(diào)劑290分已經(jīng)有8人回復(fù)
321求調(diào)劑一志愿浙江工業(yè)大學(xué)生物醫(yī)藥已經(jīng)有4人回復(fù)
266分，一志愿電氣工程，本科材料，求材料專業(yè)調(diào)劑已經(jīng)有4人回復(fù)

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

請(qǐng)教一個(gè)關(guān)于吸附的問題已經(jīng)有9人回復(fù)
求matlab求解一個(gè)方程組的問題已經(jīng)有6人回復(fù)
【求助】向各位大俠求助matlab求解微分方程組遇到的一個(gè)問題已經(jīng)有21人回復(fù)
【求助】請(qǐng)教一個(gè)關(guān)于Mathematica中計(jì)算小量的問題已經(jīng)有4人回復(fù)
【討論】請(qǐng)教一個(gè)氣-液-固三相流模擬可能性的問題已經(jīng)有6人回復(fù)
【求助】偏微分方程解法已經(jīng)有19人回復(fù)

找到自己。

1樓 2010-10-11 19:52:22

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

robert2020

銀蟲 (著名寫手)

仿真EPI: 2
應(yīng)助: 0 (幼兒園)
貴賓: 1.307
金幣: 292.3
散金: 9143
紅花: 19
沙發(fā): 6
帖子: 1950
在線: 620.1小時(shí)
蟲號(hào): 861676
注冊(cè): 2009-10-02

★ ★
ctgu_zheng(金幣+2):辛苦了。。。 2010-10-12 07:52:23
chuzhenhui(金幣+6):多謝��！ 2010-10-12 15:04:01

求解程序：

引用回帖:

clear all
syms x
w=0.0378;c=9.992;k=0.9247;y=0.026;
f=(w*c*k*x)/((1-k*x)*(1-k*x+c*k*x))-y;
x=vpa(solve(f))

驗(yàn)證程序：

引用回帖:

solve('0.026-(0.0378*9.992*0.9247*x)/((1-0.9247*x)*(1-0.9247*x+9.992*0.9247*x))=0')

結(jié)果：

引用回帖:

x =

-.926322075415857
.140404159021461

贊一下(1人)

回復(fù)此樓

Truthfulness,Compassion,Tolerance,Beauty

2樓2010-10-11 22:18:04

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

lijinfeng042

木蟲 (小有名氣)

Matlab

仿真EPI: 2
應(yīng)助: 1 (幼兒園)
金幣: 2156.1
散金: 115
帖子: 291
在線: 31.5小時(shí)
蟲號(hào): 1019062
注冊(cè): 2010-05-15
性別: GG
專業(yè): 通信理論與系統(tǒng)

★
hiqun(金幣+1):感謝應(yīng)助，辛苦了 2010-10-12 13:58:42
chuzhenhui(金幣+4):多謝��！ 2010-10-12 15:04:13

CODE:

>> [result,status] = evalin(symengine,'w:=0.0378;c:=9.992;k:=0.9247;y:=0.026;solve((w*c*k*x)/((1-k*x)*(1-k*x+c*k*x))-y)')

result =

[ [x = 0.14040415902146110092877008088218], [x = -0.92632207541585789273828437767426]]

status =

     0

方法還有不少符號(hào)計(jì)算subs代換也可以 vpa的內(nèi)部也是調(diào)用了計(jì)算引擎mupadengine

贊一下(1人)

回復(fù)此樓

工作了,偶爾會(huì)上來(lái)~可以關(guān)注新浪微博 @云是風(fēng)的夢(mèng)_Matlab

3樓2010-10-12 13:50:45

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

chuzhenhui

金蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 600.7
散金: 35
帖子: 129
在線: 10.2小時(shí)
蟲號(hào): 763256
注冊(cè): 2009-05-05
性別: GG
專業(yè): 食品&包裝

引用回帖:

Originally posted by robert2020 at 2010-10-11 22:18:04:
求解程序：

驗(yàn)證程序：

結(jié)果：

多謝指教，如果我要取x>0 的解，怎么弄？有專門設(shè)定方程解范圍的函數(shù)可以用嗎。。直接設(shè)定x>0好像不可以，我用y=max(x)也不行，請(qǐng)?jiān)僦笇?dǎo)一下吧，謝謝！

贊一下

回復(fù)此樓

找到自己。

4樓2010-10-12 17:00:05

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

chuzhenhui

金蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 600.7
散金: 35
帖子: 129
在線: 10.2小時(shí)
蟲號(hào): 763256
注冊(cè): 2009-05-05
性別: GG
專業(yè): 食品&包裝

引用回帖:

Originally posted by lijinfeng042 at 2010-10-12 13:50:45:
[code]

>> [result,status] = evalin(symengine,'w:=0.0378;c:=9.992;k:=0.9247;y:=0.026;solve((w*c*k*x)/((1-k*x)*(1-k*x+c*k*x))-y)')

result =

[ [x = 0.14040415902146110092877008088218] ...

謝謝指導(dǎo)，這個(gè)方程解出來(lái)的x值是中間值，我后面的程序中還要調(diào)用x值，我要的是x的正解，怎么弄啊，想了一些辦法，行不通。。。希望大俠再指導(dǎo)下！謝謝

[ Last edited by chuzhenhui on 2010-10-12 at 17:07 ]

贊一下

回復(fù)此樓

找到自己。

5樓2010-10-12 17:04:33

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

chuzhenhui

金蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 600.7
散金: 35
帖子: 129
在線: 10.2小時(shí)
蟲號(hào): 763256
注冊(cè): 2009-05-05
性別: GG
專業(yè): 食品&包裝

★
nono2009(金幣+1):鼓勵(lì)分享。 2010-10-13 06:20:34

我暫時(shí)解決了，主要是因?yàn)閤是sym型數(shù)據(jù)，我整成double型，再用 y=max(x)就可以得到x>0的解了，但這個(gè)方法有點(diǎn)投機(jī)，如果x有兩個(gè)以上大于0的解，這樣就行不通了，所以還是要請(qǐng)教正統(tǒng)的解決方法。

贊一下(1人)

回復(fù)此樓

找到自己。

6樓2010-10-12 17:16:29

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

robert2020

銀蟲 (著名寫手)

★
小木蟲(金幣+0.5):給個(gè)紅包，謝謝回帖交流

引用回帖:

x=[-1 -2 -3 4 2 4 0];
m=find(x>0);
x(m)

7樓2010-10-12 17:38:43

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

lijinfeng042

木蟲 (小有名氣)

Matlab

仿真EPI: 2
應(yīng)助: 1 (幼兒園)
金幣: 2156.1
散金: 115
帖子: 291
在線: 31.5小時(shí)
蟲號(hào): 1019062
注冊(cè): 2010-05-15
性別: GG
專業(yè): 通信理論與系統(tǒng)

★ ★
小木蟲(金幣+0.5):給個(gè)紅包，謝謝回帖交流
hiqun(金幣+1):感謝應(yīng)助 2010-10-13 16:18:48

引用回帖:

Originally posted by chuzhenhui at 2010-10-12 17:16:29:
我暫時(shí)解決了，主要是因?yàn)閤是sym型數(shù)據(jù)，我整成double型，再用 y=max(x)就可以得到x>0的解了，但這個(gè)方法有點(diǎn)投機(jī)，如果x有兩個(gè)以上大于0的解，這樣就行不通了，所以還是要請(qǐng)教正統(tǒng)的解決方法。

就是直接加個(gè)>0的條件 assuming>0

CODE:

>> result= evalin(symengine,'w:=0.0378;c:=9.992;k:=0.9247;y:=0.026;solve((w*c*k*x)/((1-k*x)*(1-k*x+c*k*x))-y) assuming x>0') result = [x = 0.14040415902146110092877008088218]

贊一下(2人)

回復(fù)此樓

工作了,偶爾會(huì)上來(lái)~可以關(guān)注新浪微博 @云是風(fēng)的夢(mèng)_Matlab

8樓2010-10-13 14:16:08

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

chuzhenhui

金蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 600.7
散金: 35
帖子: 129
在線: 10.2小時(shí)
蟲號(hào): 763256
注冊(cè): 2009-05-05
性別: GG
專業(yè): 食品&包裝

就是直接加個(gè)>0的條件 assuming>0

[code]>> result= evalin(symengine,'w:=0.0378;c:=9.992;k:=0.9247;y:=0.026;solve((w*c*k*x)/((1-k*x)*(1-k*x+c*k*x))-y) assuming x>0') result = [x ... [/quote]

非常謝謝，又多了一些知識(shí)，關(guān)于這塊我自己也再去搗鼓搗鼓！

贊一下

回復(fù)此樓

找到自己。

9樓2010-10-18 15:18:24

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

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

[考研] 324求調(diào)劑 +5	想上學(xué)求調(diào) 2026-04-02	5/250	2026-04-02 12:59 by wxiongid
[考研] 材料化工340求調(diào)劑 +5	jhx777 2026-03-30	5/250	2026-04-02 12:45 by smileboy2006
[考研] 285求調(diào)劑 +11	AZMK 2026-04-02	11/550	2026-04-02 12:07 by 1753564080
[考研] 292分，材料與化工，申請(qǐng)調(diào)劑 +19	程晴之 2026-04-01	23/1150	2026-04-02 11:59 by 程晴之
[考研] 一志愿鄭大材料工程290求調(diào)劑 +19	Youth_ 2026-03-30	19/950	2026-04-02 11:21 by Sammy2
[考研] 348求調(diào)劑 +6	小懶蟲不懶了 2026-03-27	7/350	2026-04-02 11:18 by sanrepian
[考研] 070305高分子化學(xué)與物理 304分求調(diào)劑 +14	c297914 2026-03-28	14/700	2026-04-02 09:35 by 星空星月
[考研] 343求調(diào)劑085601 +4	要努力學(xué)習(xí)x 2026-03-29	5/250	2026-04-02 09:01 by xuhui0306
[考研] 08生物與醫(yī)藥專碩初試346找調(diào)劑 +6	dianeeee 2026-04-01	7/350	2026-04-02 08:23 by guoweigw
[考研] 286求調(diào)劑 +5	Sa67890. 2026-04-01	7/350	2026-04-01 19:50 by 6781022
[考研] 【求調(diào)劑】085601材料工程專碩 \| 總分272 \| +10	腳滑的守法公民 2026-03-27	10/500	2026-04-01 17:23 by pies112
[考研] 一志愿北交材料工程總分358 +5	cs0106 2026-04-01	7/350	2026-04-01 11:45 by wangjy2002
[考研] 0855機(jī)械初試280求調(diào)劑 +3	kazenotori 2026-03-31	3/150	2026-04-01 10:08 by JourneyLucky
[考研] 299求調(diào)劑 +8	嗯嗯嗯嗯2 2026-03-27	8/400	2026-03-31 18:23 by lizhi8172
[考研] 315求調(diào)劑 +6	akie... 2026-03-28	7/350	2026-03-31 16:48 by asdfzly
[考研] 293分求調(diào)劑，外語(yǔ)為俄語(yǔ) +5	加一一九 2026-03-31	5/250	2026-03-31 09:39 by zhshch
[考研] 370求調(diào)劑 +3	080700調(diào)劑 2026-03-30	3/150	2026-03-31 01:09 by A_Zhe
[有機(jī)交流] 甲基亞磺磺酸鈉和甲基磺酸酯反應(yīng)機(jī)理 10+3	kaobao456 2026-03-29	4/200	2026-03-30 23:16 by nBu鋰
[考研] 11408軟件工程求調(diào)劑 +3	Qiu學(xué)ing 2026-03-28	3/150	2026-03-28 21:50 by zhq0425
[考研] 316求調(diào)劑 +7	江辭666 2026-03-26	7/350	2026-03-28 21:28 by sanrepian