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

mudan

1. for x1=1:13
   # [& g* a! r5 x; |/ l) Z4 M! Y
2.     for x2=1:13
   \" @  n! Q3 @6 P' Y\" c$ m5 h  r
3.         x3=(30000-4000*x1-3000*x2)/2000;
     U/ O. y1 D$ i% i: |
4.         m=13-x3-x1-x2;
   8 n: X: o: }; {; F# f6 Y8 c# k
5.         if sign(x3)>=0 &  rem(x3,1)==0 & sign(m)==0  
   % Q  Z% B) i7 \, X1 a/ ^
6.             disp([x1,x2,x3])
   & @; _; X/ ?% ~& @7 d6 l+ L9 a/ ^$ d
7.         end
   : K% `% z9 M7 W2 h5 H, R) k, A% ?
8.     end% K) q$ w' c* T
   
9. end

复制代码

madio

rem这个函数可以换为mod，容易知道一些！
matlab中的rem()和mod()函数matlab的rem()和mod()函数
: h% V' ?) L6 R+ _1 n5 Nrem(x,y):求整除x/y的余数
0 S: z. V4 y( e2 S" z2 S. y
mod(x,y):求模

0 m/ T4 `' u/ ^& _

; G1 a9 X9 P! X) U. D- q8 ~3 `$ `rem(x,y)=x-y.*fix(x./y);  (fix()向0取整)

( v3 }4 k% H5 Q% Emod(x,y)=x-y.*floor(x./y); （floor()向左取整）
. C6 R) o' s6 \: \
" r4 \5 B3 n- ?3 w如果x和y的符号相同（同为‘+’，同为‘-’），那么rem(x,y)=mod(x,y)（正数与正数，负数与负数，取整结果两个函数效果一样）
8 V3 M9 ^! [. n$ D1 h
如果x和y的符号相反，那么mod(x,y)=rem(x,y)+y(正数与负数的取整，看你希望得到什么结果就选择用什么函数)

当正数与负数取余时，当得到的余数结果的符号希望跟除数(x)一样，用rem()函数；当得到的余数结果的符号希望跟被除数(y)一样，用mod()函数
- p( Q' x& P9 k& E9 x& g  K  r1 }
9 l8 q) h% \; G) ~  H. y) L6 S9 g例如：
- H1 L- F9 `6 g( j; w
4 O: G) _$ N! F# G7 B2 Brem(3,2)=1;mod(3,2)=1;
5 N* D: V& r, f+ Q% L" ?
rem(-3,-2)=-1;mod(-3,-2)=-1;

rem(3,-2)=1;mod(3,-2)=-1;
8 k# [* ]7 `! K& {, k0 q- r$ l! ^
& B2 }6 U+ a3 e: n5 ?1 Xrem(-3,2)=-1;mod(-3,2)=1;


" ]# `/ c/ p9 B) E  v
mod(3,2)=rem(3,2)=1;
3 K+ O7 ?' t+ j5 I! \+ o3 O1 _
# q% j9 F/ m; [- Mmod(3,-2)=rem(3,-2)-2=1-2=-1;
" a6 L6 |3 o( l7 C- A( b( P
* e1 [3 Z& P+ s$ \" umod(-3,2)=rem(-3,2)+2=-1+2=1;

mudan

madio 发表于 2012-10-30 06:38
  a& i% X' ^2 }rem这个函数可以换为mod，容易知道一些！
matlab中的rem()和mod()函数matlab的rem()和mod()函数
2 `2 k- y+ c- Y$ X1 o0 ?7 zrem(x,y) ...

3 g3 H4 q* Z' r讲的很细致，坛主能不能推荐一本matlab程序实现的书？侧重实例的，谢谢

wy617958197

多谢分享资料！