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

mudan

1. for x1=1:13
   $ t. E( Q7 q. b  l. _
2.     for x2=1:13* E1 k& A1 F) V/ H8 V4 y, ?& |
   
3.         x3=(30000-4000*x1-3000*x2)/2000;
   2 v9 b& x) m+ }1 E, `( E# T
4.         m=13-x3-x1-x2;- W0 _, Q( g0 Y
   
5.         if sign(x3)>=0 &  rem(x3,1)==0 & sign(m)==0  3 v, I* X9 j8 _6 a' E
   
6.             disp([x1,x2,x3])
   8 j5 c7 q  o/ l4 z0 ^
7.         end/ l5 {6 S9 T) `; `' b; W2 Y
   
8.     end
   ! ^7 d* w- J9 T5 a$ k8 s: I& f2 d7 h
9. end

复制代码

madio

rem这个函数可以换为mod，容易知道一些！
matlab中的rem()和mod()函数matlab的rem()和mod()函数
$ r) \/ }  i) Q3 Arem(x,y):求整除x/y的余数
, t2 |# D1 u. ~* B
mod(x,y):求模
6 o8 U) I! v, L! q
) o% q( I2 p( w7 g7 M

. G1 l. Y" L2 u/ Z" z( Prem(x,y)=x-y.*fix(x./y);  (fix()向0取整)
7 a. O5 g7 w2 f$ M+ A
mod(x,y)=x-y.*floor(x./y); （floor()向左取整）
. F4 K/ m$ v. p8 v0 F
7 g+ [5 `6 j5 N& G1 H如果x和y的符号相同（同为‘+’，同为‘-’），那么rem(x,y)=mod(x,y)（正数与正数，负数与负数，取整结果两个函数效果一样）
8 d* T$ c9 E' Q4 `
! F9 s% H3 |$ e# ^如果x和y的符号相反，那么mod(x,y)=rem(x,y)+y(正数与负数的取整，看你希望得到什么结果就选择用什么函数)

+ Q8 B6 d+ }; h/ [" i  i  ?当正数与负数取余时，当得到的余数结果的符号希望跟除数(x)一样，用rem()函数；当得到的余数结果的符号希望跟被除数(y)一样，用mod()函数

例如：

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

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

rem(3,-2)=1;mod(3,-2)=-1;
$ q" R* O# b9 Y( i% t( L4 Z% l
rem(-3,2)=-1;mod(-3,2)=1;


9 H( a9 j# U; G; @$ `. L7 N
0 [8 w  e2 Y1 o" R/ x5 i% e: l  N) X' Tmod(3,2)=rem(3,2)=1;

mod(3,-2)=rem(3,-2)-2=1-2=-1;
5 W  W! X, ?+ y' F8 g8 ?) _
; h5 I* \: [1 A. s8 fmod(-3,2)=rem(-3,2)+2=-1+2=1;
$ h; t& M1 y) {! I! r4 C0 q

mudan

madio 发表于 2012-10-30 06:38
$ d% o' s7 w# w' a" C* A# yrem这个函数可以换为mod，容易知道一些！
. V' a3 V* T% W* u- umatlab中的rem()和mod()函数matlab的rem()和mod()函数
rem(x,y) ...

% {. o- ], r: {8 ?' z; w3 d& K讲的很细致，坛主能不能推荐一本matlab程序实现的书？侧重实例的，谢谢

wy617958197

多谢分享资料！