数学建模社区-数学中国

标题: matlab编程 [打印本页]

---

作者: E304731541    时间: 2009-8-20 19:47
标题: matlab编程
d=sym(3*x^3-2*x+1);
diff(d,2)
ezplot('cos(x)',[0,pi])//4.2
ezplot('(cos(t))^3','(sin(t))^3',[0,2*pi])//4.3
/ y" B5 g* n$ C( C+ j& X& j" Eezplot('exp(x)+sin(x*y)',[-2,0.5,0,2])//4.4
fminbnd('x+3*(x^2+cos(x))',-1,1)  //直接求
f=sym('x+3*(x^2+cos(x))');//由导数为零来求解
  ff=solve(diff(f));
( f- l- v* j6 W& z0 Z, Q* Tsubs(f,x,ff),subs(f,x,-1),subs(f,x,-1)
syms x;
1 R& |  V+ a; z$ i8 g7 b8 \int(1/sin(x)^2*cos(x)^2)
5 f7 R4 ~' P/ jsyms x z;
$ A+ t" i! t) k, p4 B& Yint(x/(1+z^2),x)
+ ^; h1 a4 p. s9 S4 o7 KA=[1,2,3;4,5,6;7,8,9;10,11,12]
) N/ t( \/ }+ B, \1 ?A =
: z8 Q6 _: K, w2 I     1     2     3
     4     5     6
     7     8     9
) i% Y5 H4 y1 c6 h, Q5 q) F# H8 f    10    11    12
B=[0,1,1;2,2,1;1,1,1];
; q6 N- h9 j2 GA*B
A^3//错误，矩阵的幂必须矩阵的行列相同才行
rref(A)
# t5 |- W- X! O- u" n% Tx=linspace(0,2*pi,30);
y=sin(x);
z=cos(x);
plot(x,y,'r',x,z,'pg')
8 g" t) E- Q/ H. J' h7 e2 Ux=linspace(0,2*pi,30);
0 d. u% I" n1 g% E$ Y: My=sin(x);
/ l( {( G) Z+ n  v: @z=cos(x);
% h& H' V6 e1 t& W" F+ f/ pxlabel('自变量x')
8 g! R4 }; Q8 F1 b$ Bylabel('函数y')
title('示意图')
* @8 ^8 r: i4 o+ ?9 b( V2 Q  pplot(x,y,'r',x,z,'pg'),grid on;
, I! n9 E/ l9 @- e$ m8 Ex=linspace(0,2*pi,30);
: c6 h+ T! p% l! I6 ey=sin(x);
z=cos(x);
  x" P$ L- Z% r5 Y! \, {plot(x,z,'pg');
- @0 _8 m9 F6 A8 a4 E+ ?hold on;
! @1 ~1 P4 t* j, D, t& [1 [+ f8 \plot(x,y,'r')
x=linspace(0,2*pi,100);
- @: W# w0 O! [' b  \y=sin(x);
! A' D* A6 Y' _7 ?. \# Az=cos(x);
plot(x,y);
! A0 W. h+ l) U1 dtitle('sin(x)')
+ p! f& [3 Q. J2 C& u4 v# |$ lpause;
% W" }* ?# P$ y$ v3 G& z! Q* B& Qfigure(2);
plot(x,z,'pg');
title('cos(x)')
+ c7 Z4 w& F% {x=linspace(0,2*pi,30);
y=sin(x);
plot(x,y,'r');
zoom on
p=[1 -4 9 -10];
& K# k. Q( t' f9 zr=roots(p)//次数小于等于4，若大于4则用该命令不一定能求出所有根的准确形式。
r=solve('x^3-4*x^2+9*x-10=0')//简单方法
! E3 ?8 [  y+ Y8 S7 `" H2 l; Ur=solve('x^2-a*x-4*b=0','x')
; I' p0 b/ ]. N/ Xd=solve('p*sin(x)-r','x')
6 s' o9 F/ r' l1 m8 n( Z# P3 x( Tr=fzero('x-(cos(x))^2',1)//其中必须输入的是方程对应的函数名
+ q7 Y! J% k5 ?: r[x,y]=solve('x+y=1','x-11*y=5','x','y')

---

作者: mrx    时间: 2009-8-21 06:07
matlab 里的主成分分析.rar (864 Bytes, 下载次数: 4)

2009-8-21 06:07 上传

点击文件名下载附件
下载积分: 体力 -2 点

关于校车安排问题.rar

2009-8-21 06:07 上传

点击文件名下载附件

下载积分: 体力 -2 点

669 Bytes, 下载次数: 4, 下载积分: 体力 -2 点

火车票全程对号销售系统优化研究.rar

2009-8-21 06:07 上传

点击文件名下载附件

下载积分: 体力 -2 点

634 Bytes, 下载次数: 6, 下载积分: 体力 -2 点

金融风暴下快餐连锁企业的促销推广策略.rar

2009-8-21 06:07 上传

点击文件名下载附件

下载积分: 体力 -2 点

1.05 KB, 下载次数: 7, 下载积分: 体力 -2 点

蚁群算法的基本原理.rar

2009-8-21 06:07 上传

点击文件名下载附件

下载积分: 体力 -2 点

2.55 KB, 下载次数: 8, 下载积分: 体力 -2 点

---

作者: qingqing20072    时间: 2009-8-21 10:28
hahahahahahhahahahhao a hao a

---

作者: qingqing20072    时间: 2009-8-21 10:29
ahaoahhahahahhahahhahahah

---

作者: qingqing20072    时间: 2009-8-21 10:32
hahahahahahhaaaaaaaaaaaa

---

作者: xxj198712    时间: 2009-8-31 14:12
hahahahahahhaaaaaaaaaaaa

---

作者: zhula    时间: 2009-8-31 14:18
2# mrx
谢谢啊！！

---

作者: zhula    时间: 2009-8-31 14:20
2# mrx
太感谢了

---

作者: zhula    时间: 2009-8-31 14:21
2# mrx
$ Z) c* {3 S$ k0 \2 e最近在做有关促销的文章，可以借鉴，３Ｑ

---

作者: byit    时间: 2009-8-31 14:30
Thanks you~~

---

作者: 阿KEN    时间: 2009-8-31 14:31
不太懂？？？？

---

作者: 遗忘    时间: 2009-8-31 14:57
非常感谢啊！

---

作者: 248641018    时间: 2009-8-31 21:54
真的很不错哦~~这些程序解了燃眉之急啊！！

---

作者: 248641018    时间: 2009-8-31 22:01
哎呀没有矩阵币呀~~~

---

作者: amlen2007    时间: 2009-8-31 22:43
不错。。。。。。。。。

---

作者: cool_gr    时间: 2009-9-1 15:21
金融风暴是个热点啊

---

作者: nicecatbox    时间: 2009-9-1 20:42
不错不错~~

---

作者: 12as12as    时间: 2009-9-1 21:05
太感谢了THANK

---

作者: alair005    时间: 2012-2-6 15:46
支持一,下楼主辛苦了5891630098609316

---

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

Powered by Discuz! X2.5