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

mudan

1. for x1=1:13
   ' E) e8 m/ o1 s( M
2.     for x2=1:13
   , h! W1 Z  i8 J$ k7 h
3.         x3=(30000-4000*x1-3000*x2)/2000;9 c( ?5 O! ]! F' b7 v& p  P
   
4.         m=13-x3-x1-x2;' B8 A% O) z4 O) ?( S3 B
   
5.         if sign(x3)>=0 &  rem(x3,1)==0 & sign(m)==0  2 `5 _! A! e# Q% C, C2 Z! n: T
   
6.             disp([x1,x2,x3])
   ' I9 T: o9 b% R9 D+ g3 a0 `
7.         end
   % |6 T6 a; {6 x/ P& V
8.     end
   : H/ U& k4 t. x
9. end

复制代码

wy617958197

, V3 }, `& e" ]8 d! }; q多谢分享资料！

mudan

madio 发表于 2012-10-30 06:38
! J4 L9 {; v5 |' c0 w, crem这个函数可以换为mod，容易知道一些！
matlab中的rem()和mod()函数matlab的rem()和mod()函数
; c7 X! {1 K" P( i9 H2 Qrem(x,y) ...

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

madio

rem这个函数可以换为mod，容易知道一些！
matlab中的rem()和mod()函数matlab的rem()和mod()函数
rem(x,y):求整除x/y的余数
, x" Z$ C( @$ C  i
mod(x,y):求模
' Y  o8 z  L9 e6 {& _9 D

% C' X: T0 N, j$ Q2 N" i
! a' a/ Y8 z9 |rem(x,y)=x-y.*fix(x./y);  (fix()向0取整)

mod(x,y)=x-y.*floor(x./y); （floor()向左取整）

& k( z% y6 V; I6 e. m/ W7 Q如果x和y的符号相同（同为‘+’，同为‘-’），那么rem(x,y)=mod(x,y)（正数与正数，负数与负数，取整结果两个函数效果一样）
- s+ Q: T2 Z6 \0 m
8 j; ?- K; h" Q  |- D' `$ {如果x和y的符号相反，那么mod(x,y)=rem(x,y)+y(正数与负数的取整，看你希望得到什么结果就选择用什么函数)

当正数与负数取余时，当得到的余数结果的符号希望跟除数(x)一样，用rem()函数；当得到的余数结果的符号希望跟被除数(y)一样，用mod()函数
% u' K) k8 R4 S1 s: ~+ H+ O
8 D8 y# v& G8 l0 _例如：

4 y/ l0 R2 o& p- Prem(3,2)=1;mod(3,2)=1;

rem(-3,-2)=-1;mod(-3,-2)=-1;
# [) p/ \, T$ P" E; h
0 S4 D9 l9 R( Y; k# _' frem(3,-2)=1;mod(3,-2)=-1;
4 X# s7 m8 V  l% k1 V4 D" Z  a
rem(-3,2)=-1;mod(-3,2)=1;
+ x  V+ Y9 g1 E! d8 l1 C8 [2 v" _


mod(3,2)=rem(3,2)=1;
/ q; c! ~3 E6 z2 r% s; \( Z
( b$ D7 A4 [0 K8 a& dmod(3,-2)=rem(3,-2)-2=1-2=-1;

* [% h* p" F3 B1 n$ e, {# xmod(-3,2)=rem(-3,2)+2=-1+2=1;
5 ^0 |8 `/ ]$ G