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

mudan

1. for x1=1:13
   & h7 _) Y' m1 B$ Y& q
2.     for x2=1:13
   \" {& M6 G0 d8 a) j, g1 S: E
3.         x3=(30000-4000*x1-3000*x2)/2000;
   8 Y8 |% Z\" P\" `8 s# z, k  J  J
4.         m=13-x3-x1-x2;
   8 B% L( M$ ^* Y
5.         if sign(x3)>=0 &  rem(x3,1)==0 & sign(m)==0  
   $ f& @$ V4 @: v2 M% B
6.             disp([x1,x2,x3])4 X) y' j% x0 ~  R3 l; i
   
7.         end
   ) f; }7 G7 R. M, i4 V7 f/ L
8.     end
   / U5 M* Y: {& H/ q; `
9. end

复制代码

madio

rem这个函数可以换为mod，容易知道一些！
matlab中的rem()和mod()函数matlab的rem()和mod()函数
% y  v. j1 l  u1 b) h4 {rem(x,y):求整除x/y的余数
/ {  `. g" m% R; K; A
3 p  K, |, e8 s2 U" Q7 tmod(x,y):求模

8 X7 {+ X2 s- Q; a" p4 x

rem(x,y)=x-y.*fix(x./y);  (fix()向0取整)
3 n7 m. _# o1 |% a5 q
mod(x,y)=x-y.*floor(x./y); （floor()向左取整）

如果x和y的符号相同（同为‘+’，同为‘-’），那么rem(x,y)=mod(x,y)（正数与正数，负数与负数，取整结果两个函数效果一样）

/ M; B+ T  v3 A0 ?如果x和y的符号相反，那么mod(x,y)=rem(x,y)+y(正数与负数的取整，看你希望得到什么结果就选择用什么函数)

' S- G& S* [& w' z7 Z$ M当正数与负数取余时，当得到的余数结果的符号希望跟除数(x)一样，用rem()函数；当得到的余数结果的符号希望跟被除数(y)一样，用mod()函数

例如：

rem(3,2)=1;mod(3,2)=1;
, E: }5 k% v, X
rem(-3,-2)=-1;mod(-3,-2)=-1;
- g; @# `- z" z- M- F. W' Z# f' P
rem(3,-2)=1;mod(3,-2)=-1;

rem(-3,2)=-1;mod(-3,2)=1;
+ i0 l  ]3 g' W% I& F5 Z, L3 t, N; E
! F! f/ Y% H3 a# @, a
0 _: M* Y+ h5 f& C) X; v7 m& W5 E
mod(3,2)=rem(3,2)=1;

mod(3,-2)=rem(3,-2)-2=1-2=-1;
+ F" \' D& c6 u/ N- b; {* I
mod(-3,2)=rem(-3,2)+2=-1+2=1;
5 R0 U9 E3 j- u5 G9 |/ m2 e

mudan

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

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

wy617958197

多谢分享资料！