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

mudan

1. for x1=1:13
   , e) [! ^$ b/ Y2 H+ ^2 K% O/ @. T
2.     for x2=1:138 w! O* E, g8 R) D* D
   
3.         x3=(30000-4000*x1-3000*x2)/2000;* P. E3 P$ A# R7 |5 w
   
4.         m=13-x3-x1-x2;$ H& q7 {0 V' a; u7 _$ Z
   
5.         if sign(x3)>=0 &  rem(x3,1)==0 & sign(m)==0  \" J& g; o9 ~. g; f) B
   
6.             disp([x1,x2,x3])4 Y. N* Y$ J6 D6 C: n
   
7.         end
   $ O* N/ G$ C- [  Y1 x
8.     end/ Z9 v* l$ b* i9 ^: A* l\" Y
   
9. end

复制代码

madio

rem这个函数可以换为mod，容易知道一些！
matlab中的rem()和mod()函数matlab的rem()和mod()函数
/ i% ]# H5 a& q. p) lrem(x,y):求整除x/y的余数
( C+ o' v, H& U* R& A5 X2 m
mod(x,y):求模

' J0 u$ {" g  Y9 O( `
$ P& q2 n! K; K( i# W- u
+ p- M( l# |# Crem(x,y)=x-y.*fix(x./y);  (fix()向0取整)
+ ]$ X: e. ?9 E
mod(x,y)=x-y.*floor(x./y); （floor()向左取整）

4 j1 f" b+ i. T" q% w- |+ Q1 b0 E& e' l如果x和y的符号相同（同为‘+’，同为‘-’），那么rem(x,y)=mod(x,y)（正数与正数，负数与负数，取整结果两个函数效果一样）
& D! R2 y- }. K
如果x和y的符号相反，那么mod(x,y)=rem(x,y)+y(正数与负数的取整，看你希望得到什么结果就选择用什么函数)

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

( N$ s7 G, b" g$ d7 `. S: E例如：

4 q$ K/ J7 D2 ~/ u0 Prem(3,2)=1;mod(3,2)=1;

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

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

rem(-3,2)=-1;mod(-3,2)=1;
& @" m' x, R+ l% ]4 j; U) s0 t


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

' z  }9 M) _) N8 Lmod(3,-2)=rem(3,-2)-2=1-2=-1;

( i' B+ `( L0 P1 u  |2 G" E+ zmod(-3,2)=rem(-3,2)+2=-1+2=1;
/ Q% u9 C' A1 k: s

mudan

madio 发表于 2012-10-30 06:38
rem这个函数可以换为mod，容易知道一些！
matlab中的rem()和mod()函数matlab的rem()和mod()函数
rem(x,y) ...

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

wy617958197

& v) ]1 C% A6 V) x" t8 M, b多谢分享资料！