应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » MATLAB/Mathematica » 对数据向量求导
41/1返回列表
查看: 2502  |  回复: 3
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

小蛮漫漫

新虫 (初入文坛)

[求助] 对数据向量求导

m = load('123.txt');
x = m(:,1);
y = m(:,2);
                               %%     123.txt中有两列数据，左边一列从30-30000，步长是30，总共1000个数据；
                              %%     右边是对应的x取不同值的时候对应的值。
怎么求 y 在x = 3000处的导数啊，求高手指教！！！！
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2013-06-28 18:43:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ajitai

铁杆木虫 (职业作家)
你必须先拟合出函数来，才能够求导吧。没有函数，只有数据可能不行吧。
赞一下 回复此楼
2楼2013-06-30 08:52:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sukiyq

木虫 (小有名气)

【答案】应助回帖

感谢参与，应助指数 +1
数值差分吧,比如：

一阶中心差分 f'(x) = [f(x+s) - f(x-s)]/2/s
一阶向前差分 f'(x) = f(x+s) - f(x)
赞一下 回复此楼
比菜鸟强一点点
3楼2013-07-01 11:24:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sukiyq

木虫 (小有名气)

【答案】应助回帖

数值差分吧,比如：

一阶中心差分 f'(x) = [f(x+s) - f(x-s)]/2/s
一阶向前差分 f'(x) = [f(x+s) - f(x)]/s

用泰勒展式还可以得到其他形式的精度更高的导数表达式，比如

f'(x) = (8f(x+s) - 8f(x-s) - f(x+2s) +f(x-2s))/12s

其中s是步长
赞一下 回复此楼
比菜鸟强一点点
4楼2013-07-01 11:34:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 小蛮漫漫 的主题更新
41/1返回列表
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