找回密码
 注册
关于网站域名变更的通知
查看: 409|回复: 2
打印 上一主题 下一主题

请教一个程序

[复制链接]

该用户从未签到

跳转到指定楼层
1#
发表于 2020-1-9 10:13 | 只看该作者 回帖奖励 |倒序浏览 |阅读模式

EDA365欢迎您登录!

您需要 登录 才可以下载或查看,没有帐号?注册

x
h = 0.001;       % step sizeX = -pi:h:pi;    % domainf = sin(X);      % rangeY = diff(f)/h;   % first derivativeZ = diff(Y)/h;   % second derivativeplot(X(:,1:length(Y)),Y,'r',X,f,'b', X(:,1:length(Z)),Z,'k')请问plot中的(:,1:length(Y))是什么含义?对于:和length(Y)含义不是很懂* @/ x8 y" `2 Q2 v6 g; {* t
* C0 f4 j5 s7 G' D; A! |  t. U

4 y& e2 {5 `* z

该用户从未签到

2#
发表于 2020-1-9 18:05 | 只看该作者
是X(:,1:length(Y)),在本例中,X是一个行向量,:表明取X向量的第一行,length(Y)计算Y向量的长度,1:length(Y)表明X向量从第一个值开始取长度为length(Y)的向量。在这里,使用一次diff函数会导致计算结果向量比X向量少一个数,使用X(:,1:length(Y))是为了确保绘图的时候横纵坐标的数据长度一致。

该用户从未签到

3#
发表于 2020-1-9 21:40 | 只看该作者
先说明length(y):0 j+ x) k: D+ y2 W* c: P. X( Y
1、若y为行向量/列向量,则返回该向量的长度
" m' A; A9 \# A5 {7 P8 x6 b5 D2、若y为m*n的矩阵,则返回n与m二者中的最大值
1 ~6 r$ l7 }, l2 ~0 [原程序对应情况一
9 b+ }1 d  Z( O5 A" q. F  k9 u; M
接着‘:’符号:- t* u  W& l9 E4 A+ E; }
例:代码:1:6
- v$ E( y( r6 s8 p  y9 M6 y       实际:[1 2 3 4 5 6]    %默认间隔为1的1至6的向量
5 ~3 F7 E+ l7 ^" s( W1 o       代码:1:0.5:4
; M8 ?7 _" j& w1 s4 J        实际:[1 1.5 2 2.5 3 3.5  4 ]    %间隔为0.5的1至4的向量
. N; C7 Y6 o# Q
; O8 e* }$ }4 D2 n4 q! N9 vmatlab中diff函数使用相邻两位的数据求解近似导数,所以若x有n个数值,y就只能有n-1个$ y) h4 x% b. S: C4 m; u
而plot(a,b)这个函数必须保证输入向量a与b长度相等% ~& g/ u) I; [  D, l8 E, l
为保证输出,需要用y的长度来限制x取用的长度) s- l4 z( \7 v1 h# v7 l
x(1:length(y))是取x中第1到n-1个元素作为向量
您需要登录后才可以回帖 登录 | 注册

本版积分规则

关闭

推荐内容上一条 /1 下一条

EDA365公众号

关于我们|手机版|EDA365电子论坛网 ( 粤ICP备18020198号-1 )

GMT+8, 2025-11-1 20:27 , Processed in 0.156250 second(s), 24 queries , Gzip On.

深圳市墨知创新科技有限公司

地址:深圳市南山区科技生态园2栋A座805 电话:19926409050

快速回复 返回顶部 返回列表