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

mudan

1. for x1=1:13
   2 |0 T/ Z' T: T2 F# f! [: b
2.     for x2=1:13
   ; a' }7 J/ P4 f7 c0 A6 h
3.         x3=(30000-4000*x1-3000*x2)/2000;7 |5 q# |7 p* G$ k. e
   
4.         m=13-x3-x1-x2;
   $ g\" o( x6 b0 T9 _9 s
5.         if sign(x3)>=0 &  rem(x3,1)==0 & sign(m)==0  ) @3 f. f# C9 _/ O
   
6.             disp([x1,x2,x3])
   \" n3 O/ x  p% m' N* v
7.         end2 A: \4 X) h. q! [1 q4 H! R
   
8.     end3 a! [- A! \3 O
   
9. end

复制代码

madio

rem这个函数可以换为mod，容易知道一些！
matlab中的rem()和mod()函数matlab的rem()和mod()函数
; p0 o# y8 @$ q4 f  F2 srem(x,y):求整除x/y的余数

) V2 ?& H  D% l0 T; ymod(x,y):求模

! Y; h6 c0 U; ?, f9 e2 S3 b

rem(x,y)=x-y.*fix(x./y);  (fix()向0取整)

2 Z2 L- x* p. w0 N5 h: l+ ymod(x,y)=x-y.*floor(x./y); （floor()向左取整）
2 ]- C/ W! y9 ]- H
如果x和y的符号相同（同为‘+’，同为‘-’），那么rem(x,y)=mod(x,y)（正数与正数，负数与负数，取整结果两个函数效果一样）

如果x和y的符号相反，那么mod(x,y)=rem(x,y)+y(正数与负数的取整，看你希望得到什么结果就选择用什么函数)
" D2 K- _8 |1 X- c$ d
当正数与负数取余时，当得到的余数结果的符号希望跟除数(x)一样，用rem()函数；当得到的余数结果的符号希望跟被除数(y)一样，用mod()函数
+ G1 |3 f$ ~' m/ r/ @/ \3 }
例如：

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

rem(-3,-2)=-1;mod(-3,-2)=-1;
) D, v1 u8 m9 n
1 L) d0 S, W. G: o! Lrem(3,-2)=1;mod(3,-2)=-1;
: }/ B2 W* k' m4 [$ ~0 Y7 y# T
rem(-3,2)=-1;mod(-3,2)=1;
  p* r" g' Y" m! ~. ~1 k$ ~
! A8 p7 w5 }1 a* [8 s0 c

/ q- v* P: N. R! jmod(3,2)=rem(3,2)=1;

6 q- F3 n( _% E) w) |; _& bmod(3,-2)=rem(3,-2)-2=1-2=-1;

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

mudan

madio 发表于 2012-10-30 06:38
  w/ q0 F7 l3 v9 Hrem这个函数可以换为mod，容易知道一些！
% g* D% \4 @; h2 Smatlab中的rem()和mod()函数matlab的rem()和mod()函数
  d7 s# Z, R5 w$ Y$ v4 H: `# rrem(x,y) ...

+ V* F" N" V# c9 `讲的很细致，坛主能不能推荐一本matlab程序实现的书？侧重实例的，谢谢

wy617958197

多谢分享资料！