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

mudan

1. for x1=1:13* h% W: k7 U% i
   
2.     for x2=1:135 t3 G# }9 z' L6 @3 B+ X
   
3.         x3=(30000-4000*x1-3000*x2)/2000;- W* ?' \) c9 W7 M+ {
   
4.         m=13-x3-x1-x2;
   / Y! L6 \: I( B\" Z( ~: g\" s' d
5.         if sign(x3)>=0 &  rem(x3,1)==0 & sign(m)==0  
   ! t7 l0 ?6 G) i; {
6.             disp([x1,x2,x3])
   ' h9 e5 n4 S! D/ l/ ~$ P
7.         end& k9 C( {6 H2 u4 H6 _
   
8.     end2 U. g; i8 G* G! m: b2 g
   
9. end

复制代码

madio

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

0 p( \4 W% B& e3 K3 Ymod(x,y):求模
3 c! K+ O, E  R

& Y- L3 o; S5 g9 `3 x! k
rem(x,y)=x-y.*fix(x./y);  (fix()向0取整)
3 v: p8 f0 T! H- @: y
mod(x,y)=x-y.*floor(x./y); （floor()向左取整）
1 [" @4 j" o( n( n
- n+ X" Z) T0 Y如果x和y的符号相同（同为‘+’，同为‘-’），那么rem(x,y)=mod(x,y)（正数与正数，负数与负数，取整结果两个函数效果一样）

. Q4 j- J. Q9 M/ y3 y如果x和y的符号相反，那么mod(x,y)=rem(x,y)+y(正数与负数的取整，看你希望得到什么结果就选择用什么函数)
& _& d1 u/ O7 u: T: K
: `& S" F) ^# M( s/ z" i% _当正数与负数取余时，当得到的余数结果的符号希望跟除数(x)一样，用rem()函数；当得到的余数结果的符号希望跟被除数(y)一样，用mod()函数

例如：
" k  t  b  R( A/ v
3 W$ [4 I- k+ H* P+ _7 |rem(3,2)=1;mod(3,2)=1;

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

7 |) k7 L4 l3 _  c0 J! Krem(3,-2)=1;mod(3,-2)=-1;
" Z. r1 E5 p  F' M& Q
9 @. x' f' J1 }& T' K" j9 Rrem(-3,2)=-1;mod(-3,2)=1;
  }) d& W4 `2 _3 \1 N( x


4 E; A5 m6 L2 J. p5 Vmod(3,2)=rem(3,2)=1;

mod(3,-2)=rem(3,-2)-2=1-2=-1;
6 @6 y' c! K2 S# e; F: @% c
& [- Q6 ?, x: K) V. m* \mod(-3,2)=rem(-3,2)+2=-1+2=1;
1 K! K6 C3 @* b) B- d  r

mudan

madio 发表于 2012-10-30 06:38
9 I7 k7 \  l1 R' y0 u4 v6 h+ crem这个函数可以换为mod，容易知道一些！
matlab中的rem()和mod()函数matlab的rem()和mod()函数
rem(x,y) ...

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

wy617958197

多谢分享资料！