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）。

回復(fù)此樓

» 猜你喜歡

323求調(diào)劑已經(jīng)有6人回復(fù)
一志愿北京化工大學(xué) 070300 學(xué)碩 336分求調(diào)劑已經(jīng)有4人回復(fù)
352求調(diào)劑已經(jīng)有3人回復(fù)
一志愿東華大學(xué)化學(xué)070300，求調(diào)劑已經(jīng)有8人回復(fù)
277材料科學(xué)與工程080500求調(diào)劑已經(jīng)有7人回復(fù)
317求調(diào)劑已經(jīng)有18人回復(fù)
293求調(diào)劑已經(jīng)有5人回復(fù)
280分求調(diào)劑一志愿085802 已經(jīng)有7人回復(fù)
0854電子信息求調(diào)劑已經(jīng)有3人回復(fù)
263求調(diào)劑已經(jīng)有4人回復(fù)

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

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

奮斗

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

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

sogooda

金蟲 (正式寫手)

應(yīng)助: 19 (小學(xué)生)
金幣: 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)，一下子就出來了。

贊一下

回復(fù)此樓

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

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

sogooda

金蟲 (正式寫手)

應(yīng)助: 19 (小學(xué)生)
金幣: 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

贊一下

回復(fù)此樓

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

已閱回復(fù)此樓關(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)，一下子就出來了。

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

贊一下

回復(fù)此樓

奮斗

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

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

sogooda

金蟲 (正式寫手)

應(yīng)助: 19 (小學(xué)生)
金幣: 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:
謝謝，我只是舉例，實(shí)際情況是方程以及變量都很多，沒辦法窮舉的。但方程的解就只有0或1兩種情況，有沒有什么較好的辦法？

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

贊一下

回復(fù)此樓

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

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

snowwhite1

木蟲 (小有名氣)

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

沒有其他辦法了嗎？

贊一下

回復(fù)此樓

奮斗

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

已閱回復(fù)此樓關(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

贊一下

回復(fù)此樓

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

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

feign_te

金蟲 (小有名氣)

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

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

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

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

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

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

贊一下

回復(fù)此樓

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

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

feign_te

金蟲 (小有名氣)

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

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

贊一下

回復(fù)此樓

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

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

feign_te

金蟲 (小有名氣)

應(yīng)助: 15 (小學(xué)生)
金幣: 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，一個變量組就相當(dāng)于一個2進(jìn)制數(shù)字。
共有2^n組變量，要檢驗(yàn)2^n次約束方程，1000萬次時n=23，得考慮時間問題了。

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

約束方程可以有無窮多個項(xiàng)，，

贊一下

回復(fù)此樓

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

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

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

[考研] 一志愿武理材料工程348求調(diào)劑 +5	￣^￣゜汗 2026-03-19	7/350	2026-03-22 19:44 by 公瑾逍遙
[考研] 289求調(diào)劑 +7	懷瑾握瑜l 2026-03-20	7/350	2026-03-22 15:57 by ColorlessPI
[考研] 一志愿中南大學(xué)化學(xué)學(xué)碩0703總分337求調(diào)劑 +3	niko- 2026-03-22	3/150	2026-03-22 15:15 by 楊楊楊紫
[考研] 384求調(diào)劑 +3	子系博 2026-03-22	4/200	2026-03-22 11:04 by 搏擊518
[考研] 材料學(xué)碩333求調(diào)劑 +3	北道巷 2026-03-18	3/150	2026-03-21 18:17 by 學(xué)員8dgXkO
[考研] 313求調(diào)劑 +4	肆叁貳壹22 2026-03-19	4/200	2026-03-21 17:33 by ColorlessPI
[考研] 材料 271求調(diào)劑 +5	展信悅_ 2026-03-21	5/250	2026-03-21 17:29 by 學(xué)員8dgXkO
[考研] 296求調(diào)劑 +4	www_q 2026-03-20	4/200	2026-03-21 17:26 by 學(xué)員8dgXkO
[考研] 求調(diào)劑 +3	.m.. 2026-03-21	4/200	2026-03-21 16:25 by barlinike
[考研] 279求調(diào)劑 +5	紅衣隱官 2026-03-21	5/250	2026-03-21 14:59 by lature00
[考研] 華東師范大學(xué)-071000生物學(xué)-293分-求調(diào)劑 +3	研究生何瑤明 2026-03-18	3/150	2026-03-21 01:30 by JourneyLucky
[考研] 一志愿南昌大學(xué)，327分，材料與化工085600 +9	Ncdx123456 2026-03-19	9/450	2026-03-20 23:41 by lovewei0727
[考研] 中南大學(xué)化學(xué)學(xué)碩337求調(diào)劑 +3	niko- 2026-03-19	6/300	2026-03-20 21:58 by luoyongfeng
[考研] 材料與化工專碩調(diào)劑 +7	heming3743 2026-03-16	7/350	2026-03-20 19:31 by zhukairuo
[考研] 一志愿中國海洋大學(xué)，生物學(xué)，301分，求調(diào)劑 +5	1孫悟空 2026-03-17	6/300	2026-03-19 23:46 by zcl123
[考研] 085601材料工程專碩求調(diào)劑 +10	慕寒mio 2026-03-16	10/500	2026-03-19 15:26 by 丁丁*
[考研] 312求調(diào)劑 +8	陌宸希 2026-03-16	9/450	2026-03-18 12:39 by Linda Hu
[考研] 一志愿蘇州大學(xué)材料工程（085601）專碩有科研經(jīng)歷三項(xiàng)國獎兩個實(shí)用型專利一項(xiàng)省級立項(xiàng) +6	大火山小火山 2026-03-16	8/400	2026-03-17 15:05 by 無懈可擊111
[論文投稿] 有沒有大佬發(fā)小論文能帶我個二作 +3	增銳漏人 2026-03-17	4/200	2026-03-17 09:26 by xs74101122
[考研] 11408 一志愿西電，277分求調(diào)劑 +3	zhouzhen654 2026-03-16	3/150	2026-03-17 07:03 by laoshidan