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

mudan

1. for x1=1:13* s- o! m' ]: U% Y, o$ {* l
   
2.     for x2=1:13' p( r! J! `' e; x) S
   
3.         x3=(30000-4000*x1-3000*x2)/2000;  o) c# L$ t6 p* D% ?- W
   
4.         m=13-x3-x1-x2;5 s7 y, a$ @  x) }0 p
   
5.         if sign(x3)>=0 &  rem(x3,1)==0 & sign(m)==0  
   2 S! G: \1 m& f1 G4 O
6.             disp([x1,x2,x3])
   # }8 i- z( p. S
7.         end! M# I& s! s0 E: V& r6 A, y
   
8.     end
   1 K: u\" D8 @0 n. z+ {1 v# }8 y
9. end

复制代码

wy617958197

多谢分享资料！

mudan

madio 发表于 2012-10-30 06:38
rem这个函数可以换为mod，容易知道一些！
) z- t1 H- u: ~1 e- y* x+ V- v* Dmatlab中的rem()和mod()函数matlab的rem()和mod()函数
rem(x,y) ...

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

madio

rem这个函数可以换为mod，容易知道一些！
matlab中的rem()和mod()函数matlab的rem()和mod()函数
# g9 D* P" K: }& Qrem(x,y):求整除x/y的余数
9 W4 j- V6 `' e# r9 b+ V
mod(x,y):求模

4 O1 p% j3 e. D# j
( ?$ c6 Z( ^/ d% ^0 Z
7 j$ z* `7 G' E/ frem(x,y)=x-y.*fix(x./y);  (fix()向0取整)
. X$ ~9 W/ t: X4 y
- B* Z  `/ ^* A, W5 ]& Cmod(x,y)=x-y.*floor(x./y); （floor()向左取整）

% b8 M3 A$ |/ \/ f如果x和y的符号相同（同为‘+’，同为‘-’），那么rem(x,y)=mod(x,y)（正数与正数，负数与负数，取整结果两个函数效果一样）

* E6 q. g8 `4 d2 d" |1 f7 s如果x和y的符号相反，那么mod(x,y)=rem(x,y)+y(正数与负数的取整，看你希望得到什么结果就选择用什么函数)

当正数与负数取余时，当得到的余数结果的符号希望跟除数(x)一样，用rem()函数；当得到的余数结果的符号希望跟被除数(y)一样，用mod()函数
9 [  e* n( Z+ Y
2 w0 N/ v+ c6 ^4 F5 a5 i例如：

rem(3,2)=1;mod(3,2)=1;
6 W$ W- t( `' T! `
7 l/ g" I2 X3 Crem(-3,-2)=-1;mod(-3,-2)=-1;
/ _4 l) _% \. F" K8 f8 e
rem(3,-2)=1;mod(3,-2)=-1;
7 _% N1 {. u8 E( m) f  B
rem(-3,2)=-1;mod(-3,2)=1;
% ^9 S/ C: K9 s- s3 `5 g

5 x. a0 O  }9 w1 H5 L
" {: t* L+ H8 u- O3 Kmod(3,2)=rem(3,2)=1;

mod(3,-2)=rem(3,-2)-2=1-2=-1;
) o: Y. e- |& N' M& x7 S
1 @, c: V8 X3 v; a$ Q9 g. w% `  t+ _mod(-3,2)=rem(-3,2)+2=-1+2=1;