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

mudan

1. for x1=1:138 e4 @2 @( B8 Z2 n% b
   
2.     for x2=1:13
   1 p# o  `: T; ~+ m8 r8 N
3.         x3=(30000-4000*x1-3000*x2)/2000;
   2 S& P/ {' P+ V. y: I6 Y5 J
4.         m=13-x3-x1-x2;
   7 p8 M! i, t6 Y4 x2 }\" @) D& m% I
5.         if sign(x3)>=0 &  rem(x3,1)==0 & sign(m)==0  + S* E\" b/ O4 q% l6 z9 s+ `9 o
   
6.             disp([x1,x2,x3])2 M8 S3 g, m' ^  o
   
7.         end\" L) G9 \7 t% J! N; m: @
   
8.     end
     H! |% j8 U9 @+ X  e
9. end

复制代码

madio

rem这个函数可以换为mod，容易知道一些！
matlab中的rem()和mod()函数matlab的rem()和mod()函数
rem(x,y):求整除x/y的余数
- o5 J& {( E$ @* t7 m/ L
  Z. z2 b( w( `& Bmod(x,y):求模
* n6 {- R( U8 s# ]4 d
! m; Z( g& |: V1 W

# J! @4 c, }: a2 B6 `3 x+ crem(x,y)=x-y.*fix(x./y);  (fix()向0取整)
% q$ m8 g4 C  H5 r, E
4 |1 ]9 o9 G$ c3 e) bmod(x,y)=x-y.*floor(x./y); （floor()向左取整）
2 b& ?5 o, ?' }. t! f
如果x和y的符号相同（同为‘+’，同为‘-’），那么rem(x,y)=mod(x,y)（正数与正数，负数与负数，取整结果两个函数效果一样）
$ B0 ?/ s( k2 t2 p3 |: h! S
如果x和y的符号相反，那么mod(x,y)=rem(x,y)+y(正数与负数的取整，看你希望得到什么结果就选择用什么函数)
) l3 h- j, |, k$ d
当正数与负数取余时，当得到的余数结果的符号希望跟除数(x)一样，用rem()函数；当得到的余数结果的符号希望跟被除数(y)一样，用mod()函数
/ Z4 ~0 p  h4 i4 q) O1 n+ N
例如：

3 Y' s1 _- k; ]rem(3,2)=1;mod(3,2)=1;

rem(-3,-2)=-1;mod(-3,-2)=-1;
" R  \/ [4 {8 b) U) C
$ @. x# ]9 R$ S' x- rrem(3,-2)=1;mod(3,-2)=-1;
- V3 @- Q9 T, P& B: s$ b3 n8 G8 O: h
# M7 T% ~  `- O# k0 A& ~+ {) Prem(-3,2)=-1;mod(-3,2)=1;

+ t: I. A+ \" T. ]2 B
. m7 y# f+ _) v1 h7 F) C1 L
mod(3,2)=rem(3,2)=1;
' b6 Z1 k' F. E
mod(3,-2)=rem(3,-2)-2=1-2=-1;

% k$ p* i" X" d& ?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()函数
7 N% j7 u0 [9 G! I8 Prem(x,y) ...

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

wy617958197

多谢分享资料！