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

mudan

1. for x1=1:13% ]\" T& ]8 b2 x+ @
   
2.     for x2=1:13: v3 C5 X- T& ?) B
   
3.         x3=(30000-4000*x1-3000*x2)/2000;
   ' N3 o' g& l; f9 ^# C; {9 r, K6 _  C
4.         m=13-x3-x1-x2;2 r; j9 ~% v1 Z% p8 ?# K
   
5.         if sign(x3)>=0 &  rem(x3,1)==0 & sign(m)==0  
   * d\" L. p, G7 J* E4 B: V4 `% E
6.             disp([x1,x2,x3])
   9 S' R+ G6 R) ?( Q: x. j% Z
7.         end
   * @' v% n2 }# `( i1 q' ]
8.     end\" f# {% M0 u: N; B. I5 l( G
   
9. end

复制代码

madio

rem这个函数可以换为mod，容易知道一些！
- }0 I. X2 D" u3 |matlab中的rem()和mod()函数matlab的rem()和mod()函数
/ q3 j3 F' {3 ~' K$ z0 h- |4 Arem(x,y):求整除x/y的余数
1 _# U2 C# w  _- b6 H8 t
mod(x,y):求模
4 @$ ^) S# J8 v( b' I3 [( s


rem(x,y)=x-y.*fix(x./y);  (fix()向0取整)

" P% v9 |( ~) omod(x,y)=x-y.*floor(x./y); （floor()向左取整）

如果x和y的符号相同（同为‘+’，同为‘-’），那么rem(x,y)=mod(x,y)（正数与正数，负数与负数，取整结果两个函数效果一样）
! u3 P9 B, }- j- O+ U% F( b
" _4 y* q  j4 F如果x和y的符号相反，那么mod(x,y)=rem(x,y)+y(正数与负数的取整，看你希望得到什么结果就选择用什么函数)
+ {  [. H1 x0 }- @% V# b5 G
当正数与负数取余时，当得到的余数结果的符号希望跟除数(x)一样，用rem()函数；当得到的余数结果的符号希望跟被除数(y)一样，用mod()函数
* `: K  P7 X/ q
3 |7 ?. ~# ]5 Z' e例如：
0 F  ?- V9 ?. c0 A$ L4 }4 J
* }# p4 V' V7 s& r* g* {5 [% Krem(3,2)=1;mod(3,2)=1;
7 A3 C( F+ Y" C4 U, _! e- Q) X& Q, P
3 w! [* ], _6 K; {7 I% urem(-3,-2)=-1;mod(-3,-2)=-1;

1 \. z3 l% e, j  Xrem(3,-2)=1;mod(3,-2)=-1;

1 @# @, F1 r" t+ Q7 P0 Z) krem(-3,2)=-1;mod(-3,2)=1;
; o& X/ V# ]. R" g4 e1 t

  i9 `* w0 e1 a9 S# ~0 s3 L
. p/ ?7 Q; V: J# Cmod(3,2)=rem(3,2)=1;
2 C* @8 w& H! q. g. g( A
mod(3,-2)=rem(3,-2)-2=1-2=-1;

mod(-3,2)=rem(-3,2)+2=-1+2=1;
8 ^% ^4 s$ G+ ^1 o! p8 t& E

mudan

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

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

wy617958197

4 t# x/ y, n; ~+ E4 r; ^7 u# j多谢分享资料！