一阶导数 matlab,三阶样条插值(一阶导数边界条件) matlab程序
该楼层疑似违规已被系统折叠 隐藏此楼查看此楼
function
S=Threch1(X,Y,dy0,dyn,xi)
% X
为已知数据的横坐标
%Y
为已知数据的纵坐标
%xi
插值点处的横坐标
%S
求得的三次样条插值函数的值
%dy0
左端点处的一阶导数
% dyn
右端点处的一阶导数
n=length(X)-1;
d=zeros(n+1,1);
h=zeros(1,n-1);
f1=zeros(1,n-1);
f2=zeros(1,n-2);
for
i=1:n
%
求函数的一阶差商
h(i)=X(i+1)-X(i);
f1(i)=(Y(i+1)-Y(i))/h(i);
end
for
i=2:n
%
求函数的二阶差商
f2(i)=(f1(i)-f1(i-1))/(X(i+1)-X(i-1));
d(i)=6*f2(i);
end
d(1)=6*(f1(1)-dy0)/h(1);
d(n+1)=6*(dyn-f1(n-1))/h(n-1);
%¸
赋初值
A=zeros(n+1,n+1);
B=zeros(1,n-1);
C=zeros(1,n-1);
for
i=1:n-1
B(i)=h(i)/(h(i)+h(i+1));
C(i)=1-B(i);
end
A(1,2)=1;
A(n+1,n)=1;
for
i=1:n+1
A(i,i)=2;
end
for
i=2:n
A(i,i-1)=B(i-1);
A(i,i+1)=C(i-1);
end
M=A\d;
syms
x
;
for
i=1:n
Sx(i)=collect(Y(i)+(f1(i)-(M(i)/3+M(i+1)/6)*h(i))*(x-X(i))
...
+M(i)/2*(x-X(i))^2+(M(i+1)-M(i))/(6*h(i))*(x-X(i))^3);
digits(4);
Sx(i)=vpa(Sx(i));%
三样条插值函数表达式
end
for
i=1:n
disp(
'S(x)='
);
fprintf(
'%s (%d,%d)\n'
,char(Sx(i)),X(i),X(i+1));
end
for
i=1:n
if
xi>=X(i)&&xi<=X(i+1)
S=Y(i)+(f1(i)-(M(i)/3+M(i+1)/6)*h(i))*(xi-X(i))+M(i)/2*(xi-X(i))^2+(M
(i+1)-M(i))/(6*h(i))*(xi-X(i))^3;
end
end
disp(
'xi S'
);
fprintf(
'%d,%d\n'
,xi,S);
return
一阶导数 matlab,三阶样条插值(一阶导数边界条件) matlab程序相关推荐
- matlab如何表示一阶导数,三阶样条插值(一阶导数边界条件) matlab程序
该楼层疑似违规已被系统折叠 隐藏此楼查看此楼 function S=Threch1(X,Y,dy0,dyn,xi) % X 为已知数据的横坐标 %Y 为已知数据的纵坐标 %xi 插值点处的横坐标 %S ...
- matlab怎么求三阶导数,三阶样条插值(一阶导数边界条件) matlab程序
该楼层疑似违规已被系统折叠 隐藏此楼查看此楼 function S=Threch1(X,Y,dy0,dyn,xi) % X 为已知数据的横坐标 %Y 为已知数据的纵坐标 %xi 插值点处的横坐标 %S ...
- matlab spline边界约束,三阶样条插值(一阶导数边界条件) matlab程序
该楼层疑似违规已被系统折叠 隐藏此楼查看此楼 function S=Threch1(X,Y,dy0,dyn,xi) % X 为已知数据的横坐标 %Y 为已知数据的纵坐标 %xi 插值点处的横坐标 %S ...
- matlab数组求一阶导数,用matlab求符号函数数组或函数矩阵的导数实例教程
用matlab求符号函数数组或函数矩阵的导数实例教程 在符号矩阵中,矩阵的每个元素都可以存放符号函数,这是对符号函数数组的求到运算是对数组元素逐个进行的.接下来就以实例向大家介绍一下用matlab对符 ...
- 聚类分析matlab检验,「matlab聚类分析」聚类分析的Matlab 程序—系统聚类(附有案例分析) - 金橙教程网...
matlab聚类分析 聚类分析的Matlab 程序-系统聚类 (1)计算数据集每对元素之间的距离,对应函数为pdistw. 调用格式:Y=pdist(X),Y=pdist(X,'metric'), Y ...
- morlet小波matlab代码,morlet小波matlab程序
论-Matlab 7.0复Morlet小波分析操作实例.? 92? 2015年5 月 科技创新 中文科技期刊数据库 ( 文摘版 )自然科学 Matlab 7.0复MorletP]~ 波分析 操作 .. ...
- matlab中存档算法代码,MATLAB 智能算法超级学习手册中程序代码
[实例简介] MATLAB 智能算法超级学习手册中程序代码 [实例截图] [核心代码] dc90ef43-7920-434e-bdb8-0636c31c0b44 └── MATLAB 智能算法超级学习 ...
- MATLAB之易经卜卦程序+GUI
MATLAB之易经卜卦程序+GUI 日月为易,刚柔相推. 是故易有太极,是生两仪,两仪生四象,四象生八卦,八卦定吉凶,吉凶生大业.是故法象莫大乎天地,变通莫大乎四时,悬象著明莫大乎日月. 本文 ...
- 光伏并网逆变器设计方案,附有相关的matlab电路仿真文件,以及DSP的程序代码
光伏并网逆变器设计方案,附有相关的matlab电路仿真文件,以及DSP的程序代码,方案.仿真文件.代码三者结合使用效果好,事半功倍 ID:6910671781885764
- matlab怎么求三次微分,matlab课设三阶微分方程多种方法求解.doc
matlab课设三阶微分方程多种方法求解 目录 一.课程设计题目及意义 -------- 1 页 二.课程设计任务及要求 --------2 页 三.课程设计详细过程及结果 --------3至10页 ...
最新文章
- eclipsse 关闭 方法提示_MacOS10.15.7关闭AppStore右上角小红标提示及系统更新右上角小红标提示的方法...
- 78. Subsets 子集
- Java-Concurrent 线程池ThreadPoolExecutor使用
- oracle估算数据增长,如何估算oracle 数据库,数据库对象历史增长情况
- tp6 实现汉字转拼音
- nutch ajax mysql_Nutch爬取Ajax请求的动态网页
- 微生物组-扩增子16S分析和可视化(2023.2)
- PB调用DLL的常见问题及处理方法
- 读书笔记——OpenGL超级宝典
- 二手房中介管理系统设计
- 【数据可视化应用】Python反距离权重(IDW)插值计算及可视化绘制
- canvas rotate() 中心旋转的实际运用
- 粉碎机无法粉碎文件_粉碎站点地图
- any_value()
- IDEA启动main函数时,Command line is too long Shorten command line for XXX or also问题解决
- 使用ROT13加密解密
- python小问题:依赖包/conda-forge/conda无法安装包/librosa/
- 让网站更新快照的方法
- 基于Java毕业设计新锐台球厅管理系统源码+系统+mysql+lw文档+部署软件
- 计算机定时关机教程,电脑定时关机怎么设置
热门文章
- App和小程序开发成本对比
- WiFi速率控制:吟游诗人速率控制算法(minstrel rate control algorithm for 802.11)
- Cannot format given Object as a Number
- ps cc2019 安装教程
- MacBook 更新Big Sur后,虚拟机无法运行时显示该主机 CPU 类型不支持虚拟化性能计数器,模块“VPMC”启动失败,未能启动虚拟机
- 随机中文姓名 php,PHP生成随机中文姓名 阿星小栈
- win7无法自动检测修复计算机,无法启动,启动修复正在检查系统问题
- 1138: C语言合法标识符
- squid高性能代理缓存服务器
- 南京软件测试女薪资2019,2019年第一扎!南京最新平均工资最新出炉,看完不敢过年!...