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

mudan

1. for x1=1:13
   2 t7 E7 f3 p9 o7 e
2.     for x2=1:13
   , l& T( V: Q+ H( G
3.         x3=(30000-4000*x1-3000*x2)/2000;
   : h3 P: L! u: E
4.         m=13-x3-x1-x2;. w. p0 g1 r, S* g  B7 T
   
5.         if sign(x3)>=0 &  rem(x3,1)==0 & sign(m)==0  
   4 A8 R0 O2 M. J% L0 N
6.             disp([x1,x2,x3])
     f& k4 T1 C* \6 y/ v- M: v
7.         end6 C8 ]; q/ q\" X; B
   
8.     end
     Q5 \/ _) B# Q\" a2 Q
9. end

复制代码

madio

rem这个函数可以换为mod，容易知道一些！
matlab中的rem()和mod()函数matlab的rem()和mod()函数
4 C7 @0 \2 w) V6 `2 k! X2 Nrem(x,y):求整除x/y的余数

mod(x,y):求模


' A  }. a( ^! V) C, r# ?
rem(x,y)=x-y.*fix(x./y);  (fix()向0取整)
! D$ q+ t0 V. [" {5 D0 D
mod(x,y)=x-y.*floor(x./y); （floor()向左取整）

' p4 {& R; n! e+ N/ i& @如果x和y的符号相同（同为‘+’，同为‘-’），那么rem(x,y)=mod(x,y)（正数与正数，负数与负数，取整结果两个函数效果一样）

2 E% K8 ^+ v. g2 z; `如果x和y的符号相反，那么mod(x,y)=rem(x,y)+y(正数与负数的取整，看你希望得到什么结果就选择用什么函数)
( _; O$ l6 _2 u( F, K4 B' |- q
3 P# W" H1 D  P! k6 j当正数与负数取余时，当得到的余数结果的符号希望跟除数(x)一样，用rem()函数；当得到的余数结果的符号希望跟被除数(y)一样，用mod()函数

例如：

$ h$ \/ B  I7 m* ~! x3 _) ?% Grem(3,2)=1;mod(3,2)=1;
) w' }; u) r, B4 x9 p9 {8 i
1 Y  M3 l$ W: ^/ e' `" C  erem(-3,-2)=-1;mod(-3,-2)=-1;
' f: k3 M" V% \* g0 q1 E
3 V' C1 f6 r8 hrem(3,-2)=1;mod(3,-2)=-1;

! ^0 v! Y; a6 j5 z) h" d3 [rem(-3,2)=-1;mod(-3,2)=1;


7 W6 ?" g8 k: A4 S8 X' j8 ^+ I
( Z6 _8 Z, y2 H( Dmod(3,2)=rem(3,2)=1;
1 F' d+ b- E" F8 I
mod(3,-2)=rem(3,-2)-2=1-2=-1;
* M% u1 `+ }" G  d1 p
( a  S4 l) n5 x- w! gmod(-3,2)=rem(-3,2)+2=-1+2=1;
% k1 s& D, ~# u4 r

mudan

madio 发表于 2012-10-30 06:38
& q# W( b9 t9 ]/ f, E' prem这个函数可以换为mod，容易知道一些！
2 U* {% x( d4 L( g- d/ a. K0 Jmatlab中的rem()和mod()函数matlab的rem()和mod()函数
rem(x,y) ...

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

wy617958197

多谢分享资料！