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

mudan

1. for x1=1:13
   2 N0 P+ A5 M6 I1 ?4 y
2.     for x2=1:13& ?+ J: D: E$ Y0 b2 [
   
3.         x3=(30000-4000*x1-3000*x2)/2000;
     c' d/ G5 R3 E
4.         m=13-x3-x1-x2;
   $ W9 \! j! l* G\" z
5.         if sign(x3)>=0 &  rem(x3,1)==0 & sign(m)==0  8 h) z% d+ G' N' p3 S
   
6.             disp([x1,x2,x3])) U# ^7 e# }5 e2 X6 F
   
7.         end4 E\" P' `6 w: i7 T: U
   
8.     end
   3 z  O\" e: H7 L; ~  h: C* w
9. end

复制代码

madio

rem这个函数可以换为mod，容易知道一些！
0 M, v; N1 N# a+ Y" omatlab中的rem()和mod()函数matlab的rem()和mod()函数
rem(x,y):求整除x/y的余数

mod(x,y):求模
' M) V6 T9 y$ I) F- m# G, y* n
) F5 c4 _3 f! f0 p. y- l: A

* Y/ D  F1 |1 e$ p7 d1 s9 rrem(x,y)=x-y.*fix(x./y);  (fix()向0取整)
7 i1 n5 i5 P/ G
mod(x,y)=x-y.*floor(x./y); （floor()向左取整）

如果x和y的符号相同（同为‘+’，同为‘-’），那么rem(x,y)=mod(x,y)（正数与正数，负数与负数，取整结果两个函数效果一样）
+ }* r- y, `; V9 C1 ?; i0 V9 h; ^& A
" M7 t+ q) Q$ c; [) A3 f如果x和y的符号相反，那么mod(x,y)=rem(x,y)+y(正数与负数的取整，看你希望得到什么结果就选择用什么函数)
3 z. R, m# U( c
8 T5 k) M, @; w- f: R4 h当正数与负数取余时，当得到的余数结果的符号希望跟除数(x)一样，用rem()函数；当得到的余数结果的符号希望跟被除数(y)一样，用mod()函数
3 c3 {' |2 I1 r- [* q
例如：

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

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

& D8 \" p8 K4 }rem(3,-2)=1;mod(3,-2)=-1;
& n& u! [8 X5 z" B2 u1 E
8 u; L" Y. l* T6 h: Y/ Y, _& |rem(-3,2)=-1;mod(-3,2)=1;
3 s% q1 w7 E. A1 x4 W" m. p) `
, |$ L. i' r& o( H4 r0 L  H

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

mod(3,-2)=rem(3,-2)-2=1-2=-1;
4 z" ]+ d; {1 C* U
/ P+ u5 F: M; H# W! Z0 \. k* E3 qmod(-3,2)=rem(-3,2)+2=-1+2=1;

mudan

madio 发表于 2012-10-30 06:38
8 S2 f, I/ l7 F; s9 \% Yrem这个函数可以换为mod，容易知道一些！
matlab中的rem()和mod()函数matlab的rem()和mod()函数
: I( p1 s# E& ]4 h+ m/ j0 Drem(x,y) ...

* l/ v& E5 l; f讲的很细致，坛主能不能推荐一本matlab程序实现的书？侧重实例的，谢谢

wy617958197

9 I" C& I3 }  ~* S多谢分享资料！