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

mudan

1. for x1=1:13+ ^7 ~! c) @/ x- J
   
2.     for x2=1:137 t3 [4 S( @' \6 y8 A! P% D
   
3.         x3=(30000-4000*x1-3000*x2)/2000;
   6 h+ T  `3 w; L& y0 C% V  f# w
4.         m=13-x3-x1-x2;
   9 `1 J' A3 u+ ?+ g1 n) H1 o# y* |
5.         if sign(x3)>=0 &  rem(x3,1)==0 & sign(m)==0  \" c, A' \4 ~8 x9 {* x7 t+ |
   
6.             disp([x1,x2,x3])
   ; `$ D3 Z9 o0 e8 N4 o! ?& K
7.         end
   4 T: M9 j+ x8 ]3 A! K$ v8 f' ^  W9 T/ M
8.     end
   . i, i+ i& o' {$ m
9. end

复制代码

madio

rem这个函数可以换为mod，容易知道一些！
matlab中的rem()和mod()函数matlab的rem()和mod()函数
rem(x,y):求整除x/y的余数

mod(x,y):求模
( q5 d! j$ ~* `  r" U+ D* d- R
' a- W6 Q/ s1 ]4 h

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

& ?1 y9 b( k$ [! Z* Vmod(x,y)=x-y.*floor(x./y); （floor()向左取整）

6 B4 r1 x. n- K' p% s+ H' x& m; D如果x和y的符号相同（同为‘+’，同为‘-’），那么rem(x,y)=mod(x,y)（正数与正数，负数与负数，取整结果两个函数效果一样）
: N! _7 `$ E) f$ s# ~1 c+ H
9 Y2 e4 j/ p0 B: h; ]( J如果x和y的符号相反，那么mod(x,y)=rem(x,y)+y(正数与负数的取整，看你希望得到什么结果就选择用什么函数)
2 V# ]8 \' q, H& L* H, e/ P+ r
$ Y" U4 b' }4 b$ A  A0 v6 c1 e当正数与负数取余时，当得到的余数结果的符号希望跟除数(x)一样，用rem()函数；当得到的余数结果的符号希望跟被除数(y)一样，用mod()函数
  U4 d) @1 j# q8 F
例如：
) N% L$ }/ }0 E5 J+ @  m* w
0 e3 t4 a0 @+ T# S0 D& z9 frem(3,2)=1;mod(3,2)=1;

5 Y% _; O+ U2 |9 K7 B  erem(-3,-2)=-1;mod(-3,-2)=-1;
! z# p. g  v) J. ]; J8 U
rem(3,-2)=1;mod(3,-2)=-1;

# k" A# _7 ]& orem(-3,2)=-1;mod(-3,2)=1;
/ v6 M6 |' z( N0 ^
' y! u" j% G7 ]4 ~9 F/ `: Y" n. Z

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

: k" w+ s1 \, ?8 w" Emod(3,-2)=rem(3,-2)-2=1-2=-1;
! Q# v% y  \4 v( c% R; X
7 r- X5 w+ ^' a0 X- |0 U. J) p5 Dmod(-3,2)=rem(-3,2)+2=-1+2=1;

mudan

madio 发表于 2012-10-30 06:38
8 e' d5 y( P! U8 {5 O% trem这个函数可以换为mod，容易知道一些！
matlab中的rem()和mod()函数matlab的rem()和mod()函数
rem(x,y) ...

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

wy617958197

多谢分享资料！