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

mudan

1. for x1=1:13
   ) C) Q. o1 ^5 R* p9 u
2.     for x2=1:13
   0 j1 _% l, ~3 k1 u% p) t' Y0 m
3.         x3=(30000-4000*x1-3000*x2)/2000;0 a0 t. U* h# O4 w9 ~
   
4.         m=13-x3-x1-x2;& ^, |6 y! J. P7 J
   
5.         if sign(x3)>=0 &  rem(x3,1)==0 & sign(m)==0  
   ; W' Q* R. O' x  u8 c2 @
6.             disp([x1,x2,x3])
     V6 F; C( G. ~$ f* b( m! _
7.         end, K) _8 w' n4 m8 e& s- K- `3 I+ \
   
8.     end
   ; @9 H7 f$ x6 y# b# [! I\" d
9. end

复制代码

madio

rem这个函数可以换为mod，容易知道一些！
matlab中的rem()和mod()函数matlab的rem()和mod()函数
rem(x,y):求整除x/y的余数
% l1 e) l. V8 u: Q) q6 e# J
5 m1 ~& `) j/ v& P) Zmod(x,y):求模
+ `, P7 j1 C6 u1 H( h$ {3 M

# J/ F* a' `2 M$ {0 x/ f
rem(x,y)=x-y.*fix(x./y);  (fix()向0取整)

* x' f+ b' H; h: U. y& D& nmod(x,y)=x-y.*floor(x./y); （floor()向左取整）

如果x和y的符号相同（同为‘+’，同为‘-’），那么rem(x,y)=mod(x,y)（正数与正数，负数与负数，取整结果两个函数效果一样）
$ p) }$ A: v) _9 C0 n/ s
如果x和y的符号相反，那么mod(x,y)=rem(x,y)+y(正数与负数的取整，看你希望得到什么结果就选择用什么函数)
9 e  w* ~" v3 g9 x9 m
当正数与负数取余时，当得到的余数结果的符号希望跟除数(x)一样，用rem()函数；当得到的余数结果的符号希望跟被除数(y)一样，用mod()函数

0 {. I9 G: h# w4 `' e+ o例如：
: W/ I  Z5 j, f9 A
' J, I% b( j  Q$ }1 I8 }6 c, Krem(3,2)=1;mod(3,2)=1;
0 O, h8 h% J5 I& N' {
$ j1 h+ O, p' V3 g1 t: J' D0 }rem(-3,-2)=-1;mod(-3,-2)=-1;

/ P9 g3 g. A0 V* i4 }  F: nrem(3,-2)=1;mod(3,-2)=-1;
) a% M; e  Y1 ?2 j' @
- p, J+ I1 J  S+ v1 k  nrem(-3,2)=-1;mod(-3,2)=1;


3 F" f, y; L2 o& i
0 Z6 K$ L4 d5 omod(3,2)=rem(3,2)=1;
2 Q, I) z" H" l2 }1 d+ Z7 S
mod(3,-2)=rem(3,-2)-2=1-2=-1;

# p8 E8 @4 P3 y6 X# K! h" hmod(-3,2)=rem(-3,2)+2=-1+2=1;
  M. Z' Q7 o# B2 z2 P! m2 q2 l

mudan

madio 发表于 2012-10-30 06:38
rem这个函数可以换为mod，容易知道一些！
matlab中的rem()和mod()函数matlab的rem()和mod()函数
& w5 W  q% e0 p+ \6 i) ^rem(x,y) ...

/ B2 i# b3 t: h( D; M; o2 W讲的很细致，坛主能不能推荐一本matlab程序实现的书？侧重实例的，谢谢

wy617958197

多谢分享资料！