三次样条插值及三弯矩法完整（Matlab实现），2024年最新阿里程序员面试

% dy1=5; % 边界条件代入

% dy9=29.16; % 边界条件代入

% d(1)=6/h(1)*((y(2)-y(1))/h(1)-dy1); % 边界条件代入

% d(9)=6/h(8)*(dy9-(y(9)-y(8))/h(8)); % 边界条件代入

% u(i-1)=h(i-1)/(h(i-1)+h(i)); % 公式

% v(i)=1-u(i-1); % 公式

% d(i)=6*fejcs(i-1); % 公式

%end

%%三弯矩方程组：AM=d

%A=2*eye(9,9)+diag(v,1)+diag(u,-1)

%d1=d’

%AN=inv(A);

%M=AN*d1

%M2=M(2,1);

%M3=M(3,1);

%M7=M(7,1);

%M8=M(8,1);

%syms z c

%S2(z)=((((x(3)-z)^{3)/(6*h(3)))*M2)+(((z-x(2))}3)/(6h(3))M3)+((y(2)-(M2(((h(3))^{2)/6)))*(x(3)-z)/h(3))+((y(3)-(M3*((h(3))}2/6)))(z-x(2))/h(3));

%z=-0.02;

%fzl=vpa(S2(z))

%S7©=((((x(8)-c)^{3)/(6*h(8)))*M7)+(((c-x(7))}3)/(6h(8))M8)+((y(7)-(M7(((h(8))^{2)/6)))*(x(8)-c)/h(8))+((y(8)-(M8*((h(8))}2/6)))(c-x(7))/h(8));

%c=2.56;

%fcl=vpa(S7©)

%%方法二

df1=5;df2=29.16;%边界条件

n=length(x);%计算维度

h=zeros(1,n);%初始化赋值

dy=zeros(1,n);%初始化赋值

u=zeros(1,n);%初始化赋值

lamda=zeros(1,n);%初始化赋值

d=zeros(1,n);%初始化赋值

l=zeros(1,n);%初始化赋值

r=zeros(1,n);%初始化赋值

z=zeros(1,n);%初始化赋值

M=zeros(1,n);%初始化赋值

for i=2:n

h(i)=x(i)-x(i-1); %步长

end

for i=1:n-1

dy(i)=(y(i)-y(i+1))/(x(i)-x(i+1)); %一阶差商

end

for i=2:n-1

u(i)=h(i)/(h(i)+h(i+1));%求u

lamda(i)=1-u(i);%求lamda

d(i)=6*(dy(i)-dy(i-1))/(h(i)+h(i+1));%求d

end

d(1)=6*(dy(1)-df1)/h(2);

d(n)=6*(df2-dy(8))/h(n);

u(n)=1; %u的最后一个元素为1

lamda(1)=1; %lamda第1个元素为1

r(1)=2;

for i=2:n %追赶法求解M

l(i)=u(i)/r(i-1);

r(i)=2-l(i)*lamda(i-1);

end

z(1)=d(1);

for i=2:n

z(i)=d(i)-l(i)*z(i-1);

end

M(n)=z(n)/r(n);

for i=n-1👎1

M(i)=(z(i)-lamda(i)*M(i+1))/r(i)

end

c=input(‘请输入c：’);

for k=2:8

if x(k-1)<=c & c<=x(k)

S=M(k-1)/6/h(k)(x(k)-c)^{3+M(k)/6/h(k)*(c-x(k-1))}3+1/h(k)(y(k)-M(k)h(k)^{2/6)*(c-x(k-1))+1/h(k)*(y(k-1)-M(k-1)*h(k)}2/6)(x(k)-c)%第一公式

%S=((((x(k)-c)^{3)/(6*h(k)))*M(k-1))+(((c-x(k-1))}3)/(6h(k))M(k))+((y(k)-(M(k-1)(((h(k))^{2)/6)))*(x(k)-c)/h(k))+((y(k)-(M(k)*((h(k))}2/6)))(c-x(k-1))/h(k))%第二公式

return;

end

自我介绍一下，小编13年上海交大毕业，曾经在小公司待过，也去过华为、OPPO等大厂，18年进入阿里一直到现在。

深知大多数Python工程师，想要提升技能，往往是自己摸索成长或者是报班学习，但对于培训机构动则几千的学费，着实压力不小。自己不成体系的自学效果低效又漫长，而且极易碰到天花板技术停滞不前！

因此收集整理了一份《2024年Python开发全套学习资料》，初衷也很简单，就是希望能够帮助到想自学提升又不知道该从何学起的朋友，同时减轻大家的负担。

既有适合小白学习的零基础资料，也有适合3年以上经验的小伙伴深入学习提升的进阶课程，基本涵盖了95%以上Python开发知识点，真正体系化！

由于文件比较大，这里只是将部分目录大纲截图出来，每个节点里面都包含大厂面经、学习笔记、源码讲义、实战项目、讲解视频，并且后续会持续更新

如果你觉得这些内容对你有帮助，可以添加V获取：vip1024c （备注Python）

最后

🍅 硬核资料：关注即可领取PPT模板、简历模板、行业经典书籍PDF。
🍅 技术互助：技术群大佬指点迷津，你的问题可能不是问题，求资源在群里喊一声。
🍅 面试题库：由技术群里的小伙伴们共同投稿，热乎的大厂面试真题，持续更新中。
🍅 知识体系：含编程语言、算法、大数据生态圈组件（Mysql、Hive、Spark、Flink）、数据仓库、Python、前端等等。

一个人可以走的很快，但一群人才能走的更远。不论你是正从事IT行业的老鸟或是对IT行业感兴趣的新人，都欢迎扫码加入我们的的圈子（技术交流、学习资源、职场吐槽、大厂内推、面试辅导），让我们一起学习成长！

ark、Flink）、数据仓库、Python、前端等等。

一个人可以走的很快，但一群人才能走的更远。不论你是正从事IT行业的老鸟或是对IT行业感兴趣的新人，都欢迎扫码加入我们的的圈子（技术交流、学习资源、职场吐槽、大厂内推、面试辅导），让我们一起学习成长！
[外链图片转存中…(img-OCnBxkzH-1712774919476)]