Matlab二维坐标下极坐标向直角坐标的转换

今天要处理一些极坐标的数据，把他们转化成为直角坐标，本来这事不复杂，在python里面就一个命令的事情就可以把储存在极坐标下(矩阵大小M*N)的数据转化为任意大小的矩阵来储存直角坐标下的数据，其原理就是简单的线性插值。

http://liao.cpython.org/scipytutorial11/

但是坑爹的matlab的插值命令girddata(x,y,u,X,Y)要求x,y,u必须是一维数组，X,Y必须是二维数组，导致不像python一样直接能转换。写了一个简单的循环算法来把极坐标数据转换为直角坐标，用了循环嵌套和计数器，所以代码效率极低，但是也不是很想去优化，能用就行了。放在这里给大家借鉴，欢迎提出代码修改意见。

clear all
data = csvread('REP_2Dvelocity.csv');
r1 = reshape(data(:,1),[201,51])';
ux1 = reshape(data(:,3),[201,51])';
r1 = reshape(data(:,1),[201,51</