QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 5680|回复: 18
打印 上一主题 下一主题

[代码资源] matlab编程

[复制链接]
字体大小: 正常 放大

5

主题

5

听众

270

积分

升级  85%

该用户从未签到

跳转到指定楼层
1#
发表于 2009-8-20 19:47 |只看该作者 |倒序浏览
|招呼Ta 关注Ta
d=sym(3*x^3-2*x+1);4 ?2 i/ X5 H* z& s* ]
diff(d,2)
5 R2 k+ \) T" ?" W2 X  @& M, {7 Uezplot('cos(x)',[0,pi])//4.2* }! R. g0 ~- z1 W4 m
ezplot('(cos(t))^3','(sin(t))^3',[0,2*pi])//4.3
( l, @* U% K$ R7 Y' \ezplot('exp(x)+sin(x*y)',[-2,0.5,0,2])//4.4- D) X+ _; A0 i$ k9 Q/ b
fminbnd('x+3*(x^2+cos(x))',-1,1)  //直接求& N" y; I/ S' [! v2 R. [
f=sym('x+3*(x^2+cos(x))');//由导数为零来求解& y& ~/ r; n9 {" t" Q/ A
  ff=solve(diff(f));% e3 {1 m  A0 b  e- @2 y, t
subs(f,x,ff),subs(f,x,-1),subs(f,x,-1)
% Q  l" f" A8 j  N5 o3 K% U  Jsyms x;) v1 g6 N# J9 t) C
int(1/sin(x)^2*cos(x)^2)2 M/ _9 u  N+ E8 }* I  a+ r
syms x z;
# C: S( a1 b; X/ X3 D; |int(x/(1+z^2),x)
  |8 }! ]/ ?; r, G$ w3 k9 JA=[1,2,3;4,5,6;7,8,9;10,11,12]1 v9 `2 z( F* `1 D
A =$ G" i2 l+ Q7 t  `+ k9 Z) o
     1     2     3& ~0 l* Z1 O$ Q& E" b' w# b' b
     4     5     6
