snowwhite1

木蟲 (小有名氣)

應(yīng)助: 5 (幼兒園)
金幣: 4484.2
帖子: 275
在線: 104.1小時
蟲號: 1182357
注冊: 2011-01-03
性別: GG
專業(yè): 數(shù)論

[求助] 如何使用matlab求解非線性方程組的所有整數(shù)解？

例如：下面的方程組
x1^2=x1
x2^2=x2
x3^2=x3
x4^2=x4
x1^2=x1x2x3x4
使用fsolve需要預(yù)設(shè)初值，并且只能求出一個解，有沒有什么較好的方法能求出所有解？（這里所有的解都只能取值0或1）。

回復此樓

» 猜你喜歡

281求調(diào)劑（0805）已經(jīng)有8人回復
環(huán)境領(lǐng)域全國重點實驗室招收博士1-2名已經(jīng)有3人回復
材料專碩306英一數(shù)二已經(jīng)有10人回復
301求調(diào)劑已經(jīng)有6人回復
一志愿天津大學化學工藝專業(yè)（081702）315分求調(diào)劑已經(jīng)有7人回復
302求調(diào)劑已經(jīng)有6人回復
26博士申請已經(jīng)有3人回復
268求調(diào)劑已經(jīng)有3人回復
311求調(diào)劑已經(jīng)有10人回復
被我言中：新模板不強調(diào)格式了，假專家開始管格式了已經(jīng)有4人回復

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

求助解復雜非線性方程組的好的方法已經(jīng)有24人回復
請教一個matlab求解非線性方程組的問題已經(jīng)有9人回復
用matlab求解非線性方程組說無解，一定是方程組本身無解，還是有可能程序有問題呢？已經(jīng)有11人回復
線性方程有倆條斜率？已經(jīng)有3人回復
MATLAB求解非線性方程組已經(jīng)有5人回復
求Matlab解非線性方程代碼已經(jīng)有12人回復
第5章_解線性方程組的直接方法.ppt 已經(jīng)有13人回復
求Matlab高手解決線性方程組的迭代求解問題已經(jīng)有17人回復
使用fsolve求解非線性方程問題已經(jīng)有9人回復
matlab求解非線性方程組，錯誤提示怎么解決已經(jīng)有5人回復
請教用matlab求解一個非線性偏微分方程組的數(shù)值解已經(jīng)有8人回復
求助matlab---fsolve解非線性方程組已經(jīng)有6人回復
關(guān)于matlab線性規(guī)劃背包問題的解法求救呀~~方程知道但不知道如何解答···· 已經(jīng)有4人回復
1stopt或matlab如何求解以下的非線性方程并擬合出相應(yīng)曲線？已經(jīng)有13人回復
matlab求解非線性方程組已經(jīng)有16人回復
求高人指點用matlab求解非線性方程組，解決了追加100金幣；已經(jīng)有11人回復
MATLAB非線性方程組求解~~~反應(yīng)釜模擬已經(jīng)有9人回復
【求助】matlab求解非線性方程組，并畫圖處理。要求y,z是實數(shù)解！已經(jīng)有18人回復
【求助】不定方程：求x^2+y^2=2009的整數(shù)解已經(jīng)有3人回復
【討論】滑模、自抗擾、反饋線性化這三種控制方法的關(guān)系及前景已經(jīng)有14人回復
【求助】求教matlab解非線性方程組已經(jīng)有9人回復

奮斗

1樓 2012-03-21 18:18:28

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

sogooda

金蟲 (正式寫手)

應(yīng)助: 19 (小學生)
金幣: 1497.4
紅花: 2
帖子: 737
在線: 94.6小時
蟲號: 438143
注冊: 2007-09-16
專業(yè): 防災(zāi)工程

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

★
感謝參與，應(yīng)助指數(shù) +1
臭水溝: 金幣+1, 謝謝交流~~ 2012-03-26 10:32:27

一共就4個數(shù)，每個數(shù)還都只能取0和1兩種情況，那就窮舉法就行了。
換句話說，編個4重循環(huán)，一下子就出來了。

贊一下

回復此樓

2樓2012-03-25 21:26:47

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

sogooda

金蟲 (正式寫手)

應(yīng)助: 19 (小學生)
金幣: 1497.4
紅花: 2
帖子: 737
在線: 94.6小時
蟲號: 438143
注冊: 2007-09-16
專業(yè): 防災(zāi)工程

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

★ ★
臭水溝: 金幣+2, 謝謝交流~~ 2012-03-26 10:32:42

CODE:

clc

format compact

