| 查看: 774 | 回复: 5 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[交流]
【求助】急求帮忙解决MATLAB 样条差值求解微分代码出现的问题，使代码正常运行
|
|||
|
function y0=hillsplines3(x0,x,y) if size(x,1)~=1 x=x'; end if size(y,1)~=1 y=y'; end [a,t]=hillsplinex(x,y);a=a';t=t'; for i=1:length(t)-4 phi(i)=powerplus3(x0,t(i:i+4)); end y0=sum(a.*phi); function [ax,tx]=hillsplinex(t,y) h=diff(t);delta=diff(y)./h;n=length(h)+1 d1=pchipendpoint(h(1),h(2),delta(1),delta(2)); dn=pchipendpoint(h(n-1),h(n-2),delta(n-1),delta(n-2)); r=[d1;y';dn]; tx=[t(1)-[1:3]*h(1),t,t(length(t)+[1:3])*h(length(h))]; a=zeros(length(tx)-4,length(tx)-4); for i=1:3 a(1,i)=depowerplus3(t(1),tx(i:i+4)); end for i=length(tx)-6:length(tx)-4 a(length(tx)-4,i)=depowerplus3(t(length(t)),tx(i:i+4)); end for j=2:length(tx)-5 for i=j-1:j+1 a(j,i)=powerplus3(t(j-1),tx(i:i+4)); end end ax=a\r;tx=tx'; function y=powerplus3(x,t) c=t; for i=1:5 c(i)=[]; beta(i)=24/prod(t(i)-c); c=t; end powerplus=abs(x-t).^3; y=0.5*sum(beta.*powerplus); function y=depowerplus3(x,t) c=t; for i=1:5 c(i)=[]; beta(i)=24/prod(t(i)-c); c=t; end powerplus=3*sign(x-t).*(x-t).^2; y=0.5*sum(beta.*powerplus); function d=pchipendpoint(h1,h2,del1,del2) d=((2*h1+h2)*del1-h1*del2)/(h1+h2); if sign(d)~=sign(del1) d=0; elseif(sign(del1)~=sign(del2))&(abs(d)>abs(3*del1)) d=3*del1; end [ Last edited by huameitang05 on 2011-1-4 at 17:21 ] |
回复此楼» 猜你喜欢
分享一下我之前已中青C的计划书的filecode
已经有5人回复
-
不应该看fileCode
已经有8人回复
-
FileCode能看出啥？
已经有26人回复
-
关于Filecode分析方法
已经有11人回复
-
重要来源：本周末出结果
已经有9人回复
-
我的国基提前知道中了，可是同事的操作让我实在接受不了，怎么会有这样的人
已经有14人回复
-
结合人工智能，周易传统文化，filecode打分制来了，3分以上希望很大。
已经有4人回复
-
2019年青年基金涵评意见，大家看看几个A，几个B?
已经有11人回复
-
应该是93bebmhtak前后十一个字符比较关键
已经有36人回复
-
好奇怪的filecode
已经有6人回复
» 抢金币啦！回帖就可以得到:
-
北京理工大学集成电路学院-柔性电子器件与智造研究所-2027级博士招生
+1/468
-
中国科学院空天院精密干涉测量团队招聘光学/算法/机械方向博士后/特别研究助理
+1/82
-
中国科学技术大学-苏州国家实验室联培招生-本科直博/硕博连读/申请考核三种方式
+1/79
-
物理学报录用，散币感谢合作者和审稿人
+2/40
-
天津大学环境学院童银栋教授招收博士后（生源要素循环/温室效应/新污染物方向）
+1/30
-
哈工大（深圳）何自开团队诚招2027级化学工程与技术专业博士生
+1/26
-
韩国釜山大学储能电池与功能材料实验室诚招博后（电池/电化学/器件制备/MOF合成）
+1/14
-
南京航空航天大学招收27年力学硕士和博士生研究生
+1/11
-
硕士毕业9年的我，读博的想法日益强烈。。。。
+1/9
-
深圳理工大学招收生物医用高分子材料方向客座研究生
+1/6
-
北京理工大学-集成电路与电子学院杰青团队-招博士后
+1/6
-
中国科学院地理资源所陆地表层格局与模拟院重点实验室招聘特别研究助理（博士后）
+1/5
-
中山大学李国强课题组招收博士生（柔性电子方向）
+1/4
-
海南大学化学院-科研助理招聘 （可读博，可与华中科技大学，天津大学联培）+博后招聘
+1/4
-
北理工集成电路杰青团队 | 诚招科助理
+1/4
-
UTS 招收2027年入学博士生
+1/4
-
中山大学计算数学方向招2027届考核制博士生
+1/3
-
国家级人才团队诚聘生物信息相关专业的科研助理
+1/3
-
苏州大学能源学院管大秦课题组诚聘博士后及招收博士/硕士研究生
+1/2
-
说同济非升即走那是毛毛雨，青岛理工大学（原青岛建工学院）那才叫变态非升即走呢！
+1/1
6楼2011-02-28 21:36:36
2楼2011-01-04 19:06:32
3楼2011-01-04 19:31:42
4楼2011-01-04 19:47:16
