我只想知道2的函数的正微分和中心微分的矩阵构成代码
我需要创建矩阵不使用循环
1次观看(过去30天)
显示旧的注释
题目:数值微分
使用Matlab
下表为时间t及其对应的x位置
使用正微分和中心微分格式,通过构造一阶导数的微分矩阵,并将矩阵和向量相乘(不允许使用循环)来计算速度和加速度,因为它是基于矩阵的计算。
正微分格式= x(i+1)-x(i) / t(i+1)-t(i)
***期望使用稀疏***然而也开放用于其他功能
请注意,结束位置与开始位置相同,因此t1和t10的速度应该相等
(请为此提供Matlab代码)
时间t =
0
0.0775
0.2346
0.
3015
0.3682
0.4345
0.5006
0.5665
0.6321
0.6974
X位置
3.1909
3.1818 3.1547 3.1104 -0.0499 -0.0499 3.1104 3.1547 3.1818 3.1909