for x1=0:1

    for x2=0:1

        for x3=0:1

            for x4=0:1

                if x1^2==x1 & ...

                   x2^2==x2 & ...

                   x3^2==x3 & ...

                   x4^2==x4 & ...

                   x1^2==x1*x2*x3*x4

                disp([x1,x2,x3,x4]);

                end

            end

        end

    end

end

可得結(jié)果
   0    0    0    0
   0    0    0    1
   0    0    1    0
   0    0    1    1
   0    1    0    0
   0    1    0    1
   0    1    1    0
   0    1    1    1
   1    1    1    1

贊一下

回復此樓

3樓2012-03-25 21:32:30

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

snowwhite1

木蟲 (小有名氣)

應(yīng)助: 5 (幼兒園)
金幣: 4484.2
帖子: 275
在線: 104.1小時
蟲號: 1182357
注冊: 2011-01-03
性別: GG
專業(yè): 數(shù)論

引用回帖:

2樓: Originally posted by sogooda at 2012-03-25 21:26:47:
一共就4個數(shù)，每個數(shù)還都只能取0和1兩種情況，那就窮舉法就行了。
換句話說，編個4重循環(huán)，一下子就出來了。

謝謝，我只是舉例，實際情況是方程以及變量都很多，沒辦法窮舉的。但方程的解就只有0或1兩種情況，有沒有什么較好的辦法？

贊一下

回復此樓

奮斗

4樓2012-03-26 00:14:08

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

sogooda

金蟲 (正式寫手)

應(yīng)助: 19 (小學生)
金幣: 1497.4
紅花: 2
帖子: 737
在線: 94.6小時
蟲號: 438143
注冊: 2007-09-16
專業(yè): 防災(zāi)工程

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

★
臭水溝: 金幣+1, 謝謝交流~~ 2012-03-26 10:32:49

引用回帖:

4樓: Originally posted by snowwhite1 at 2012-03-26 00:14:08:
謝謝，我只是舉例，實際情況是方程以及變量都很多，沒辦法窮舉的。但方程的解就只有0或1兩種情況，有沒有什么較好的辦法？

暫時還沒有別的想法。
但是我還是固執(zhí)的認為窮舉法對于你的問題來說可能確實是一個很好的方法。
如果問題的規(guī)模很大，可以考慮一部分一部分的解決。

贊一下

回復此樓

5樓2012-03-26 08:53:55

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

snowwhite1

木蟲 (小有名氣)

應(yīng)助: 5 (幼兒園)
金幣: 4484.2
帖子: 275
在線: 104.1小時
蟲號: 1182357
注冊: 2011-01-03
性別: GG
專業(yè): 數(shù)論

沒有其他辦法了嗎？

贊一下

回復此樓

奮斗

6樓2012-03-27 23:13:22

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

dingd

鐵桿木蟲 (職業(yè)作家)

應(yīng)助: 1641 (講師)
金幣: 15037.3
散金: 101
紅花: 234
帖子: 3410
在線: 1223.7小時
蟲號: 291104
注冊: 2006-10-28

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

★ ★
xiegangmai: 金幣+2, 謝謝參與！ 2013-04-16 23:24:30

1stOpt求解：

CODE:

MultiRun = 30;

BinParameter x(4);

Function

x1^2=x1;

x2^2=x2;

x3^2=x3;

x4^2=x4;

x1^2=x1*x2*x3*x4;

結(jié)果：
No. x1 x2 x3 x4
1 0 1 1 1
2 0 0 0 1
3 0 1 0 0
4 0 0 1 1
5 0 0 1 0
6 0 1 0 1
7 1 1 1 1
8 0 0 0 0
9 0 1 1 0

贊一下

回復此樓

7樓2013-04-16 09:24:06

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

feign_te

金蟲 (小有名氣)

應(yīng)助: 15 (小學生)
金幣: 746
散金: 10
紅花: 2
帖子: 67
在線: 70小時
蟲號: 2550259
注冊: 2013-07-17
性別: GG
專業(yè): 質(zhì)譜分析

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

大概有多少個變量？
比如有n個變量，每個變量時0或者1，一個變量組就相當于一個2進制數(shù)字。
共有2^n組變量，要檢驗2^n次約束方程，1000萬次時n=23，得考慮時間問題了。

n很大的話可能要分析約束方程，提取公因子，整體考慮
比如這題約束方程為x1*（1-x2*x3*x4)=0

如果約束方程很復雜的話，從約束方程的系數(shù)入手，說不定有一般解法。
約束方程只有一個時，有2^n項，借助線性代數(shù)，可能會發(fā)現(xiàn)什么。

不知道7樓解法的極限，如果能很隨意的求出四五十個變量，二三十個約束方程下的解，那么幾乎可以斷定有一般解了。

