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

mudan

1. for x1=1:13
   ! {# s) F, S: ^9 @5 y; I# V
2.     for x2=1:139 ^: {2 g4 l& G
   
3.         x3=(30000-4000*x1-3000*x2)/2000;3 [8 I' j2 K: X; x\" z4 Y
   
4.         m=13-x3-x1-x2;
   7 i7 }: c- v: y1 P, Y$ m
5.         if sign(x3)>=0 &  rem(x3,1)==0 & sign(m)==0  , b1 D' z- Q( g' Z
   
6.             disp([x1,x2,x3])
   4 b0 T. r4 R2 M8 r4 n( L$ E\" W
7.         end
   & g; h3 `+ H# ~0 \\" |/ l/ |0 G3 `1 A
8.     end
   4 l. F, Y4 `, n; l
9. end

复制代码

madio

rem这个函数可以换为mod，容易知道一些！
matlab中的rem()和mod()函数matlab的rem()和mod()函数
rem(x,y):求整除x/y的余数
3 T( `$ i+ R2 w
mod(x,y):求模
$ W* N4 e2 `6 f$ ^; V0 x# s; {3 b

1 h8 N. \  l: Z6 {9 Z* g
rem(x,y)=x-y.*fix(x./y);  (fix()向0取整)

- ^. U  n2 T  |1 S4 J( ?5 O* Z5 fmod(x,y)=x-y.*floor(x./y); （floor()向左取整）

如果x和y的符号相同（同为‘+’，同为‘-’），那么rem(x,y)=mod(x,y)（正数与正数，负数与负数，取整结果两个函数效果一样）
. G8 V7 w  @- r, H4 e0 _: x# r% Z
如果x和y的符号相反，那么mod(x,y)=rem(x,y)+y(正数与负数的取整，看你希望得到什么结果就选择用什么函数)

; O' Z8 _. F. \6 d) q- H3 X当正数与负数取余时，当得到的余数结果的符号希望跟除数(x)一样，用rem()函数；当得到的余数结果的符号希望跟被除数(y)一样，用mod()函数

例如：
; l; E# S4 E. w
. P$ T* y! x3 m# E. \2 b# jrem(3,2)=1;mod(3,2)=1;
" ]8 z+ N1 b% d! m. d) ?, C
, j1 @) f, I7 |6 x  zrem(-3,-2)=-1;mod(-3,-2)=-1;

, w6 ^( N" F" {, a4 ?- Grem(3,-2)=1;mod(3,-2)=-1;

+ b9 H8 z3 r6 t4 E+ X/ W" W4 c( Rrem(-3,2)=-1;mod(-3,2)=1;
, [4 R) W" c  e/ k
: V( `) K0 q! [/ @" m1 q2 K

7 G9 u# r, E& \mod(3,2)=rem(3,2)=1;

! ~) @9 j: O3 Rmod(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
rem这个函数可以换为mod，容易知道一些！
; Z* G, a+ _2 e1 X( h- E- w$ X* nmatlab中的rem()和mod()函数matlab的rem()和mod()函数
9 N0 [5 `, x- s2 B( Wrem(x,y) ...

1 l% R; S. F2 O! L/ v% U9 J讲的很细致，坛主能不能推荐一本matlab程序实现的书？侧重实例的，谢谢

wy617958197

$ R* B9 f% A. r7 ?0 Q多谢分享资料！