数学建模社区-数学中国

标题: 求解4000*x1+3000*x2+2000*x3=30000的自然数解程序 [打印本页]

---

作者: mudan    时间: 2012-10-29 16:54
标题: 求解4000*x1+3000*x2+2000*x3=30000的自然数解程序

1. for x1=1:13
   
2.     for x2=1:13
   9 V9 A& J4 A2 `$ s5 K1 a
3.         x3=(30000-4000*x1-3000*x2)/2000;
   
4.         m=13-x3-x1-x2;
   
5.         if sign(x3)>=0 &  rem(x3,1)==0 & sign(m)==0  
   : ~+ ]* I# y2 T& J2 [2 i! Q: O
6.             disp([x1,x2,x3])
   " d: U7 J! E$ I  T
7.         end
   
8.     end
   
9. end

复制代码

---

作者: madio    时间: 2012-10-30 06:38
rem这个函数可以换为mod，容易知道一些！
% R8 B* v) _' T4 t* vmatlab中的rem()和mod()函数matlab的rem()和mod()函数
3 e* B3 V* h1 D6 arem(x,y):求整除x/y的余数

% u% D+ h! ^) ]' Q! w8 o" Ymod(x,y):求模

8 [* F, _; e. m% s
$ Y$ P2 c, J1 t, m! f% ]' ~
rem(x,y)=x-y.*fix(x./y);  (fix()向0取整)

0 g- Q' E  Y( U" [3 P6 qmod(x,y)=x-y.*floor(x./y); （floor()向左取整）

8 R3 d; c# u3 r4 x# X如果x和y的符号相同（同为‘+’，同为‘-’），那么rem(x,y)=mod(x,y)（正数与正数，负数与负数，取整结果两个函数效果一样）
4 A9 }& N! v( e- `6 K; Y
( R; Y) V( q; S- v+ y如果x和y的符号相反，那么mod(x,y)=rem(x,y)+y(正数与负数的取整，看你希望得到什么结果就选择用什么函数)
! X, J; L, {- s/ y* [
当正数与负数取余时，当得到的余数结果的符号希望跟除数(x)一样，用rem()函数；当得到的余数结果的符号希望跟被除数(y)一样，用mod()函数

例如：

: \5 I9 H# M- z/ E& Nrem(3,2)=1;mod(3,2)=1;

# ^3 g6 m) g7 f5 J* X1 A7 zrem(-3,-2)=-1;mod(-3,-2)=-1;

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

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

+ @6 C$ a8 h+ |/ t$ [

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

0 j0 M  {" C$ _- fmod(3,-2)=rem(3,-2)-2=1-2=-1;

) j& J* O9 z3 ]6 jmod(-3,2)=rem(-3,2)+2=-1+2=1;
# J+ `. t# e6 J1 F

---

作者: mudan    时间: 2012-10-30 17:43

madio 发表于 2012-10-30 06:38
6 K5 b8 M# x* v" y& Z2 d5 d+ drem这个函数可以换为mod，容易知道一些！
matlab中的rem()和mod()函数matlab的rem()和mod()函数
rem(x,y) ...

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

---

作者: wy617958197    时间: 2014-9-2 13:07

多谢分享资料！

---

欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)

Powered by Discuz! X2.5