qulv 利用MATLAB处理图形图像

利用MATLAB处理图形图像，求曲线的曲率

%求曲率
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