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

mudan

1. for x1=1:13
   ) h% T; w! B+ |\" p3 `8 _
2.     for x2=1:13  N2 a7 J) w) s  z) n6 p
   
3.         x3=(30000-4000*x1-3000*x2)/2000;
   $ N2 @- f) C8 B
4.         m=13-x3-x1-x2;
   - b% f! T/ ?; X& h1 @7 l  P4 h. [
5.         if sign(x3)>=0 &  rem(x3,1)==0 & sign(m)==0  
   $ |0 F) l; c; v9 R- {# F7 f$ g
6.             disp([x1,x2,x3])
   + v) m0 o4 O- s0 p
7.         end4 W5 D& N+ j) E* T
   
8.     end& G' U0 s5 A  Z2 J: }( `: y; Q# O
   
9. end

复制代码

madio

rem这个函数可以换为mod，容易知道一些！
5 c" @2 Q0 A8 n3 c4 Z+ Y  I+ ]matlab中的rem()和mod()函数matlab的rem()和mod()函数
rem(x,y):求整除x/y的余数

mod(x,y):求模



rem(x,y)=x-y.*fix(x./y);  (fix()向0取整)

* `! T/ T2 A8 e; C7 K$ e, z, Wmod(x,y)=x-y.*floor(x./y); （floor()向左取整）
" {$ U8 J% [5 `) a* o
/ X& F& h8 r) G" @如果x和y的符号相同（同为‘+’，同为‘-’），那么rem(x,y)=mod(x,y)（正数与正数，负数与负数，取整结果两个函数效果一样）

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

: s" v% x& S+ k. j+ L# ^7 G& S当正数与负数取余时，当得到的余数结果的符号希望跟除数(x)一样，用rem()函数；当得到的余数结果的符号希望跟被除数(y)一样，用mod()函数
5 w: L' c5 Y/ B/ z, d& h
9 n9 {% g1 k6 A: ~( }) B  K2 h& A1 y例如：

9 w5 E  Z  e. Hrem(3,2)=1;mod(3,2)=1;

) O1 d& B, t; Brem(-3,-2)=-1;mod(-3,-2)=-1;

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

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


5 g( K- ^( F6 }& D
' U: ?, b8 i) v: t5 qmod(3,2)=rem(3,2)=1;

mod(3,-2)=rem(3,-2)-2=1-2=-1;
# K+ x6 K# K! }( b) W0 M
7 a! P: f7 s- Wmod(-3,2)=rem(-3,2)+2=-1+2=1;
1 t1 C7 D* r$ j' y

mudan

madio 发表于 2012-10-30 06:38
rem这个函数可以换为mod，容易知道一些！
9 G1 i4 ?2 J5 bmatlab中的rem()和mod()函数matlab的rem()和mod()函数
+ n$ |5 M1 e' t- _rem(x,y) ...

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

wy617958197

" i3 |7 r' g6 e* S; j: d5 V. I多谢分享资料！