雨人007

從MATLAB中找出的源程序，直接不懂！
>>type griddata
function [xi,yi,zi] = gdatav4(x,y,z,xi,yi)
%GDATAV4 MATLAB 4 GRIDDATA interpolation

% Reference:  David T. Sandwell, Biharmonic spline
% interpolation of GEOS-3 and SEASAT altimeter
% data, Geophysical Research Letters, 2, 139-142,
% 1987.  Describes interpolation using value or
% gradient of value in any dimension.

xy = x( + y(*sqrt(-1);

% Determine distances between points
d = xy(:,ones(1,length(xy)));
d = abs(d - d.');
n = size(d,1);
% Replace zeros along diagonal with ones (so these don't show up in the
% find below or in the Green's function calculation).
d(1:n+1:numel(d)) = ones(1,n);

non = find(d == 0, 1);
if ~isempty(non),
  % If we've made it to here, then some points aren't distinct.  Remove
  % the non-distinct points by averaging.
  [r,c] = find(d == 0);
  k = find(r < c);
  r = r(k); c = c(k); % Extract unique (row,col) pairs
  v = (z(r) + z(c))/2; % Average non-distinct pairs

  rep = find(diff(c)==0);
  if ~isempty(rep), % More than two points need to be averaged.
runs = find(diff(diff(c)==0)==1)+1;
for i=1:length(runs),
   k = (c==c(runs(i))); % All the points in a run
   v(runs(i)) = mean(z([r(k);c(runs(i))])); % Average (again)
end
  end
  z(r) = v;
  if ~isempty(rep),
z(r(runs)) = v(runs); % Make sure average is in the dataset
  end

  % Now remove the extra points.
  z(c) = [];
  xy(c, = [];
  xy(:,c) = [];
  d(c, = [];
  d(:,c) = [];

  % Determine the non distinct points
  ndp = sort([r;c]);
  ndp(ndp(1:length(ndp)-1)==ndp(2:length(ndp))) = [];

  warning('MATLAB:griddata:NonDistinctPoints',['Averaged %d non-distinct ' ...
         'points.\n       Indices are: %s.'],length(ndp),num2str(ndp'))
end

% Determine weights for interpolation
g = (d.^2) .* (log(d)-1); % Green's function.
% Fixup value of Green's function along diagonal
g(1:size(d,1)+1:numel(d)) = zeros(size(d,1),1);
weights = g \ z(;

[m,n] = size(xi);
zi = zeros(size(xi));
jay = sqrt(-1);
xy = xy.';

% Evaluate at requested points (xi,yi).  Loop to save memory.
for i=1:m
  for j=1:n
d = abs(xi(i,j)+jay*yi(i,j) - xy);
mask = find(d == 0);
if ~isempty(mask), d(mask) = ones(length(mask),1); end
g = (d.^2) .* (log(d)-1); % Green's function.
% Value of Green's function at zero
if ~isempty(mask), g(mask) = zeros(length(mask),1); end
zi(i,j) = g * weights;
  end
end

if nargout<=1,
  xi = zi;
end

贊一下

回復(fù)此樓

想做的事情擋不�。�

5樓2014-12-25 10:22:14

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

Kevlar_J

銅蟲(chóng) (正式寫(xiě)手)

應(yīng)助: 5 (幼兒園)
金幣: 297.2
散金: 85
紅花: 5
帖子: 366
在線(xiàn): 116.9小時(shí)
蟲(chóng)號(hào): 2697572
注冊(cè): 2013-10-04
性別: GG
專(zhuān)業(yè): 電磁場(chǎng)與電路

引用回帖:

5樓: Originally posted by 雨人007 at 2014-12-25 10:22:14
從MATLAB中找出的源程序，直接不懂！
>>type griddata
function = gdatav4(x,y,z,xi,yi)
%GDATAV4 MATLAB 4 GRIDDATA interpolation

% Reference: David T. Sandwell, Biharmonic spline
% int ...

這說(shuō)的我也不懂，我覺(jué)得一句句翻譯成數(shù)學(xué)語(yǔ)言，應(yīng)該可以理解的，別急，慢慢看

[ 發(fā)自小木蟲(chóng)客戶(hù)端 ]

贊一下

回復(fù)此樓

振葉以尋根，觀瀾而溯源！www.source-grid.com

6樓2014-12-25 12:11:43

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

Kevlar_J

銅蟲(chóng) (正式寫(xiě)手)

應(yīng)助: 5 (幼兒園)
金幣: 297.2
散金: 85
紅花: 5
帖子: 366
在線(xiàn): 116.9小時(shí)
蟲(chóng)號(hào): 2697572
注冊(cè): 2013-10-04
性別: GG
專(zhuān)業(yè): 電磁場(chǎng)與電路

引用回帖:

有什么不懂到matlab愛(ài)好者論壇問(wèn)問(wèn)，哪里大神多

[ 發(fā)自小木蟲(chóng)客戶(hù)端 ]

贊一下

回復(fù)此樓

振葉以尋根，觀瀾而溯源！www.source-grid.com

7樓2014-12-25 12:12:29

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

雨人007

銀蟲(chóng) (小有名氣)

應(yīng)助: 12 (小學(xué)生)
金幣: 247.5
帖子: 212
在線(xiàn): 43.7小時(shí)
蟲(chóng)號(hào): 2672857
注冊(cè): 2013-09-22
性別: GG
專(zhuān)業(yè): 資源化工

引用回帖:

7樓: Originally posted by Kevlar_J at 2014-12-25 12:12:29
有什么不懂到matlab愛(ài)好者論壇問(wèn)問(wèn)，哪里大神多
...

多謝�。�

回復(fù)此樓

想做的事情擋不��！

8樓2014-12-26 10:02:01

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

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

[考研] 【雙一流院校新能源、環(huán)境材料，材料加工與模擬招收大量調(diào)劑】 +4	Higraduate 2026-03-22	7/350	2026-03-24 11:23 by 種大樹(shù)
[考研] 298-一志愿中國(guó)農(nóng)業(yè)大學(xué)-求調(diào)劑 +11	手機(jī)用戶(hù) 2026-03-17	12/600	2026-03-23 23:51 by 熱情沙漠
[考研] 341求調(diào)劑(一志愿湖南大學(xué)070300) +5	番茄頭--- 2026-03-22	6/300	2026-03-23 23:45 by Txy@872106
[考研] 材料292調(diào)劑 +4	橘頌思美人 2026-03-23	4/200	2026-03-23 23:16 by peike
[考研] 一志愿211 初試270分求調(diào)劑 +4	谷雨上岸 2026-03-23	5/250	2026-03-23 21:18 by 不惑可樂(lè)
[考研] 317求調(diào)劑 +12	申子申申 2026-03-19	18/900	2026-03-22 22:23 by luoyongfeng
[考研] 生物學(xué)071000 329分求調(diào)劑 +5	我愛(ài)生物生物愛(ài)?/a> 2026-03-17	5/250	2026-03-22 16:42 by tcx007
[考研] 285求調(diào)劑 +6	ytter 2026-03-22	6/300	2026-03-22 12:09 by 星空星月
[基金申請(qǐng)] 山東省面上項(xiàng)目限額評(píng)審 +4	石瑞0426 2026-03-19	4/200	2026-03-22 08:50 by Wei_ren
[考研] 286分人工智能專(zhuān)業(yè)請(qǐng)求調(diào)劑愿意跨考！ +4	lemonzzn 2026-03-17	8/400	2026-03-21 22:49 by lemonzzn
[考研] 0703化學(xué)調(diào)劑 +4	妮妮ninicgb 2026-03-21	4/200	2026-03-21 18:39 by 學(xué)員8dgXkO
[考研] 297求調(diào)劑 +3	喜歡還是不甘心 2026-03-20	3/150	2026-03-21 18:33 by 學(xué)員8dgXkO
[考研] 278求調(diào)劑 +9	煙火先于春 2026-03-17	9/450	2026-03-21 17:47 by 學(xué)員8dgXkO
[考研] 313求調(diào)劑 +4	肆叁貳壹22 2026-03-19	4/200	2026-03-21 17:33 by ColorlessPI
[考研] 求調(diào)劑 +3	Ma_xt 2026-03-17	3/150	2026-03-21 02:05 by JourneyLucky
[考研] 324分 085600材料化工求調(diào)劑 +4	llllkkkhh 2026-03-18	4/200	2026-03-21 01:24 by JourneyLucky
[考研] 一志愿中海洋材料工程專(zhuān)碩330分求調(diào)劑 +8	小材化本科 2026-03-18	8/400	2026-03-20 23:16 by JourneyLucky
[考研] 一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑 +11	Ymlll 2026-03-18	15/750	2026-03-20 19:40 by 丁丁*
[考研] 求調(diào)劑 +3	@taotao 2026-03-20	3/150	2026-03-20 19:35 by JourneyLucky
[論文投稿] 申請(qǐng)回稿延期一個(gè)月，編輯同意了。但系統(tǒng)上的時(shí)間沒(méi)變，給編輯又寫(xiě)郵件了，沒(méi)回復(fù) 10+3	wangf9518 2026-03-17	4/200	2026-03-19 23:55 by babero