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

mudan

1. for x1=1:138 t% b0 L0 ]4 V1 \' A
   
2.     for x2=1:13, M% o8 ?! _5 r0 y! E
   
3.         x3=(30000-4000*x1-3000*x2)/2000;
   % u7 N  e6 o; d( d3 x7 J
4.         m=13-x3-x1-x2;
   : \% J6 h! C9 M4 I* _3 C' B% o0 v
5.         if sign(x3)>=0 &  rem(x3,1)==0 & sign(m)==0  
   ( x( G1 n6 a' p7 u* G, P- E' I5 Y
6.             disp([x1,x2,x3])! O' s* S( k5 r8 W
   
7.         end) a1 N, u* a! K5 I: f
   
8.     end
   # Y& u. P! f, l5 @
9. end

复制代码

madio

rem这个函数可以换为mod，容易知道一些！
# u/ R# m! a& N$ x+ l# bmatlab中的rem()和mod()函数matlab的rem()和mod()函数
rem(x,y):求整除x/y的余数
( u8 K, E$ m) f6 [; ]
mod(x,y):求模
2 G! M6 U) W& D/ t  t
# r4 R3 k, g4 L* _' J
  P+ A+ Q9 Z* g  Q5 Q
- m7 i4 A  G8 c8 rrem(x,y)=x-y.*fix(x./y);  (fix()向0取整)

1 V5 R* A7 K$ M8 Pmod(x,y)=x-y.*floor(x./y); （floor()向左取整）

& C) p2 A  T8 r' N! N9 Y% ^如果x和y的符号相同（同为‘+’，同为‘-’），那么rem(x,y)=mod(x,y)（正数与正数，负数与负数，取整结果两个函数效果一样）
9 s0 @& F5 U4 @% U
如果x和y的符号相反，那么mod(x,y)=rem(x,y)+y(正数与负数的取整，看你希望得到什么结果就选择用什么函数)
9 X: [7 f4 r+ s0 B
当正数与负数取余时，当得到的余数结果的符号希望跟除数(x)一样，用rem()函数；当得到的余数结果的符号希望跟被除数(y)一样，用mod()函数

例如：
/ q3 O' z0 W2 C& @! T, R7 R
rem(3,2)=1;mod(3,2)=1;

/ C. {; u2 j' c* }, Vrem(-3,-2)=-1;mod(-3,-2)=-1;

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

7 C) C- ^, V" |8 @8 M0 ]4 yrem(-3,2)=-1;mod(-3,2)=1;

  l) v3 O8 P5 Y0 M2 A

mod(3,2)=rem(3,2)=1;
6 z% ~( F! ?3 m6 w' ~
' V0 n2 u" ^4 U0 t, F) umod(3,-2)=rem(3,-2)-2=1-2=-1;
8 S; T& w/ l9 H! h8 }0 r
4 j8 f# f% N  K& _5 w# ~2 Emod(-3,2)=rem(-3,2)+2=-1+2=1;

mudan

madio 发表于 2012-10-30 06:38
rem这个函数可以换为mod，容易知道一些！
matlab中的rem()和mod()函数matlab的rem()和mod()函数
rem(x,y) ...

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

wy617958197

; a( |, k2 O+ H9 N9 l多谢分享资料！