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

mudan

1. for x1=1:13
   / f4 h- l5 b1 }+ u
2.     for x2=1:13( p# J' ]! O  n4 D/ m& C* ~& L5 w
   
3.         x3=(30000-4000*x1-3000*x2)/2000;  ?9 i\" x& e+ Q) {4 E
   
4.         m=13-x3-x1-x2;
   . _1 Q+ Y& [9 u6 J/ \
5.         if sign(x3)>=0 &  rem(x3,1)==0 & sign(m)==0  1 y) [4 K/ q1 e; x0 U
   
6.             disp([x1,x2,x3])4 Y( K5 _* n. u. }- m6 b
   
7.         end( X3 Q3 t7 S, n9 `9 c' o; E
   
8.     end
   1 [, N3 W7 Y  O* N
9. end

复制代码

madio

rem这个函数可以换为mod，容易知道一些！
* d& {, j/ h# d- tmatlab中的rem()和mod()函数matlab的rem()和mod()函数
rem(x,y):求整除x/y的余数
2 y4 D- [# A4 \5 N' f+ m  P
. I( e; @' F/ L. y4 C" Mmod(x,y):求模
% Q- t+ h2 X' h' V8 X


  s4 f# ~/ v. M; Y/ jrem(x,y)=x-y.*fix(x./y);  (fix()向0取整)

+ Y& ~* R( o& Gmod(x,y)=x-y.*floor(x./y); （floor()向左取整）
! }9 s! x- t' ]* z' _
如果x和y的符号相同（同为‘+’，同为‘-’），那么rem(x,y)=mod(x,y)（正数与正数，负数与负数，取整结果两个函数效果一样）
  m8 ^; @+ y& n0 X, n. N
如果x和y的符号相反，那么mod(x,y)=rem(x,y)+y(正数与负数的取整，看你希望得到什么结果就选择用什么函数)

$ B5 B1 t2 m4 _' ]/ `9 m9 u当正数与负数取余时，当得到的余数结果的符号希望跟除数(x)一样，用rem()函数；当得到的余数结果的符号希望跟被除数(y)一样，用mod()函数

. X# m. K4 q2 R  N& C例如：
5 E% K% M, S; K5 Z1 p
rem(3,2)=1;mod(3,2)=1;

; L  r  q7 I# w5 I0 G; E1 Mrem(-3,-2)=-1;mod(-3,-2)=-1;

rem(3,-2)=1;mod(3,-2)=-1;
+ r, C: O/ h. o  L. R+ v+ p3 V2 k, ]
rem(-3,2)=-1;mod(-3,2)=1;


! I0 @& t4 a* t7 v
mod(3,2)=rem(3,2)=1;

3 Y( U* b+ i* N  Hmod(3,-2)=rem(3,-2)-2=1-2=-1;
7 w8 P6 ]5 x! m6 X1 V8 F; {( M
6 U% ~+ q' r' {1 t- R* c% N+ nmod(-3,2)=rem(-3,2)+2=-1+2=1;
4 K" `9 \6 O0 g2 F6 M. \

mudan

madio 发表于 2012-10-30 06:38
$ `3 n( c) m9 e' P8 r# xrem这个函数可以换为mod，容易知道一些！
+ h) T; P8 R% J5 j9 @matlab中的rem()和mod()函数matlab的rem()和mod()函数
rem(x,y) ...

; v5 c1 Y* p- Z& A讲的很细致，坛主能不能推荐一本matlab程序实现的书？侧重实例的，谢谢

wy617958197

! F& W# n7 ~5 O多谢分享资料！