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

mudan

1. for x1=1:13
   9 n1 h, B: T7 p
2.     for x2=1:13
   2 Q2 c& M# I& X7 x! ^  h# j
3.         x3=(30000-4000*x1-3000*x2)/2000;
   / L. e/ {\" u: M3 r\" A) X
4.         m=13-x3-x1-x2;6 \: a5 V$ u. k( i1 ^2 h% Q( s6 k
   
5.         if sign(x3)>=0 &  rem(x3,1)==0 & sign(m)==0  / Y- \8 s# F! q* p7 }. v6 q
   
6.             disp([x1,x2,x3])/ X* Q  ~7 e' t( A\" a% b
   
7.         end0 f0 D- B- b& |% M
   
8.     end2 S+ ]. U. K$ j. h
   
9. end

复制代码

wy617958197

3 Q& P$ U  F  y/ z多谢分享资料！

mudan

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

) i5 R$ u- a  T, Y7 J; K& J; s' y讲的很细致，坛主能不能推荐一本matlab程序实现的书？侧重实例的，谢谢

madio

rem这个函数可以换为mod，容易知道一些！
matlab中的rem()和mod()函数matlab的rem()和mod()函数
rem(x,y):求整除x/y的余数
& L: |! _0 [0 H3 u+ s
mod(x,y):求模


: M6 O5 K1 A1 t; G9 k& A
! _$ I0 E; u& ^, E$ [* Zrem(x,y)=x-y.*fix(x./y);  (fix()向0取整)
$ I+ s3 i. l, y; z7 S
. Q6 R/ h) u6 c0 `' w# Wmod(x,y)=x-y.*floor(x./y); （floor()向左取整）
' v$ K& o/ c) l+ z
如果x和y的符号相同（同为‘+’，同为‘-’），那么rem(x,y)=mod(x,y)（正数与正数，负数与负数，取整结果两个函数效果一样）
9 w% l# @4 P' |
如果x和y的符号相反，那么mod(x,y)=rem(x,y)+y(正数与负数的取整，看你希望得到什么结果就选择用什么函数)

; l: I: N: u) {  q当正数与负数取余时，当得到的余数结果的符号希望跟除数(x)一样，用rem()函数；当得到的余数结果的符号希望跟被除数(y)一样，用mod()函数
' W+ S2 f, d& ?8 l% P6 x( t
例如：

* x' o* N2 I- h7 s! Q; p7 S2 A  Qrem(3,2)=1;mod(3,2)=1;
! N7 z0 `& B2 s2 [4 H
rem(-3,-2)=-1;mod(-3,-2)=-1;

rem(3,-2)=1;mod(3,-2)=-1;
6 V  B3 B! ?9 e* _. v; W
rem(-3,2)=-1;mod(-3,2)=1;


8 t4 ]  U: b: p$ Q
$ g0 I( S5 s! j) W1 vmod(3,2)=rem(3,2)=1;

mod(3,-2)=rem(3,-2)-2=1-2=-1;
4 u) L) C" h- _/ o3 D
mod(-3,2)=rem(-3,2)+2=-1+2=1;
: d( m' C4 Q! \' B5 P3 M