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

登錄注冊

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

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ù)此樓

» 猜你喜歡

本人考085602 化學(xué)工程專碩已經(jīng)有23人回復(fù)
求調(diào)劑院校信息已經(jīng)有4人回復(fù)
085600材料與化工306 已經(jīng)有4人回復(fù)
286求調(diào)劑已經(jīng)有10人回復(fù)
328求調(diào)劑，英語六級551，有科研經(jīng)歷已經(jīng)有9人回復(fù)
一志愿北京化工大學(xué)070300 學(xué)碩336求調(diào)劑已經(jīng)有4人回復(fù)
286分人工智能專業(yè)請求調(diào)劑愿意跨考！已經(jīng)有8人回復(fù)
資源與環(huán)境調(diào)劑申請(333分) 已經(jīng)有5人回復(fù)
280求調(diào)劑已經(jīng)有12人回復(fù)
269專碩求調(diào)劑已經(jīng)有5人回復(fù)

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

查看全部散金貼

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

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

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的回帖

lijinfeng042

木蟲 (小有名氣)

引用回帖:

程序啊數(shù)值計算的書都有啊需要的話你搜一下我發(fā)的帖子就有

贊一下

回復(fù)此樓

5樓2011-01-04 19:49:24

已閱回復(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的回帖

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 huameitang05 的主題更新

返回列表

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

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

[考研] 求調(diào)劑院校信息 +4	CX 330 2026-03-21	4/200	2026-03-21 23:48 by ms629
[考研] 085600材料與化工306 +4	z1z2z3879 2026-03-21	4/200	2026-03-21 23:44 by ms629
[考研] 材料學(xué)碩301分求調(diào)劑 +7	Liyouyumairs 2026-03-21	7/350	2026-03-21 22:31 by peike
[考研] 材料 271求調(diào)劑 +5	展信悅_ 2026-03-21	5/250	2026-03-21 17:29 by 學(xué)員8dgXkO
[考研] 336求調(diào)劑 +5	rmc8866 2026-03-21	5/250	2026-03-21 17:24 by 學(xué)員8dgXkO
[考研] 277材料科學(xué)與工程080500求調(diào)劑 +6	自由煎餅果子 2026-03-16	6/300	2026-03-21 17:21 by 學(xué)員8dgXkO
[考研] 266求調(diào)劑 +3	哇呼哼呼哼 2026-03-20	3/150	2026-03-21 16:46 by barlinike
[考研] 能源材料化學(xué)課題組招收碩士研究生8-10名 +5	脫穎而出 2026-03-16	15/750	2026-03-21 10:16 by 脫穎而出
[考研] 303求調(diào)劑 +5	睿08 2026-03-17	7/350	2026-03-21 03:11 by JourneyLucky
[考研] 一志愿中國石油大學(xué)（華東）本科齊魯工業(yè)大學(xué) +3	石能偉 2026-03-17	3/150	2026-03-21 02:22 by JourneyLucky
[考研] 332求調(diào)劑 +4	ydfyh 2026-03-17	4/200	2026-03-21 02:20 by JourneyLucky
[考研] 一志愿中海洋材料工程專碩330分求調(diào)劑 +8	小材化本科 2026-03-18	8/400	2026-03-20 23:16 by JourneyLucky
[考研] 考研調(diào)劑求學(xué)校推薦 +3	伯樂29 2026-03-18	5/250	2026-03-20 22:59 by JourneyLucky
[考研] 295復(fù)試調(diào)劑 +8	簡木ChuFront 2026-03-19	8/400	2026-03-20 20:44 by zhukairuo
[考研] 295材料求調(diào)劑，一志愿武漢理工085601專碩 +5	Charlieyq 2026-03-19	5/250	2026-03-20 20:35 by JourneyLucky
[考研] 312求調(diào)劑 +8	陌宸希 2026-03-16	9/450	2026-03-18 12:39 by Linda Hu
[考博] 26申博 +4	八6八68 2026-03-16	4/200	2026-03-17 13:00 by 輕松不少隨
[論文投稿] 有沒有大佬發(fā)小論文能帶我個二作 +3	增銳漏人 2026-03-17	4/200	2026-03-17 09:26 by xs74101122
[考研] 318求調(diào)劑 +3	Yanyali 2026-03-15	3/150	2026-03-16 16:41 by houyaoxu
[考研] 327求調(diào)劑 +6	拾光任染 2026-03-15	11/550	2026-03-15 22:47 by 拾光任染

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

24小時熱門版塊排行榜

huameitang05

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

» 猜你喜歡

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

lijinfeng042

shl1025

ysh2225

lijinfeng042

honeypirl

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

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

查看全部散金貼