| 查看: 687 | 回复: 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 ] |
回复此楼» 猜你喜欢
材料科学与工程339求调剂
已经有9人回复
-
材料与化工调剂一志愿大连海事085600，349
已经有9人回复
-
292求调剂
已经有8人回复
-
275求调剂
已经有7人回复
-
材料调剂
已经有11人回复
-
288求调剂 一志愿哈工大 材料与化工
已经有25人回复
-
085602化学工程268分蹲调剂
已经有4人回复
-
282求调剂
已经有14人回复
-
301求调剂
已经有3人回复
-
总分343，求生物学调剂
已经有3人回复
» 抢金币啦！回帖就可以得到:
-
2026级药学学硕调剂招生
+1/287
-
陕西理工大学资源与环境、材料化工部分调剂名额，理工兼收
+1/185
-
东莞理工土木工程和土木水利专业硕士接受调剂，名额充足！
+1/87
-
一志愿 中国科学技术大学生物学 365分 求调剂
+1/78
-
下一站，钦州
+1/76
-
烟台大学环境与材料工程学院2026年硕士研究生招生
+1/46
-
考研调剂!
+1/40
-
华东师范大学物理化学专业2026级博士研究生招生（第二批）
+1/38
-
【林业专业硕士调剂1名】中国林业科学研究院
+1/33
-
上海大学 “生物有机电子材料及器件”团队硕士、博士研究生招聘
+1/26
-
南开大学环境学院 环境微生物方向 2026年博士招生（名额1人）
+1/20
-
南开大学环境学院 环境微生物方向 2026年博士招生（名额1人）
+1/19
-
丽水学院，材料与化工专业
+1/9
-
【博士招生】广东工业大学国家优青课题组招收2026年博士生（环境/化学/材料/微生物）
+1/8
-
调剂招生
+1/7
-
上海应用技术大学化工学院招收调剂
+1/6
-
齐鲁工业大学（山东省科学院） 谷亚威老师课题组拟招2026年调剂研究生
+1/5
-
齐齐哈尔大学李莉课题组诚招2026级考研调剂生（学硕和专硕）
+1/3
-
哈尔滨工程大学动力学院赵建辉团队招收2026博士研究生
+1/3
-
2026年江西科技师范大学Tide课题组招收化学、生物、药学、材料相关学硕或专硕调剂生
+1/1
2楼2011-01-04 19:06:32
3楼2011-01-04 19:31:42
4楼2011-01-04 19:47:16
5楼2011-01-04 19:49:24
6楼2011-02-28 21:36:36
