QQ登录

只需要一步,快速开始

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

比较通过符号计算得到的导数和通过数值差分方法计算得到的导数的结果

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

1186

主题

4

听众

2922

积分

该用户从未签到

跳转到指定楼层
1#
发表于 2024-9-26 17:21 |只看该作者 |倒序浏览
|招呼Ta 关注Ta
h=0.05; x=0:h:pi; syms x1; y=sin(x1)/(x1^2+4*x1+3);
3 Z0 [$ b- B2 ]  a- i4 l8 F, z9 e* zyy1=diff(y); f1=subs(yy1,x1,x);   % 求各阶导数的解析解与对照数据
- f/ Z$ X8 |; l' _) P4 ]yy2=diff(yy1); f2=subs(yy2,x1,x); yy3=diff(yy2); f3=subs(yy3,x1,x);) `) z- ]- ]; ^' P7 S
yy4=diff(yy3); f4=subs(yy4,x1,x);! x7 N: ?' d& A* W
& l4 W- ]1 r" ^7 j0 l$ t$ i
y=sin(x)./(x.^2+4*x+3);   % 生成已知数据点7 {+ e( H* i; o
[y1,dx1]=diff_ctr(y,h,1); subplot(221),plot(x,f1,dx1,y1,':');, \) a0 G5 H6 h% t
[y2,dx2]=diff_ctr(y,h,2); subplot(222),plot(x,f2,dx2,y2,':')- i3 I& Y+ j7 w1 b/ T; a. _6 T$ p
[y3,dx3]=diff_ctr(y,h,3); subplot(223),plot(x,f3,dx3,y3,':');# u/ ~* S; r5 X$ j1 T/ z4 F
[y4,dx4]=diff_ctr(y,h,4); subplot(224),plot(x,f4,dx4,y4,':')3 [/ ~! G: K2 [. S% P! M5 c8 y
& f+ q- H4 h8 \; B# D9 `
norm((y4-f4(4:60))./f4(4:60))
, `; @# L) y  T. U这段MATLAB代码旨在比较通过符号计算得到的导数和通过数值差分方法计算得到的导数的结果。以下是对代码的逐步解析:
+ X# i% r0 H( g0 N' q" e7 L
& S. B- ]4 L7 ?; |

$ Q- p6 A' N1 ^5 I9 h# k### 1. 定义和设置/ |6 Y! |& a( ^
```matlab4 L2 I* c- r3 ?5 s$ @& F
h = 0.05; 8 ]7 [& r1 T+ F+ P% C( B, l
x = 0:h:pi;
6 a0 b1 D' O1 M7 ^' b* e3 K8 rsyms x1; 0 _+ H& l) c* r4 I0 E# N& |, S. f
y = sin(x1) / (x1^2 + 4*x1 + 3);
' ]& }3 b( W6 o- P5 ^0 {```  ~2 H' |4 Q, u# F
- `h = 0.05;` 定义了一个小步长 `h`,用于生成数据点。
+ z, i+ A2 t, r6 u8 d9 i9 h- `x = 0:h:pi;` 生成从0到π的等间距点。
3 L5 F6 I' _4 r6 I( S: ?; r# l- `syms x1;` 定义了一个符号变量 `x1`。2 v! _& u; l; S3 I6 q; B7 _
- `y = sin(x1) / (x1^2 + 4*x1 + 3);` 定义了符号函数 `y`,这是一个关于 `x1` 的函数。9 P$ r4 N7 u! n2 j+ h

9 j" Q4 {, i* Q  N* i. c
- L' @# u; W6 |: ?. [9 b1 e; B
### 2. 计算导数
; Z2 H8 y* Y3 Q& s( B```matlab7 k% T  N+ [0 ?; B
yy1 = diff(y);
- p% ^  g; N' n* H1 zf1 = subs(yy1, x1, x);' w5 J" _# s$ t5 G
```9 ^  J' @, s4 k
- `yy1 = diff(y);` 计算函数 `y` 的第一阶导数。
: D; o1 @: ]' x- `f1 = subs(yy1, x1, x);` 使用 `subs` 将导数中的符号变量 `x1` 替换为具体的 `x` 值,得到导数的值,储存在 `f1` 中。  _$ m4 k8 u1 f8 n; E

, G+ z) b* g/ k" i9 [$ U
6 b0 `% r5 a+ ~" s* J" l0 h) f
对后续导数的计算重复相同的操作:
5 |) C$ }* d% f6 J0 X9 n' m6 v```matlab
3 C( W) D: k* z9 I8 byy2 = diff(yy1);
* V6 R- ^3 T$ A0 ]; Xf2 = subs(yy2, x1, x);
7 h* m% [" R. v+ a4 k) Eyy3 = diff(yy2);
% }1 }# R+ }) @4 ]7 Vf3 = subs(yy3, x1, x);; h7 [  d' ]7 t; X) }% l
yy4 = diff(yy3); ! S# T& j: W/ \
f4 = subs(yy4, x1, x);
* P, x. m) n4 m6 V```
* Y" d) n8 ~% }3 E* s- 上述代码计算了 `y` 的第二、第三和第四阶导数,并相应地进行了替换得到了 `f2`、`f3` 和 `f4`。8 \) E6 q7 Z2 K( e8 I9 [

/ P9 ?$ o- Z* p4 _
1 Y3 T% k2 ?/ M4 Q
### 3. 生成已知数据点
2 @' E( R# d" v```matlab
, a. h7 m: O0 P% yy = sin(x) ./ (x.^2 + 4*x + 3); , V. ?% A8 [6 }
```
& _1 t2 n/ m6 J- 计算 `y` 的值,将 `x` 中的每个元素代入,生成已知数据点。
: M4 r! U8 A. Y# j
3 {1 u3 j0 {* ^( ]# n, F5 \/ g
# \5 o* c& m0 C6 H5 |$ ]
### 4. 使用差分方法计算导数& V- c+ r0 ~$ G
以下是使用数值方法计算导数的步骤,假设 `diff_ctr` 函数已经定义并实现:
. ^/ t5 Y- ?+ v( C/ E1 ?
6 ^' Q! S% \3 W/ v' [5 _! X% S! {* V
( i% v. ?5 {: e* _  X4 c2 ~
```matlab+ ?8 H# D2 p6 Y  g' v% H, t/ L
[y1, dx1] = diff_ctr(y, h, 1); " m! i$ J: S4 O$ `3 s6 b$ y
subplot(221), plot(x, f1, dx1, y1, ':');' D6 e1 L+ T+ W& k/ j
[y2, dx2] = diff_ctr(y, h, 2); ( Y+ A4 W" I8 j% w* _* o  `: y
subplot(222), plot(x, f2, dx2, y2, ':');* v4 f* S& L1 B9 t# Q9 g0 H  |
[y3, dx3] = diff_ctr(y, h, 3); ; s. _9 |) }7 b" h( s% u
subplot(223), plot(x, f3, dx3, y3, ':');- p' j1 }" E+ r2 U  M" a  V- M
[y4, dx4] = diff_ctr(y, h, 4);
/ L+ S& p, J8 I  ]subplot(224), plot(x, f4, dx4, y4, ':');7 M# s5 x) n2 |( W. H. ]
```
: Z8 m' a7 c4 W2 G! x% O4 p- `diff_ctr(y, h, n)` 函数可能是一个自定义的函数,用来通过中心差分方法计算函数 `y` 的 `n` 阶导数。
2 i9 h  s  t5 W' K- `subplot` 和 `plot` 用于可视化结果,将每个计算的导数和数值差分结果绘制在不同的子图中。
6 g% `2 [+ v2 g' A* _
, G+ e; R+ ?' k9 S. u1 j4 O( K. X
" [4 c" O( {+ j( O3 q
### 5. 计算误差的规范化) S: I" C! y& ^
```matlab/ v  ?( s% _4 ~" N& ~0 _
norm((y4 - f4(4:60)) ./ f4(4:60))
& X5 D% c% E( @/ g# [' w```' t* C! l7 g" ]9 U% Y1 g
- `norm(...)` 用于计算向量的范数,这里比较后处理的数值导数 `y4` 和符号计算的导数 `f4(4:60)`,以计算相对误差。
; |" a" ^* F7 v& W, }- 通过这种方法,用户可以评估数值计算的精度。
6 m* [$ p6 g& E$ b' x, w8 y4 l# f7 f

) y6 y4 i$ q1 J& N### 总结
' V. Y$ a3 B7 s3 f这段代码的主要目标是通过符号计算和中心差分法计算函数的导数,并绘制出它们的比较图。最终的规范化误差提供了一种衡量数值方法精度的手段。这种方法对理解和验证数值微分的重要性以及与解析解的比较相当有用。
' a9 T5 O7 g, _+ G
1 P/ {9 k, Z: _- f8 A$ b7 }6 {# T2 N. b7 G

3 s. q8 L7 n: N' `& w
7 L$ H) w. X' I$ S/ V4 |% g7 [, n& J# d( {4 X& f0 ^7 V9 E) M1 Q
zan
转播转播0 分享淘帖0 分享分享0 收藏收藏0 支持支持0 反对反对0 微信微信
您需要登录后才可以回帖 登录 | 注册地址

qq
收缩
  • 电话咨询

  • 04714969085
fastpost

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

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

蒙公网安备 15010502000194号

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

GMT+8, 2026-4-10 16:44 , Processed in 1.495243 second(s), 51 queries .

回顶部