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

mudan

1. for x1=1:13\" @, Z) g5 C2 {) l% q1 n
   
2.     for x2=1:13
   \" G/ r9 M4 ~1 `4 i  w
3.         x3=(30000-4000*x1-3000*x2)/2000;# F% T0 \1 D3 m9 ^
   
4.         m=13-x3-x1-x2;; v( M: D& m\" R: c\" g( O
   
5.         if sign(x3)>=0 &  rem(x3,1)==0 & sign(m)==0  
   0 r' b1 \. W  e
6.             disp([x1,x2,x3])# X1 _5 ]7 C4 o5 {# {
   
7.         end' x' r1 w. u: q# j: D  @
   
8.     end. k* P/ |) F6 T' G; [- `- K; J
   
9. end

复制代码

madio

rem这个函数可以换为mod，容易知道一些！
matlab中的rem()和mod()函数matlab的rem()和mod()函数
: r; G4 y3 z  L5 d* k4 X% Vrem(x,y):求整除x/y的余数
# X- o: v0 T# s2 y8 S# C
mod(x,y):求模
! l5 N! m5 o. m! V- z3 E5 _3 \" m0 j
8 q" ~1 h  @0 ?: T. b
4 l; a" ^% H$ O8 s6 l
5 v  a) q: c% y! v+ Mrem(x,y)=x-y.*fix(x./y);  (fix()向0取整)
0 u2 _6 ]+ {" A
mod(x,y)=x-y.*floor(x./y); （floor()向左取整）

) \' ^' V& ?5 @$ Q5 b如果x和y的符号相同（同为‘+’，同为‘-’），那么rem(x,y)=mod(x,y)（正数与正数，负数与负数，取整结果两个函数效果一样）

如果x和y的符号相反，那么mod(x,y)=rem(x,y)+y(正数与负数的取整，看你希望得到什么结果就选择用什么函数)
- D$ |* d2 E; T
* t& M+ z: A  e  L7 \: z当正数与负数取余时，当得到的余数结果的符号希望跟除数(x)一样，用rem()函数；当得到的余数结果的符号希望跟被除数(y)一样，用mod()函数

例如：
$ @5 o4 M: k% V. x  f
rem(3,2)=1;mod(3,2)=1;
+ m& x7 F0 T$ [- w' p
' q8 h8 o) Q2 @. v% Orem(-3,-2)=-1;mod(-3,-2)=-1;
& A: Y4 x6 Z1 A; f) k; m
0 ^5 Q) Y/ U/ w. j' prem(3,-2)=1;mod(3,-2)=-1;
2 K8 y& C& v$ a8 H
rem(-3,2)=-1;mod(-3,2)=1;
, P) L4 W9 w, S2 V
. T3 d) X2 I. p! V) b
7 E. r* b8 R9 ~/ B0 T8 `) g6 U
- z$ V! d) p8 i0 |- q+ H2 d0 G, ?* ~mod(3,2)=rem(3,2)=1;
9 j: i; g- _6 ^4 x% n5 r, A7 g/ k
mod(3,-2)=rem(3,-2)-2=1-2=-1;

' F: Z2 Y$ |1 [5 @mod(-3,2)=rem(-3,2)+2=-1+2=1;
. M, [) B6 [3 x$ M  ]+ [1 O- X

mudan

madio 发表于 2012-10-30 06:38
$ o: P( b& p) grem这个函数可以换为mod，容易知道一些！
matlab中的rem()和mod()函数matlab的rem()和mod()函数
5 M( a5 A6 z; Qrem(x,y) ...

8 D! E7 S1 B* w0 h& o* M2 K讲的很细致，坛主能不能推荐一本matlab程序实现的书？侧重实例的，谢谢

wy617958197

多谢分享资料！