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

mudan

1. for x1=1:13
   3 I$ o) O8 }. R. |7 W
2.     for x2=1:13
   6 l& g9 g3 U* O' p6 \, H$ K8 r
3.         x3=(30000-4000*x1-3000*x2)/2000;$ i# K4 h* J! _+ u7 E  b
   
4.         m=13-x3-x1-x2;  B4 P% m) h' V6 }( s7 K
   
5.         if sign(x3)>=0 &  rem(x3,1)==0 & sign(m)==0  / L; l0 c5 M( B/ c* L, v- x
   
6.             disp([x1,x2,x3])
   : i# [. U5 Q) m
7.         end
   1 u% H1 j1 y- P7 y
8.     end
   ) `* d. \% b1 f! e  Y
9. end

复制代码

madio

rem这个函数可以换为mod，容易知道一些！
matlab中的rem()和mod()函数matlab的rem()和mod()函数
rem(x,y):求整除x/y的余数
0 \4 R: y+ z: \) }. a! b
mod(x,y):求模
" m' }0 y7 k) a, Y. v

1 X3 c% E& g' h% ~: r
  z$ ]/ j: U' e1 Qrem(x,y)=x-y.*fix(x./y);  (fix()向0取整)
( d5 i  U" A& [/ N6 G9 Q) r2 U
; A4 n1 k7 p& o$ ^- I3 Qmod(x,y)=x-y.*floor(x./y); （floor()向左取整）
& O3 C  C3 {7 ~0 h, v
: M& W: H+ g' l3 D; P如果x和y的符号相同（同为‘+’，同为‘-’），那么rem(x,y)=mod(x,y)（正数与正数，负数与负数，取整结果两个函数效果一样）

如果x和y的符号相反，那么mod(x,y)=rem(x,y)+y(正数与负数的取整，看你希望得到什么结果就选择用什么函数)

当正数与负数取余时，当得到的余数结果的符号希望跟除数(x)一样，用rem()函数；当得到的余数结果的符号希望跟被除数(y)一样，用mod()函数

! U! j2 W8 c$ T4 O3 ^: {例如：

rem(3,2)=1;mod(3,2)=1;
& H! Z; M, T- C3 i1 r4 F. z
8 u! I6 S6 \* C9 g0 \$ H8 irem(-3,-2)=-1;mod(-3,-2)=-1;
2 D* ~2 D) y' T) t* M" L* Y
; `. D- w% N/ @- Q9 grem(3,-2)=1;mod(3,-2)=-1;

rem(-3,2)=-1;mod(-3,2)=1;

. X1 v. V% P+ {& l9 Q* p

mod(3,2)=rem(3,2)=1;
: a) [- f$ R/ s8 J' u
" i0 ~6 `7 p* {2 n' B1 h4 @/ emod(3,-2)=rem(3,-2)-2=1-2=-1;
3 x" m1 [! T7 ]( i$ N7 l/ }1 Z5 @
4 j# D7 |$ Q) B8 _, P* j9 C+ @1 xmod(-3,2)=rem(-3,2)+2=-1+2=1;

mudan

madio 发表于 2012-10-30 06:38
' Z/ k( O( A$ Z8 S( P6 nrem这个函数可以换为mod，容易知道一些！
$ o( t, P/ Q# }: T: }1 Zmatlab中的rem()和mod()函数matlab的rem()和mod()函数
rem(x,y) ...

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

wy617958197

7 ~/ _' i5 c8 f. c, k多谢分享资料！