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

mudan

1. for x1=1:13
   \" S. ]% r4 \8 Z0 j
2.     for x2=1:13
   - Z5 ^/ F3 C\" w0 q  h4 f
3.         x3=(30000-4000*x1-3000*x2)/2000;6 J\" [# F, B% \
   
4.         m=13-x3-x1-x2;. b7 d5 u  ~; l\" S7 R, P; k
   
5.         if sign(x3)>=0 &  rem(x3,1)==0 & sign(m)==0  0 U3 h& x: m4 r% a5 C/ N6 `5 m! X
   
6.             disp([x1,x2,x3])
   \" w! f$ o1 _: Q7 l3 w2 t% A. w
7.         end4 @  k. @( j% z# t3 d9 X
   
8.     end
   ' ?\" P1 H2 R3 j  W2 U
9. end

复制代码

wy617958197

; o) p% T& w+ [, R多谢分享资料！

mudan

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

% M2 W# w, B8 T% a% \( a讲的很细致，坛主能不能推荐一本matlab程序实现的书？侧重实例的，谢谢

madio

rem这个函数可以换为mod，容易知道一些！
matlab中的rem()和mod()函数matlab的rem()和mod()函数
7 H5 |% L% g5 n, U. b; K) \rem(x,y):求整除x/y的余数
+ d9 d* h2 m) c6 V
mod(x,y):求模


% u) a, Z3 s# x" X) u
$ u% Y/ I/ L+ _  jrem(x,y)=x-y.*fix(x./y);  (fix()向0取整)
0 i* u5 X9 g# [6 o& r2 |: b
9 l0 R/ `. R% Z9 S) V/ _0 emod(x,y)=x-y.*floor(x./y); （floor()向左取整）
0 L5 h! R3 a6 S- w5 S
; B5 v( R+ f1 ^0 L; T3 s如果x和y的符号相同（同为‘+’，同为‘-’），那么rem(x,y)=mod(x,y)（正数与正数，负数与负数，取整结果两个函数效果一样）

如果x和y的符号相反，那么mod(x,y)=rem(x,y)+y(正数与负数的取整，看你希望得到什么结果就选择用什么函数)
: ]7 S6 g$ z+ y& ^: S* d
当正数与负数取余时，当得到的余数结果的符号希望跟除数(x)一样，用rem()函数；当得到的余数结果的符号希望跟被除数(y)一样，用mod()函数
  C$ j! N% G, H8 e" H" W
/ m4 T3 b$ r: U2 }- X# t例如：

( e+ ^# k2 w3 ]0 F8 b; ~rem(3,2)=1;mod(3,2)=1;

5 F! t  `' O, }, urem(-3,-2)=-1;mod(-3,-2)=-1;

. ^+ c$ J" x/ Z' B7 l- ]* J0 Frem(3,-2)=1;mod(3,-2)=-1;

rem(-3,2)=-1;mod(-3,2)=1;
3 S2 |% z/ d8 J- |
- L. n/ i) d; i5 \% h) i
* `6 n3 P+ I7 r: L) x
mod(3,2)=rem(3,2)=1;

; K2 \% z7 I0 B7 B! Gmod(3,-2)=rem(3,-2)-2=1-2=-1;
/ }$ k! z" }% e- h; J
mod(-3,2)=rem(-3,2)+2=-1+2=1;
  }. [+ a$ T2 T! }; I+ v! a