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

mudan

1. for x1=1:13( W9 ^' |& `& f& l, K% C% j
   
2.     for x2=1:13
   : E9 O0 ]' X( ^% x. s; a
3.         x3=(30000-4000*x1-3000*x2)/2000;
   2 i# M) C5 O* z$ t# \6 \
4.         m=13-x3-x1-x2;7 c/ U; v0 }0 \' A. ^  k7 w8 n
   
5.         if sign(x3)>=0 &  rem(x3,1)==0 & sign(m)==0  
   7 G$ v9 n  y3 I$ b5 K% N6 l, {/ X! O
6.             disp([x1,x2,x3])6 o! h5 F3 Q3 L2 b
   
7.         end# P6 L6 ^; Q- R: [5 M
   
8.     end2 [; j7 {3 _& a3 v
   
9. end

复制代码

madio

rem这个函数可以换为mod，容易知道一些！
. K/ p  S4 \3 h3 ~# Y( Q  J* Qmatlab中的rem()和mod()函数matlab的rem()和mod()函数
rem(x,y):求整除x/y的余数

2 J0 Q; h1 o7 V) Emod(x,y):求模

! s9 v5 ^) N; h" W2 Y7 |' N
/ s# v& J9 i: ?6 ^2 E! Q, e
rem(x,y)=x-y.*fix(x./y);  (fix()向0取整)
3 K( f% @1 f; W: n
2 P: u, o5 H. k6 J' c6 Qmod(x,y)=x-y.*floor(x./y); （floor()向左取整）

如果x和y的符号相同（同为‘+’，同为‘-’），那么rem(x,y)=mod(x,y)（正数与正数，负数与负数，取整结果两个函数效果一样）
) O* p. `, w0 [- M# v
: I* f( M) p" n* D; q如果x和y的符号相反，那么mod(x,y)=rem(x,y)+y(正数与负数的取整，看你希望得到什么结果就选择用什么函数)
% \0 F' I- D9 D0 \1 w
. r) k. T2 N- `7 X当正数与负数取余时，当得到的余数结果的符号希望跟除数(x)一样，用rem()函数；当得到的余数结果的符号希望跟被除数(y)一样，用mod()函数

例如：

rem(3,2)=1;mod(3,2)=1;
- ^) Q2 F9 Z* g: l3 |
0 Q2 S3 q$ T9 Z6 N  q+ grem(-3,-2)=-1;mod(-3,-2)=-1;
5 U" k+ W$ L! b; i6 n! n4 m
/ R+ q& q$ a% k( m6 m$ x- arem(3,-2)=1;mod(3,-2)=-1;
/ m) V. M! W- q- Y# x8 l
rem(-3,2)=-1;mod(-3,2)=1;
7 v0 N) Z+ s7 ]; `% v- C. Z; P
% _& q) _  x) S! M& j9 G6 g
# j6 n+ q8 P0 L6 d' G- Z# y5 c
mod(3,2)=rem(3,2)=1;
6 K' P4 t& P0 a' ?* d- o% R& v0 \
; e" w1 @, P9 v: B; @mod(3,-2)=rem(3,-2)-2=1-2=-1;
5 e; ?8 b4 i# C* V
mod(-3,2)=rem(-3,2)+2=-1+2=1;

mudan

madio 发表于 2012-10-30 06:38
rem这个函数可以换为mod，容易知道一些！
matlab中的rem()和mod()函数matlab的rem()和mod()函数
rem(x,y) ...

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

wy617958197

2 k. b7 I( ^7 a7 @多谢分享资料！