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

mudan

1. for x1=1:137 N  H  I$ C9 d& u$ N7 _
   
2.     for x2=1:13\" l* @* E0 ?6 U8 s- q
   
3.         x3=(30000-4000*x1-3000*x2)/2000;
   5 I\" F& T- q8 m2 m
4.         m=13-x3-x1-x2;
   & x9 s7 f\" w4 E6 `0 p
5.         if sign(x3)>=0 &  rem(x3,1)==0 & sign(m)==0  
   / n: a+ y- p, A$ m/ I/ L6 O
6.             disp([x1,x2,x3])
   , T* S: Z$ q0 f* C3 q8 k0 D/ L
7.         end
   ( l- g* p4 _* g! v' O9 d
8.     end9 Q6 [2 z  u- x4 a9 O
   
9. end

复制代码

madio

rem这个函数可以换为mod，容易知道一些！
matlab中的rem()和mod()函数matlab的rem()和mod()函数
rem(x,y):求整除x/y的余数
) B: K$ @( G( {
mod(x,y):求模
( e9 N6 m1 {; P; @; _


rem(x,y)=x-y.*fix(x./y);  (fix()向0取整)
4 A" l5 C" b$ O  u5 |2 C
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(正数与负数的取整，看你希望得到什么结果就选择用什么函数)

当正数与负数取余时，当得到的余数结果的符号希望跟除数(x)一样，用rem()函数；当得到的余数结果的符号希望跟被除数(y)一样，用mod()函数

例如：

rem(3,2)=1;mod(3,2)=1;
1 F4 C" y" y; \; c; J" l
rem(-3,-2)=-1;mod(-3,-2)=-1;
' J: o& a, S; C3 l- J1 W8 B9 S' s
rem(3,-2)=1;mod(3,-2)=-1;
% v  r  c9 l* ^3 Z5 c: w0 ?1 j$ [
rem(-3,2)=-1;mod(-3,2)=1;
% }7 a- M* J0 i2 m$ W( b/ O1 j& A$ o. m
$ D3 x' L1 q) i4 D8 Z

$ s+ G3 ?3 |' R% f" \& w+ ~4 F7 Omod(3,2)=rem(3,2)=1;

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

- E2 r. i& m0 d( z) e9 g& l6 kmod(-3,2)=rem(-3,2)+2=-1+2=1;
* g. k# x5 ~) g0 V0 P) b! ^$ M

mudan

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

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

wy617958197

( U  @5 i' f: [8 x5 K多谢分享资料！