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

mudan

1. for x1=1:13
   / N  a* v4 d) `$ H
2.     for x2=1:13% t; C  M3 T0 F2 i9 R
   
3.         x3=(30000-4000*x1-3000*x2)/2000;) O; C  i/ G5 _7 I0 z2 g% ?* L
   
4.         m=13-x3-x1-x2;4 ~8 W% w5 m8 g% V4 g0 b# j
   
5.         if sign(x3)>=0 &  rem(x3,1)==0 & sign(m)==0  # r\" P# S+ c7 ~! j. ]- R' T' ~
   
6.             disp([x1,x2,x3])
   + C3 f  V+ q& s( ]8 ~
7.         end
   ( s% w& S# I2 D. B( L) O
8.     end7 H0 X& W2 C; b1 g6 V
   
9. end

复制代码

madio

rem这个函数可以换为mod，容易知道一些！
matlab中的rem()和mod()函数matlab的rem()和mod()函数
, x% [+ r; N6 o8 Z( B9 A3 Yrem(x,y):求整除x/y的余数
1 F! [2 }1 |* i0 ]
# L) K- B, d! B9 v, w# k) ]6 [1 nmod(x,y):求模
  r. t+ w5 I" j6 o% @& S& X; d% ?


; @! K) p# ^+ D) f1 }rem(x,y)=x-y.*fix(x./y);  (fix()向0取整)
) ]6 L" M* k/ c
) ~. l5 s2 i/ |' Rmod(x,y)=x-y.*floor(x./y); （floor()向左取整）
: m) s" c2 I. a8 F6 J# e  j
如果x和y的符号相同（同为‘+’，同为‘-’），那么rem(x,y)=mod(x,y)（正数与正数，负数与负数，取整结果两个函数效果一样）
1 Z/ g: P! }8 N8 x2 p5 ~
5 O9 ^8 G1 d/ Q/ I. H9 o, M如果x和y的符号相反，那么mod(x,y)=rem(x,y)+y(正数与负数的取整，看你希望得到什么结果就选择用什么函数)
- H7 n3 Y! B. [& H( M. q
当正数与负数取余时，当得到的余数结果的符号希望跟除数(x)一样，用rem()函数；当得到的余数结果的符号希望跟被除数(y)一样，用mod()函数

# f! {0 T- |* z  R0 P例如：
  S+ J  R8 U7 Z- x" d7 ?) M
rem(3,2)=1;mod(3,2)=1;

) D; ~+ y& D& ^$ \$ D/ Qrem(-3,-2)=-1;mod(-3,-2)=-1;
' B1 S. w& |& ~3 M! l
, v6 Y9 V+ q; Yrem(3,-2)=1;mod(3,-2)=-1;
( \. |4 S0 {% z6 A$ I. A
rem(-3,2)=-1;mod(-3,2)=1;


) N; j% X; f2 x5 G* w$ E
4 F5 k5 `, e3 S% L2 ~mod(3,2)=rem(3,2)=1;
8 d$ y+ G5 F" w: v, ^/ v8 {* N
; E$ |& Y' Q* y# T. h* pmod(3,-2)=rem(3,-2)-2=1-2=-1;
3 I  ^4 |1 e) u' h- l: X: Q* e0 f
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()函数
( i2 X- t$ w! V0 Y  O' s# \0 |8 u$ X' Krem(x,y) ...

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

wy617958197

% Q: e9 N9 v& N; J; O多谢分享资料！