8 l7 z2 V0 ^+ }     7     8     9
/ f$ M, }6 K- G    10    11    125 X, u* o+ _* F
B=[0,1,1;2,2,1;1,1,1];! S! y/ f0 V! q) b# w1 B3 A. v2 O+ i
A*B5 X* c( y: T3 L# w8 U; O8 w
A^3//错误,矩阵的幂必须矩阵的行列相同才行# u8 q+ o. x. V
rref(A)
9 l$ O3 a) ]8 n. T. W# N# Ox=linspace(0,2*pi,30);
1 {+ t1 z  ?. p& i; Wy=sin(x);& g/ V5 s) v2 C9 E# k; z# K
z=cos(x);! x9 w4 a* ^' y! M
plot(x,y,'r',x,z,'pg')
+ a0 @9 l' |$ c" F* C% T4 b. }$ xx=linspace(0,2*pi,30);
6 h3 P0 t" X3 [2 z9 k6 dy=sin(x);
/ p4 e: v4 L' m* o& ?: t3 jz=cos(x);
9 X: ^( `  J4 ^9 {" ~xlabel('自变量x')
1 F3 v3 q+ T1 L8 ]ylabel('函数y')
' S; C- D4 H# I6 V  Mtitle('示意图')3 R; t; A" K! |% g
plot(x,y,'r',x,z,'pg'),grid on;6 ~: n7 Y$ p! }( X
x=linspace(0,2*pi,30);- n; P# W: U- r8 ?5 M/ D* b. ?
y=sin(x);/ }9 G. n8 v6 h
z=cos(x);2 b& e0 l; T% R8 s6 d
plot(x,z,'pg');
8 R6 ], N# p3 c+ jhold on;
7 }# Y+ W$ l5 D1 K% x6 ?3 {) bplot(x,y,'r')' T3 q: D1 }7 n, Z$ O' e
x=linspace(0,2*pi,100);
* e1 A$ c: C# e  }7 V3 jy=sin(x);
2 u% R3 c9 {: Z1 ez=cos(x);
7 q3 r0 A2 x( q4 mplot(x,y);
0 I  l" ~8 \( {$ o' P1 `3 ~9 g, Ctitle('sin(x)')
) Z+ K6 O* }  @/ r: o: w! S; f, m+ Fpause;
: ?0 j0 M$ U) S4 ifigure(2);
7 ?9 R! ?* e5 ]' r( c: G# L8 M6 V8 nplot(x,z,'pg');% s' d4 G# ~0 e: z
title('cos(x)')5 F4 ?8 V+ e, R# u9 H, ^' t$ M6 c
x=linspace(0,2*pi,30);
% p' E' x" P5 M& T# y, zy=sin(x);
0 L+ U5 q0 f4 Bplot(x,y,'r');+ b; Y8 Q1 G- D9 b. B% s
zoom on9 \7 d- `1 V# V
p=[1 -4 9 -10];3 I5 |. L! l( @" [3 G
r=roots(p)//次数小于等于4,若大于4则用该命令不一定能求出所有根的准确形式。
7 R+ X7 N. X; L, V% W# Nr=solve('x^3-4*x^2+9*x-10=0')//简单方法
5 H- @3 f9 ~" I/ q& T* d3 H5 ]r=solve('x^2-a*x-4*b=0','x')1 w! r" c0 X; v$ q3 h, x' |
d=solve('p*sin(x)-r','x')5 H6 b; \- Q) i/ q7 U" T
r=fzero('x-(cos(x))^2',1)//其中必须输入的是方程对应的函数名/ z9 N6 |" H4 O: ~5 d9 T
[x,y]=solve('x+y=1','x-11*y=5','x','y')
zan
转播转播0 分享淘帖0 分享分享0 收藏收藏0 支持支持0 反对反对0 微信微信
mrx 实名认证       

31

主题

7

听众

799

积分

升级  49.75%

该用户从未签到

新人进步奖 发帖功臣

群组数学趣味、游戏、IQ等

群组数学建模

群组趣味数学

群组Mathematica研究小组

群组我们的小窝

matlab 里的主成分分析.rar (864 Bytes, 下载次数: 4)

关于校车安排问题.rar

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

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

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

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

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

蚁群算法的基本原理.rar

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

回复

使用道具 举报

3

主题

4

听众

18

积分

升级  13.68%

该用户从未签到

新人进步奖

回复

使用道具 举报

3

主题

4

听众

18

积分

升级  13.68%

该用户从未签到

新人进步奖

回复

使用道具 举报

3

主题

4

听众

18

积分

升级  13.68%

该用户从未签到

新人进步奖

回复

使用道具 举报

xxj198712        

0

主题

4

听众

100

积分

升级  0%

该用户从未签到

回复

使用道具 举报

zhula        

0

主题

3

听众

76

积分

升级  74.74%

该用户从未签到

回复

使用道具 举报

zhula        

0

主题

3

听众

76

积分

升级  74.74%

该用户从未签到

回复

使用道具 举报

zhula        

0

主题

3

听众

76

积分

升级  74.74%

该用户从未签到

回复

使用道具 举报

byit        

4

主题

4

听众

519

积分

升级  73%

该用户从未签到

新人进步奖

回复

使用道具 举报

您需要登录后才可以回帖 登录 | 注册地址

qq
收缩
  • 电话咨询

  • 04714969085
fastpost

关于我们| 联系我们| 诚征英才| 对外合作| 产品服务| QQ

手机版|Archiver| |繁體中文 手机客户端  

蒙公网安备 15010502000194号

Powered by Discuz! X2.5   © 2001-2013 数学建模网-数学中国 ( 蒙ICP备14002410号-3 蒙BBS备-0002号 )     论坛法律顾问:王兆丰

GMT+8, 2026-6-11 16:19 , Processed in 0.380435 second(s), 105 queries .

回顶部