版塊導(dǎo)航: 正在加載中...

登錄注冊

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進行實名認證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號能夠正常使用，請盡快對帳號進行手機號驗證，感謝您的理解與支持！

24小時熱門版塊排行榜

返回列表

當前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖

huameitang05

銀蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 250.7
帖子: 47
在線: 59.6小時
蟲號: 99774

[交流] 【求助】急求幫忙解決MATLAB 樣條差值求解微分代碼出現(xiàn)的問題，使代碼正常運行

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 ]

回復(fù)此樓

» 猜你喜歡

333求調(diào)劑已經(jīng)有10人回復(fù)
307求調(diào)劑已經(jīng)有3人回復(fù)
考研化學(xué)308分求調(diào)劑已經(jīng)有4人回復(fù)
298求調(diào)劑已經(jīng)有7人回復(fù)
263求調(diào)劑已經(jīng)有5人回復(fù)
276求調(diào)劑。有半年電池和半年高分子實習(xí)經(jīng)歷已經(jīng)有6人回復(fù)
求調(diào)劑材料學(xué)碩080500，總分289分已經(jīng)有21人回復(fù)
北科281學(xué)碩材料求調(diào)劑已經(jīng)有8人回復(fù)
08工學(xué)調(diào)劑已經(jīng)有9人回復(fù)
291 求調(diào)劑已經(jīng)有4人回復(fù)

» 搶金幣啦！回帖就可以得到:

查看全部散金貼

1樓 2011-01-04 17:13:22

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

honeypirl

銅蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 46.8
帖子: 5
在線: 14.3小時
蟲號: 1185008

★
huameitang05(金幣+1):謝謝參與

hillsplinex是個函數(shù)嗎？

贊一下(1人)

回復(fù)此樓

6樓2011-02-28 21:36:36

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

查看全部 6 個回答

lijinfeng042

木蟲 (小有名氣)

★ ★ ★ ★ ★
huameitang05(金幣+1):謝謝參與
huameitang05(金幣+1):謝謝意見 2011-01-04 19:14:27
xiegangmai(金幣+2):辛苦了！ 2011-01-07 22:06:47
robert2020(金幣+2):辛苦了！最近比較忙，評帖不及時，敬請見諒！ 2011-01-11 10:47:08

引用回帖:

Originally posted by huameitang05 at 2011-01-04 17:13:22:
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));
en ...

具體太多代碼沒看...只是說說思路離散數(shù)據(jù)的微分問題
既然你用的是樣條插值那就說說這個方法
數(shù)據(jù)->csapi(x,y)->fnder(cs)->fnval(pp,x)

贊一下(3人)

回復(fù)此樓

2樓2011-01-04 19:06:32

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

shl1025

金蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 1337.3
帖子: 26
在線: 3.6小時
蟲號: 1143568

★
huameitang05(金幣+1):謝謝參與

祝福！程序就是難�。�

回復(fù)此樓

3樓2011-01-04 19:31:42

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

ysh2225

木蟲 (著名寫手)

★
huameitang05(金幣+1):謝謝參與
robert2020(金幣-2):為了他人的方便，請勿在求助帖中純表無意義回復(fù)！ 2011-01-11 10:47:46

版主扣我BB，這大過年的。。55555

[ Last edited by ysh2225 on 2011-1-11 at 14:46 ]

贊一下(2人)

回復(fù)此樓

4樓2011-01-04 19:47:16

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

查看全部 6 個回答

普通表情龍兔虎貓高級回復(fù) (可上傳附件)

最具人氣熱帖推薦 [查看全部]		作者	回/看	最后發(fā)表

[考研] 293求調(diào)劑 +3	濤濤Wjt 2026-03-22	5/250	2026-03-22 22:21 by jiangpengfei
[考研] 0854電子信息求調(diào)劑 +3	α____ 2026-03-22	3/150	2026-03-22 21:28 by zhq0425
[考研] 石河子大學(xué)（211、雙一流）碩博研究生長期招生公告 +3	李子目 2026-03-22	3/150	2026-03-22 21:01 by 怎么釋懷
[考研] 求調(diào)劑 +7	Auroracx 2026-03-22	7/350	2026-03-22 12:38 by 素顏傾城1988
[考研] 一志愿深大，0703化學(xué)，總分302，求調(diào)劑 +4	七月-七七 2026-03-21	4/200	2026-03-21 18:20 by 學(xué)員8dgXkO
[考研] 0703化學(xué)297求調(diào)劑 +3	Daisy☆ 2026-03-20	3/150	2026-03-21 17:45 by ColorlessPI
[考研] 機械專碩299求調(diào)劑至材料 +3	kkcoco25 2026-03-16	4/200	2026-03-21 03:52 by JourneyLucky
[考研] 303求調(diào)劑 +5	睿08 2026-03-17	7/350	2026-03-21 03:11 by JourneyLucky
[考研] 材料 336 求調(diào)劑 +3	An@. 2026-03-18	4/200	2026-03-21 01:39 by JourneyLucky
[考研] 304求調(diào)劑 +6	曼殊2266 2026-03-18	6/300	2026-03-21 00:32 by JourneyLucky
[考研] 295求調(diào)劑 +4	一志愿京區(qū)211 2026-03-18	6/300	2026-03-20 23:41 by JourneyLucky
[考研] 330求調(diào)劑 +4	小材化本科 2026-03-18	4/200	2026-03-20 23:13 by JourneyLucky
[考研] 求調(diào)劑，一志愿:南京航空航天大學(xué)大學(xué) ，080500材料科學(xué)與工程學(xué)碩，總分289分 +4	@taotao 2026-03-19	4/200	2026-03-20 22:14 by JourneyLucky
[考研] 一志愿西南交通專碩材料355 本科雙非求調(diào)劑 +5	西南交通專材355 2026-03-19	5/250	2026-03-20 21:10 by JourneyLucky
[考研] 298-一志愿中國農(nóng)業(yè)大學(xué)-求調(diào)劑 +9	手機用戶 2026-03-17	9/450	2026-03-20 14:24 by 無懈可擊111
[考研] 288求調(diào)劑，一志愿華南理工大學(xué)071005 +5	ioodiiij 2026-03-17	5/250	2026-03-19 18:22 by zcl123
[考研] 277調(diào)劑 +5	自由煎餅果子 2026-03-16	6/300	2026-03-17 19:26 by 李leezz
[考研] 有沒有道鐵/土木的想調(diào)劑南林，給自己招師弟中～ +3	TqlXswl 2026-03-16	7/350	2026-03-17 15:23 by TqlXswl
[考研] 283求調(diào)劑 +3	聽風就是雨； 2026-03-16	3/150	2026-03-17 07:41 by 熱情沙漠
[考研] 333求調(diào)劑 +3	文思客 2026-03-16	7/350	2026-03-16 18:21 by 文思客

亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产 高清 中文字幕,99re热久久亚洲综合精品成人,熟妇 一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频

24小時熱門版塊排行榜

huameitang05

[交流] 【求助】急求幫忙解決MATLAB 樣條差值求解微分代碼出現(xiàn)的問題，使代碼正常運行

» 猜你喜歡

» 搶金幣啦！回帖就可以得到: 查看全部散金貼

honeypirl

lijinfeng042

shl1025

ysh2225

亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产高清中文字幕,99re热久久亚洲综合精品成人,熟妇一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频

[交流] 【求助】急求幫忙解決MATLAB 樣條差值求解微分代碼出現(xiàn)的問題，使代碼正常運行

» 搶金幣啦！回帖就可以得到:

查看全部散金貼