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

mudan

1. for x1=1:13
   / H) i$ {1 J5 Y; q, P
2.     for x2=1:13
   4 D  j\" e0 {! J
3.         x3=(30000-4000*x1-3000*x2)/2000;
   6 l# g8 N- f; O7 h' S( \
4.         m=13-x3-x1-x2;
   + H; Y* b& |1 w4 F% N5 s& y
5.         if sign(x3)>=0 &  rem(x3,1)==0 & sign(m)==0  9 ?5 _' t% k) _! n/ Y& O
   
6.             disp([x1,x2,x3]): _7 w& A+ A. \6 R
   
7.         end8 U! I* S2 T1 J
   
8.     end
   ' e, {* ~/ C7 b/ G# g3 y7 w
9. end

复制代码

madio

rem这个函数可以换为mod，容易知道一些！
+ |$ Y, M; Z$ `2 i2 e$ a* Smatlab中的rem()和mod()函数matlab的rem()和mod()函数
rem(x,y):求整除x/y的余数

( t4 E7 B1 L4 y+ x0 n; rmod(x,y):求模
2 ^* U$ L; |, n# q
, ]" a" i! u' x, H# U$ }

. n$ t) O) k5 j2 M! v& G4 Wrem(x,y)=x-y.*fix(x./y);  (fix()向0取整)

mod(x,y)=x-y.*floor(x./y); （floor()向左取整）
  Y' h, Z2 \5 P" A; Y. S- P: K
* b# q6 u8 G* O% \7 r如果x和y的符号相同（同为‘+’，同为‘-’），那么rem(x,y)=mod(x,y)（正数与正数，负数与负数，取整结果两个函数效果一样）
7 @5 l+ C" Q8 B
; R6 ~$ [# V2 P; `如果x和y的符号相反，那么mod(x,y)=rem(x,y)+y(正数与负数的取整，看你希望得到什么结果就选择用什么函数)
4 `% n1 D' @" R7 S
当正数与负数取余时，当得到的余数结果的符号希望跟除数(x)一样，用rem()函数；当得到的余数结果的符号希望跟被除数(y)一样，用mod()函数
9 i# H. X. y* H4 h" X
例如：
7 X1 N. u3 a/ n; x; \; Y
% W# g  W3 K7 Frem(3,2)=1;mod(3,2)=1;

9 \) h8 Z. K) A- Arem(-3,-2)=-1;mod(-3,-2)=-1;

; p% }" w4 _1 nrem(3,-2)=1;mod(3,-2)=-1;

rem(-3,2)=-1;mod(-3,2)=1;
& q1 z0 ~, x# ?9 v
" r! r- J( D" f; `  ~" T2 [

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

: |- W  J, v# u) Xmod(3,-2)=rem(3,-2)-2=1-2=-1;

$ d8 l7 t' m5 ~3 Kmod(-3,2)=rem(-3,2)+2=-1+2=1;

mudan

madio 发表于 2012-10-30 06:38
rem这个函数可以换为mod，容易知道一些！
0 P$ s, N2 y& f0 cmatlab中的rem()和mod()函数matlab的rem()和mod()函数
, P, D* C0 W/ j( ]' r/ {, _" zrem(x,y) ...

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

wy617958197

多谢分享资料！