 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 6 | 1/1 | 返回列表 |
| 查看: 676 | 回復(fù): 5 | |||
[交流]
【求助】急求幫忙解決MATLAB 樣條差值求解微分代碼出現(xiàn)的問題，使代碼正常運(yù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ù)此樓» 猜你喜歡
288求調(diào)劑
已經(jīng)有5人回復(fù)
-
281求調(diào)劑（0805）
已經(jīng)有19人回復(fù)
-
材料考研調(diào)劑
已經(jīng)有3人回復(fù)
-
346求調(diào)劑[0856]
已經(jīng)有6人回復(fù)
-
294求調(diào)劑材料與化工專碩
已經(jīng)有9人回復(fù)
-
328求調(diào)劑，英語六級(jí)551，有科研經(jīng)歷
已經(jīng)有12人回復(fù)
-
一志愿西安交通大學(xué)材料工程專業(yè) 282分求調(diào)劑
已經(jīng)有5人回復(fù)
-
085601材料工程專碩求調(diào)劑
已經(jīng)有7人回復(fù)
-
321求調(diào)劑
已經(jīng)有6人回復(fù)
-
一志愿中海洋材料工程專碩330分求調(diào)劑
已經(jīng)有7人回復(fù)
» 搶金幣啦！回帖就可以得到:
-
散金幣，求好運(yùn)，祝面上順利！
+2/84
-
【代人發(fā)帖勿站內(nèi)信】【成都征男友】89年未婚博士漂亮女孩
+1/69
-
湘潭大學(xué)化工學(xué)院國家級(jí)人才團(tuán)隊(duì)2026年博士研究生招生
+1/64
-
青島高校 研究生招生
+1/50
-
催化方向推薦一個(gè)極好的博導(dǎo)
+1/35
-
澳大利亞西澳大學(xué)（UWA）張金強(qiáng)課題組招收全獎(jiǎng)博士生（化工新能源方向）
+1/32
-
天津工業(yè)大學(xué)航空航天博士招生
+1/32
-
麗水學(xué)院 材料與能源工程學(xué)院 陳江飛老師課題組碩士招生----材料與化工專業(yè)
+1/25
-
江西贛州-贛南師范大學(xué)-學(xué)校推薦
+1/16
-
【博士招生】武漢科技大學(xué)招收材料、化工、環(huán)境、冶金類2026年“申請(qǐng)-考核”制博士生
+1/8
-
新加坡國立大學(xué)藥學(xué)系化學(xué)生物學(xué)課題組招PhD
+1/8
-
南京理工大學(xué)2026級(jí)碩士生招生-張光普博導(dǎo) &碩導(dǎo)
+2/8
-
教育部長江學(xué)者和創(chuàng)新團(tuán)隊(duì)發(fā)展計(jì)劃”入選團(tuán)隊(duì)招收 材料，化學(xué)與化工博士研究生
+1/7
-
國家級(jí)人才團(tuán)隊(duì)誠聘生物信息相關(guān)專業(yè)的科研助理
+1/5
-
MTI 380調(diào)劑
+1/5
-
武漢紡織大學(xué)杰青團(tuán)隊(duì)招研究生（材料，化學(xué)，高分子，化工）
+1/4
-
中山大學(xué)醫(yī)學(xué)院趙博課題組碩士招生
+1/3
-
課題組招收環(huán)境及相關(guān)專業(yè)調(diào)劑碩士研究生（歡迎優(yōu)秀學(xué)生加入）
+1/2
-
Top-88悉尼科技大學(xué)數(shù)據(jù)科學(xué)/AI 招收2027年入學(xué) 校獎(jiǎng) 博士生1到2名(國際和本地學(xué)生）
+1/2
-
浙江大學(xué)高純分離團(tuán)隊(duì)誠聘博士后
+1/2
lijinfeng042
木蟲 (小有名氣)
- 仿真EPI: 2
- 應(yīng)助: 1 (幼兒園)
- 金幣: 2156.1
- 帖子: 291
- 在線: 31.5小時(shí)
- 蟲號(hào): 1019062
★ ★ ★ ★ ★
huameitang05(金幣+1):謝謝參與
huameitang05(金幣+1):謝謝意見 2011-01-04 19:14:27
xiegangmai(金幣+2):辛苦了！ 2011-01-07 22:06:47
robert2020(金幣+2):辛苦了！最近比較忙，評(píng)帖不及時(shí)，敬請(qǐng)見諒！ 2011-01-11 10:47:08
huameitang05(金幣+1):謝謝參與
huameitang05(金幣+1):謝謝意見 2011-01-04 19:14:27
xiegangmai(金幣+2):辛苦了！ 2011-01-07 22:06:47
robert2020(金幣+2):辛苦了！最近比較忙，評(píng)帖不及時(shí)，敬請(qǐng)見諒！ 2011-01-11 10:47:08
|
具體太多代碼沒看...只是說說思路 離散數(shù)據(jù)的微分問題 既然你用的是樣條插值 那就說說這個(gè)方法 數(shù)據(jù)->csapi(x,y)->fnder(cs)->fnval(pp,x) |
2樓2011-01-04 19:06:32
3樓2011-01-04 19:31:42
4樓2011-01-04 19:47:16
lijinfeng042
木蟲 (小有名氣)
- 仿真EPI: 2
- 應(yīng)助: 1 (幼兒園)
- 金幣: 2156.1
- 帖子: 291
- 在線: 31.5小時(shí)
- 蟲號(hào): 1019062
5樓2011-01-04 19:49:24
6樓2011-02-28 21:36:36
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 huameitang05 的主題更新
| 6 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 346求調(diào)劑[0856] +3 | WayneLim327 2026-03-16 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 本科鄭州大學(xué)物理學(xué)院，一志愿華科070200學(xué)碩，346求調(diào)劑 +4 | 我不是一根蔥 2026-03-18 | 4/200 |
|
|
|
[考研] 354求調(diào)劑 +4 | Tyoumou 2026-03-18 | 7/350 |
|
|
|
[考研] 一志愿中國海洋大學(xué)，生物學(xué)，301分，求調(diào)劑 +4 | 1孫悟空 2026-03-17 | 4/200 |
|
|
|
[考研] 311求調(diào)劑 +6 | 26研0 2026-03-15 | 6/300 |
|
|
|
[考研] 材料專碩274一志愿陜西師范大學(xué)求調(diào)劑 +6 | 薛云鵬 2026-03-13 | 6/300 |
|
|
|
[考研] 工科材料085601 279求調(diào)劑 +6 | 困于星晨 2026-03-17 | 6/300 |
|
|
|
[考研] 生物學(xué)071000 329分求調(diào)劑 +3 | 我愛生物生物愛?/a> 2026-03-17 | 3/150 |
|
|
|
[考研] 268求調(diào)劑 +6 | 簡單點(diǎn)0 2026-03-17 | 6/300 |
|
|
|
[考研] 考研化學(xué)學(xué)碩調(diào)劑，一志愿985 +4 | 張vvvv 2026-03-15 | 6/300 |
|
|
|
[考研] 材料專碩326求調(diào)劑 +6 | 墨煜姒莘 2026-03-15 | 7/350 |
|
|
|
[考研] 本人考085602 化學(xué)工程 專碩 +16 | 不知道叫什么！ 2026-03-15 | 18/900 |
|
|
|
[考研] 一志愿南京大學(xué)，080500材料科學(xué)與工程，調(diào)劑 +4 | Jy? 2026-03-16 | 4/200 |
|
|
|
[考研] 考研調(diào)劑 +3 | 淇ya_~ 2026-03-17 | 5/250 |
|
|
|
[考研] 274求調(diào)劑 +5 | 時(shí)間點(diǎn) 2026-03-13 | 5/250 |
|
|
|
[考研] [導(dǎo)師推薦]西南科技大學(xué)國防/材料導(dǎo)師推薦 +3 | 尖角小荷 2026-03-16 | 6/300 |
|
|
|
[考研] 277材料科學(xué)與工程080500求調(diào)劑 +3 | 自由煎餅果子 2026-03-16 | 3/150 |
|
|
|
[基金申請(qǐng)] 現(xiàn)在如何回避去年的某一個(gè)專家，不知道名字 +3 | zk200107 2026-03-12 | 6/300 |
|
|
|
[考研] 材料301分求調(diào)劑 +5 | Liyouyumairs 2026-03-12 | 5/250 |
|
|
|
[考研] 0817化學(xué)工程與技術(shù)考研312分調(diào)劑 +3 | T123 tt 2026-03-12 | 3/150 |
|