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

mudan

1. for x1=1:13
   . \9 B$ f- `6 C5 Z\" R6 X
2.     for x2=1:13. y6 V5 K) w. B
   
3.         x3=(30000-4000*x1-3000*x2)/2000;2 o) H+ \8 a8 z- }6 W
   
4.         m=13-x3-x1-x2;
   : V' n* U, ~, j9 ?0 [, ?. L
5.         if sign(x3)>=0 &  rem(x3,1)==0 & sign(m)==0  3 M) ]/ w  d8 z: ~! p6 A6 e6 b8 X% q
   
6.             disp([x1,x2,x3])
   8 l8 K8 _$ Y6 h* s4 M4 B
7.         end6 b( J8 i' c1 N6 n3 o
   
8.     end
   4 Q* }& {4 a' f. {/ S
9. end

复制代码

madio

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

* H+ T" y  i+ D# ~* S9 b2 t* ]# Cmod(x,y):求模
' U3 ~" W. h7 S& X$ ~" w
+ A# k, Z* V5 Z

rem(x,y)=x-y.*fix(x./y);  (fix()向0取整)
; ]  {, ^8 o/ H/ K& g8 E6 D
* R/ z9 ?# Z0 U& {2 f6 xmod(x,y)=x-y.*floor(x./y); （floor()向左取整）

# W6 k5 c+ A: T) l/ w2 l如果x和y的符号相同（同为‘+’，同为‘-’），那么rem(x,y)=mod(x,y)（正数与正数，负数与负数，取整结果两个函数效果一样）
+ w1 O/ G) D7 Y7 V5 P2 R
如果x和y的符号相反，那么mod(x,y)=rem(x,y)+y(正数与负数的取整，看你希望得到什么结果就选择用什么函数)
' G; H' f5 T  G3 A+ {
! l7 c3 T5 y. D$ {当正数与负数取余时，当得到的余数结果的符号希望跟除数(x)一样，用rem()函数；当得到的余数结果的符号希望跟被除数(y)一样，用mod()函数

0 }9 T- ~/ H$ J# q例如：
  |+ N" \. r8 ]. K+ h& X9 p
" F: v5 v  `5 l* y. e! k6 a# Irem(3,2)=1;mod(3,2)=1;

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

rem(3,-2)=1;mod(3,-2)=-1;
$ p4 x# a6 i; i5 S4 p
2 [! J3 @( z- F% o/ rrem(-3,2)=-1;mod(-3,2)=1;
; b/ h. e3 s- j  ?+ M; r9 U' K& [

! R# Q5 g6 p8 m  J
6 X5 G1 U1 u- b9 T' Z  Emod(3,2)=rem(3,2)=1;

- y1 g! f, B3 U# H( b2 ]mod(3,-2)=rem(3,-2)-2=1-2=-1;
3 Y2 V- h6 j7 w5 w9 n3 q
. ^( E8 P( O. L$ f4 @0 N6 C4 Vmod(-3,2)=rem(-3,2)+2=-1+2=1;
5 R4 h2 N6 K: O) \

mudan

madio 发表于 2012-10-30 06:38
: C& {  m6 d2 E7 Arem这个函数可以换为mod，容易知道一些！
matlab中的rem()和mod()函数matlab的rem()和mod()函数
4 ^& l* e6 F: yrem(x,y) ...

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

wy617958197

# t# b$ g' ^& W多谢分享资料！