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

mudan

1. for x1=1:131 Y' Y, m& Z+ J! [% h\" {6 ^* W: h
   
2.     for x2=1:13
   \" g' d3 y* V2 F# ~
3.         x3=(30000-4000*x1-3000*x2)/2000;
   4 M7 p: F; Y# X, e
4.         m=13-x3-x1-x2;! b' Q1 j  R' A& m6 j/ M% G
   
5.         if sign(x3)>=0 &  rem(x3,1)==0 & sign(m)==0  
   . m0 V1 T& O- w+ U7 X; c
6.             disp([x1,x2,x3])
   + ^) g. u, }- H1 b) E7 q2 q7 U
7.         end1 w1 @) d# ]$ O( S/ U4 C! ]
   
8.     end
   $ M9 y, G( M+ f* W
9. end

复制代码

madio

rem这个函数可以换为mod，容易知道一些！
- f8 I' \: r5 R' C$ L8 u1 Cmatlab中的rem()和mod()函数matlab的rem()和mod()函数
rem(x,y):求整除x/y的余数

- Q, O# i& Z& a# _' w+ Fmod(x,y):求模
8 l, m9 |4 b* |3 ^; p+ [6 F
! Q2 [! x) N2 i5 G

rem(x,y)=x-y.*fix(x./y);  (fix()向0取整)
+ |$ Q2 B% n9 l) C: g2 I# M  u
mod(x,y)=x-y.*floor(x./y); （floor()向左取整）
, w; x) {+ j# s! j/ [& Y) H; @
9 P- r0 R( g8 G9 Y如果x和y的符号相同（同为‘+’，同为‘-’），那么rem(x,y)=mod(x,y)（正数与正数，负数与负数，取整结果两个函数效果一样）

" r5 b/ D# D$ n如果x和y的符号相反，那么mod(x,y)=rem(x,y)+y(正数与负数的取整，看你希望得到什么结果就选择用什么函数)
4 M; t2 R- _" R+ t7 r: k$ u. o
当正数与负数取余时，当得到的余数结果的符号希望跟除数(x)一样，用rem()函数；当得到的余数结果的符号希望跟被除数(y)一样，用mod()函数

例如：

" f4 t0 q7 @+ t3 _rem(3,2)=1;mod(3,2)=1;
& b5 o$ A( A% D# Q( H
% H5 G  F. Y& \" h! N4 @9 Rrem(-3,-2)=-1;mod(-3,-2)=-1;

; e1 D! u; E7 ~: W1 wrem(3,-2)=1;mod(3,-2)=-1;

rem(-3,2)=-1;mod(-3,2)=1;
: D  D8 y& W2 g% i* d: b

/ O. g8 J+ C0 c: Q% v1 a
mod(3,2)=rem(3,2)=1;
5 P  R3 s4 [0 N" @( ]4 G# c
5 W" u* y, t7 r* amod(3,-2)=rem(3,-2)-2=1-2=-1;

6 _7 h- H5 {; B7 ?# |# kmod(-3,2)=rem(-3,2)+2=-1+2=1;

mudan

madio 发表于 2012-10-30 06:38
/ J# G+ z! T. \" S& n9 ~rem这个函数可以换为mod，容易知道一些！
1 R/ G7 i: y+ ematlab中的rem()和mod()函数matlab的rem()和mod()函数
1 D# ^+ l- O0 o! Mrem(x,y) ...

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

wy617958197

8 t! s: u" n/ D3 L: M多谢分享资料！