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

mudan

1. for x1=1:134 _! i/ t- k( ]
   
2.     for x2=1:13% f' P) T! U  z\" @
   
3.         x3=(30000-4000*x1-3000*x2)/2000;
   & t6 V& J4 M$ D2 S. x2 k
4.         m=13-x3-x1-x2;
   * v: L: u( k! [# s- w7 e
5.         if sign(x3)>=0 &  rem(x3,1)==0 & sign(m)==0  
   0 y$ P& ?\" V0 I/ m! i7 D4 D, F' ~
6.             disp([x1,x2,x3])6 H# }( H5 }7 v5 L) \
   
7.         end2 j' ?+ U- B4 V: l5 L
   
8.     end& n6 T8 |8 P: c/ m
   
9. end

复制代码

madio

rem这个函数可以换为mod，容易知道一些！
matlab中的rem()和mod()函数matlab的rem()和mod()函数
- x. |, C0 |. |1 |! `rem(x,y):求整除x/y的余数

$ F, E; I1 O' nmod(x,y):求模


& q( u# `7 A" l$ Q7 t/ g5 [
rem(x,y)=x-y.*fix(x./y);  (fix()向0取整)

mod(x,y)=x-y.*floor(x./y); （floor()向左取整）

如果x和y的符号相同（同为‘+’，同为‘-’），那么rem(x,y)=mod(x,y)（正数与正数，负数与负数，取整结果两个函数效果一样）

如果x和y的符号相反，那么mod(x,y)=rem(x,y)+y(正数与负数的取整，看你希望得到什么结果就选择用什么函数)
! b' z, `. c$ i% ]
+ D$ B6 E8 b$ G4 k7 F% ]当正数与负数取余时，当得到的余数结果的符号希望跟除数(x)一样，用rem()函数；当得到的余数结果的符号希望跟被除数(y)一样，用mod()函数

例如：

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

rem(-3,-2)=-1;mod(-3,-2)=-1;
# @% j) T2 J5 o8 T
, k! T# b3 E6 W* ]rem(3,-2)=1;mod(3,-2)=-1;
. m& e4 h8 X! |* H& q8 D5 T
- N% o, _) |3 b5 Nrem(-3,2)=-1;mod(-3,2)=1;
- ?8 a: P- t2 X( C, V& W3 K

1 G8 B% }9 l' |5 d
mod(3,2)=rem(3,2)=1;
( G0 O7 t  E1 T: s; l
mod(3,-2)=rem(3,-2)-2=1-2=-1;
8 N. a  W( X  Z9 g0 E
mod(-3,2)=rem(-3,2)+2=-1+2=1;

mudan

madio 发表于 2012-10-30 06:38
rem这个函数可以换为mod，容易知道一些！
) X0 J. a% f- F. o! m! Z9 dmatlab中的rem()和mod()函数matlab的rem()和mod()函数
rem(x,y) ...

1 h( k3 F! v" {0 Q9 r. K讲的很细致，坛主能不能推荐一本matlab程序实现的书？侧重实例的，谢谢

wy617958197

多谢分享资料！