%求曲率
clc;?clear?all;?close?all;
x0?=?linspace（0?1）;
y0?=?sin（x0）.*cos（x0）;
h?=?abs（diff（[x0（2）?x0（1）]））;
%?模擬一階導
figure;?box?on;?hold?on;
ythe1?=?cos（x0）.^2?-?sin（x0）.^2;?%理論一階導
yapp1?=?gradient（y0?h）;?%matlab數值近似
plot（x0?ythe1?‘.‘）;
plot（x0?yapp1?‘r‘）;
legend（‘理論值‘?‘模擬值‘）;
title（‘模擬一階導‘）;
%?模擬二階導
figure;?box?on;?hold?on;
ythe2?=?（-4）*cos（x0）.*sin（x0）;?%理論二階導
yapp2?=?2*2*del2（

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件?????????873??2012-09-24?08:36??qulv.m