求解4000*x1+3000*x2+2000*x3=30000的自然数解程序

mudan

1. for x1=1:13
   $ G# u0 T( a7 {& P
2.     for x2=1:131 Q# Y2 r, `' O9 ~\" H) _5 K; \6 w
   
3.         x3=(30000-4000*x1-3000*x2)/2000;$ G# A4 L0 M\" H5 b
   
4.         m=13-x3-x1-x2;
   9 ~5 i: K; E- ^! w3 k5 J$ r
5.         if sign(x3)>=0 &  rem(x3,1)==0 & sign(m)==0  ) o% x6 `* B& g5 \+ W8 G
   
6.             disp([x1,x2,x3])
   \" W& x5 z6 b+ B. V. C$ ^7 c
7.         end
   7 L1 W$ N! S* b9 a  R( d
8.     end* O- l9 m* G6 \2 h- L- Q
   
9. end

复制代码

madio

rem这个函数可以换为mod，容易知道一些！
1 O: k6 H! R, jmatlab中的rem()和mod()函数matlab的rem()和mod()函数
6 u6 G* ?5 l. H# H+ `: M5 G# J; {" ~rem(x,y):求整除x/y的余数

! x% T+ S, E1 i# w0 U* P  J! V0 Ymod(x,y):求模
& I* U. `( |/ z/ y0 {
% r2 `4 b, J5 l3 E  U

rem(x,y)=x-y.*fix(x./y);  (fix()向0取整)

mod(x,y)=x-y.*floor(x./y); （floor()向左取整）

; h( i0 E* u7 S! x; O+ d- \如果x和y的符号相同（同为‘+’，同为‘-’），那么rem(x,y)=mod(x,y)（正数与正数，负数与负数，取整结果两个函数效果一样）

如果x和y的符号相反，那么mod(x,y)=rem(x,y)+y(正数与负数的取整，看你希望得到什么结果就选择用什么函数)

8 w2 u1 d, k0 K) r当正数与负数取余时，当得到的余数结果的符号希望跟除数(x)一样，用rem()函数；当得到的余数结果的符号希望跟被除数(y)一样，用mod()函数

例如：
7 |1 |5 e' j0 k9 `% \* m
9 `0 v5 b6 {( P" @9 erem(3,2)=1;mod(3,2)=1;
8 M1 {- s4 E7 i- h; G
% h  g! @* Y3 E2 Z1 {rem(-3,-2)=-1;mod(-3,-2)=-1;
5 o* }; H0 J! q8 D2 V; L3 T
rem(3,-2)=1;mod(3,-2)=-1;
* w1 X  E" i# \
rem(-3,2)=-1;mod(-3,2)=1;
: P" h0 C- d/ J5 z. D3 d
- L# c9 o2 W- \% j

mod(3,2)=rem(3,2)=1;

mod(3,-2)=rem(3,-2)-2=1-2=-1;
" o0 ~- D) v, m* u; t
8 ?( K# M4 C& X0 h6 z$ E* C' `mod(-3,2)=rem(-3,2)+2=-1+2=1;

mudan

madio 发表于 2012-10-30 06:38
. `0 V, T# E. j! b5 lrem这个函数可以换为mod，容易知道一些！
8 t/ ?* ]7 @' H  E& |matlab中的rem()和mod()函数matlab的rem()和mod()函数
+ [% F$ m7 Y1 X$ |2 Jrem(x,y) ...

讲的很细致，坛主能不能推荐一本matlab程序实现的书？侧重实例的，谢谢

wy617958197

多谢分享资料！