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

mudan

1. for x1=1:13\" X- V2 n/ M9 b$ J9 _
   
2.     for x2=1:13  w! R: h: D/ J3 D
   
3.         x3=(30000-4000*x1-3000*x2)/2000;
   % R5 j( G4 I5 R9 ?' c
4.         m=13-x3-x1-x2;
   $ S# h. e! e3 }; H
5.         if sign(x3)>=0 &  rem(x3,1)==0 & sign(m)==0  
   ( W+ G: r& \' n* N
6.             disp([x1,x2,x3])/ G# _: v- q+ a
   
7.         end. N' ~: c2 C# J! B3 d& ^# F
   
8.     end# K6 f' u5 y/ J
   
9. end

复制代码

wy617958197

多谢分享资料！

mudan

madio 发表于 2012-10-30 06:38
rem这个函数可以换为mod，容易知道一些！
% I. z  A$ b( C" {+ C; V9 Bmatlab中的rem()和mod()函数matlab的rem()和mod()函数
rem(x,y) ...

+ x" g6 V& P0 @7 g讲的很细致，坛主能不能推荐一本matlab程序实现的书？侧重实例的，谢谢

madio

rem这个函数可以换为mod，容易知道一些！
matlab中的rem()和mod()函数matlab的rem()和mod()函数
rem(x,y):求整除x/y的余数
2 M6 i5 V0 F, U7 b5 R
+ h2 F! z! Y- v* L- {0 o: dmod(x,y):求模
% l, H* F- m7 W. r% [+ [# {
- c  X0 j1 H9 [2 Y! }+ s, a7 _0 E! w+ C

. v3 K$ p2 f$ f" t2 t6 vrem(x,y)=x-y.*fix(x./y);  (fix()向0取整)

/ k: J9 L6 U+ }; ]( D( nmod(x,y)=x-y.*floor(x./y); （floor()向左取整）
: R4 L6 u. k" _4 e% H+ J( \
如果x和y的符号相同（同为‘+’，同为‘-’），那么rem(x,y)=mod(x,y)（正数与正数，负数与负数，取整结果两个函数效果一样）

! ?/ V) A2 l& Y9 K5 V* b* G如果x和y的符号相反，那么mod(x,y)=rem(x,y)+y(正数与负数的取整，看你希望得到什么结果就选择用什么函数)
0 ^4 g7 h, Y0 E6 b" h5 V5 B2 k
当正数与负数取余时，当得到的余数结果的符号希望跟除数(x)一样，用rem()函数；当得到的余数结果的符号希望跟被除数(y)一样，用mod()函数
- ?% H3 j( H5 Y: s8 E
例如：
% k. n, j) o$ K7 I5 B2 p% v
rem(3,2)=1;mod(3,2)=1;
  e* ^% C4 w6 A! H$ n; h
9 t$ G8 S3 ^0 J8 Wrem(-3,-2)=-1;mod(-3,-2)=-1;
) r4 P' _9 r# k3 d; I: _3 h
rem(3,-2)=1;mod(3,-2)=-1;

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

% C" D! P- A( c- X$ s( o

' N( d4 y# _3 Bmod(3,2)=rem(3,2)=1;

mod(3,-2)=rem(3,-2)-2=1-2=-1;
3 y. [5 m3 q$ ~
mod(-3,2)=rem(-3,2)+2=-1+2=1;
4 I+ \) }& t5 u9 _$ r5 M