贊一下

回復此樓

8樓2013-12-01 22:18:41

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

feign_te

金蟲 (小有名氣)

應(yīng)助: 15 (小學生)
金幣: 746
散金: 10
紅花: 2
帖子: 67
在線: 70小時
蟲號: 2550259
注冊: 2013-07-17
性別: GG
專業(yè): 質(zhì)譜分析

不過還是感覺求全部解不現(xiàn)實，解的數(shù)量接近2^(變量數(shù)-約束方程數(shù))

贊一下

回復此樓

9樓2013-12-01 22:25:16

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

feign_te

金蟲 (小有名氣)

應(yīng)助: 15 (小學生)
金幣: 746
散金: 10
紅花: 2
帖子: 67
在線: 70小時
蟲號: 2550259
注冊: 2013-07-17
性別: GG
專業(yè): 質(zhì)譜分析

引用回帖:

8樓: Originally posted by feign_te at 2013-12-01 22:18:41
大概有多少個變量？
比如有n個變量，每個變量時0或者1，一個變量組就相當于一個2進制數(shù)字。
共有2^n組變量，要檢驗2^n次約束方程，1000萬次時n=23，得考慮時間問題了。

n很大的話可能要分析約束方程，提取公因 ...

約束方程可以有無窮多個項，，

贊一下

回復此樓

10樓2013-12-02 07:54:42

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

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

[考研] 268求調(diào)劑 +3	簡單點0 2026-03-17	3/150	2026-03-17 23:15 by wangjihu
[考研] 311求調(diào)劑 +9	冬十三 2026-03-15	10/500	2026-03-17 22:53 by lbsjt
[考研] 304求調(diào)劑 +8	小熊joy 2026-03-14	8/400	2026-03-17 17:29 by ruiyingmiao
[考研] 274求調(diào)劑0856材料化工 +13	z2839474511 2026-03-11	14/700	2026-03-17 16:51 by share_joy
[考研] 求調(diào)劑，總分315，考的生物醫(yī)藥，一志愿湖南師范大學。調(diào)劑到任何專業(yè)都可以 +4	小丁想進步 2026-03-11	5/250	2026-03-17 16:05 by 外星文明
[考研] 08工科 320總分求調(diào)劑 +4	梨花珞晚風 2026-03-17	4/200	2026-03-17 13:38 by houyaoxu
[考研] 289求調(diào)劑 +6	步川酷紫123 2026-03-11	6/300	2026-03-17 10:23 by Sammy2
[考研] 283求調(diào)劑 +3	聽風就是雨； 2026-03-16	3/150	2026-03-17 07:41 by 熱情沙漠
[考研] 11408 一志愿西電，277分求調(diào)劑 +3	zhouzhen654 2026-03-16	3/150	2026-03-17 07:03 by laoshidan
[考研] 一志愿985，本科211，0817化學工程與技術(shù)319求調(diào)劑 +5	Liwangman 2026-03-15	5/250	2026-03-16 17:10 by 我的船我的海
[考研] 0703化學調(diào)劑 290分有科研經(jīng)歷，論文在投 +7	膩膩gk 2026-03-14	7/350	2026-03-16 10:12 by houyaoxu
[考研] 26考研一志愿中國石油大學(華東)305分求調(diào)劑 +3	嘉年新程 2026-03-15	3/150	2026-03-15 13:58 by 哈哈哈哈嘿嘿嘿
[考研] 330求調(diào)劑 +3	?醬給調(diào)劑跪了 2026-03-13	3/150	2026-03-14 10:13 by JourneyLucky
[考研] 招收0805（材料）調(diào)劑 +3	18595523086 2026-03-13	3/150	2026-03-14 00:33 by 123%、
[考研] 四川大學085601材料工程專碩初試294求調(diào)劑 +4	祝我們好在冬天 2026-03-11	4/200	2026-03-13 21:39 by peike
[考研] 26調(diào)劑/材料科學與工程/總分295/求收留 +9	2026調(diào)劑俠 2026-03-12	9/450	2026-03-13 20:46 by 18595523086
[考研] 求調(diào)劑 +5	一定有學上- 2026-03-12	5/250	2026-03-13 18:31 by ms629
[考研] 一志愿211化學學碩310分求調(diào)劑 +8	努力奮斗112 2026-03-12	9/450	2026-03-13 15:41 by JourneyLucky
[考研] 材料301分求調(diào)劑 +5	Liyouyumairs 2026-03-12	5/250	2026-03-13 14:42 by JourneyLucky
[考研] 化工學碩306求調(diào)劑 +9	42838695 2026-03-12	9/450	2026-03-13 10:16 by houyaoxu