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

mudan

1. for x1=1:13$ W* P; H9 l- N; F4 M, t* d
   
2.     for x2=1:13
   , K! `* F\" @1 T1 S8 b  V6 s2 J
3.         x3=(30000-4000*x1-3000*x2)/2000;
   # V1 ?3 x9 J\" T+ X/ T1 M4 R
4.         m=13-x3-x1-x2;
   0 t\" C9 b( z: W+ L/ F3 S# o! y# J
5.         if sign(x3)>=0 &  rem(x3,1)==0 & sign(m)==0  8 q- j6 M: u, A# s; k8 L
   
6.             disp([x1,x2,x3])8 z# G2 q2 N! z* f
   
7.         end
   8 m* a7 I: b3 R
8.     end
   9 u6 S. v& G5 O1 b5 o; c
9. end

复制代码

madio

rem这个函数可以换为mod，容易知道一些！
% h$ X1 G6 l2 T" k' \, xmatlab中的rem()和mod()函数matlab的rem()和mod()函数
$ d% g+ L, X; \7 z& h- ~  A% Q4 krem(x,y):求整除x/y的余数

1 v" |% d1 F" Rmod(x,y):求模
- v) N7 A1 K' d% h9 H


0 G" d8 n& j7 M, r( Urem(x,y)=x-y.*fix(x./y);  (fix()向0取整)
8 u! Z2 D# x9 p% N! y4 T! G$ Z
% z7 p/ e: S) X6 s3 l* [mod(x,y)=x-y.*floor(x./y); （floor()向左取整）

& i' c6 h/ c! b5 t' G. W: `2 g如果x和y的符号相同（同为‘+’，同为‘-’），那么rem(x,y)=mod(x,y)（正数与正数，负数与负数，取整结果两个函数效果一样）
" }. g3 c7 Q$ \) E% d- r
; ]# Y( b- w, Z; [* X; t( a3 g( R如果x和y的符号相反，那么mod(x,y)=rem(x,y)+y(正数与负数的取整，看你希望得到什么结果就选择用什么函数)
. I- H% O( g. o; b" j& W, j
; r  G. R% C1 R6 [9 w; w, W, T当正数与负数取余时，当得到的余数结果的符号希望跟除数(x)一样，用rem()函数；当得到的余数结果的符号希望跟被除数(y)一样，用mod()函数

例如：

* G7 u! c4 s$ N/ _8 n) srem(3,2)=1;mod(3,2)=1;
! u3 j' P) W" s/ K0 ^
rem(-3,-2)=-1;mod(-3,-2)=-1;
4 P2 S& ]$ Q3 ^: S8 |5 q
rem(3,-2)=1;mod(3,-2)=-1;
" {5 a5 y: ]3 n0 J9 }
rem(-3,2)=-1;mod(-3,2)=1;
5 E+ G& _9 n* p) E


mod(3,2)=rem(3,2)=1;
* `8 b0 M% w6 ~# ]. m$ V
6 j. g0 \/ b; Z" o6 Nmod(3,-2)=rem(3,-2)-2=1-2=-1;

# m: K) O# T3 ^( e$ W" Qmod(-3,2)=rem(-3,2)+2=-1+2=1;

mudan

madio 发表于 2012-10-30 06:38
# J: T+ K2 _  d/ [rem这个函数可以换为mod，容易知道一些！
) m. O: Z: v9 U2 H* i& b- y7 bmatlab中的rem()和mod()函数matlab的rem()和mod()函数
5 b3 S2 H9 ]# ]& n1 q+ arem(x,y) ...

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

wy617958197

多谢分享资料！