help me in using "interp2"

2015 9 19
1 个回答
回答已采纳
更新时间：2015 9 20
3 次查看（30 天）

0 个投票

I want to expand my matrix using interpolation. I followed the instructions but it shows some errors:
clc
clear
rw=0.10795;
n = linspace(-pi,pi,30);
m = linspace(rw,2.5,30);
[theta,r] = meshgrid(n,m);
[X,Y] = pol2cart(theta,r);
V = peaks(X,Y);
figure
surf(X,Y,V)
title('Original Sampling');
nn = linspace(-pi,pi,360);
mm = linspace(rw,2.5,200);
[thetaaa,rrr] = meshgrid(nn,mm);
[Xq,Yq] = pol2cart(thetaaa,rrr);
Vq = interp2(X,Y,V,Xq,Yq,'cubic');
figure
surf(Xq,Yq,Vq);
title('Cubic Interpolation Over Finer Grid');

请先登录，再回答此问题。

 采纳的回答

Walter Roberson
Walter Roberson 2015-9-19

1 个投票

your input grid, X, Y, must be cuboid.
Leave your inputs the same, but change the command from interp2 to griddata (with the same arguments).
Also you will want to add 'edgecolor','none' to the surf() call.
griddata() is going to tell you that duplicate values have been averaged. I suspect -pi and +pi map to the same point with pol2cart. If you were to linspace() with 31 points and then delete the last then the overlap should disappear.

3 个评论
显示 1更早的评论 隐藏 1更早的评论

Mahmood
Mahmood 2015-9-19
编辑：Walter Roberson 2015-9-19
Thanks a lot for your help, but it still not working. I correct the -pi and cuboid but I don't think the cubic method is the problem. It is the method for interpolation. Can you please present the corrected code to me? thanks a lot
clc
clear
rw=0.10795;
n = linspace(-pi,pi,31);
n(1) = [];
m = linspace(rw,2.5,30);
[theta,r] = meshgrid(n,m);
[X,Y] = pol2cart(theta,r);
V = peaks(X,Y);
figure
surf(X,Y,V)
title('Original Sampling');
nn = linspace(-pi,pi,360);
nn(1) = [];
mm = linspace(rw,2.5,200);
[thetaaa,rrr] = meshgrid(nn,mm);
[Xq,Yq] = pol2cart(thetaaa,rrr);
Vq = interp2(X,Y,V,Xq,Yq);
figure
surf(Xq,Yq,Vq);
Walter Roberson
Walter Roberson 2015-9-19
编辑：Walter Roberson 2015-9-19
rw=0.10795;
n = linspace(-pi,pi,31);
n(1) = [];
m = linspace(rw,2.5,30);
[theta,r] = meshgrid(n,m);
[X,Y] = pol2cart(theta,r);
V = peaks(X,Y);
figure
surf(X,Y,V, 'edgecolor','none')
title('Original Sampling');
nn = linspace(-pi,pi,360);
mm = linspace(rw,2.5,200);
[thetaaa,rrr] = meshgrid(nn,mm);
[Xq,Yq] = pol2cart(thetaaa,rrr);
Vq = griddata(X,Y,V,Xq,Yq,'cubic');
figure
surf(Xq,Yq,Vq, 'edgecolor', 'none');
title('Cubic Interpolation Over Finer Grid');
Mahmood
Mahmood 2015-9-20
Thanks a lot for your help. I really appreciate that.

请先登录，再进行评论。

更多回答（0 个）

类别

在 帮助中心 和 File Exchange 中查找有关 Interpolation 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by