现在有一组风速风向数据(u1,u2),u1是风速0m/s-40m/s,u2是风向0°-360°
风速用的是gamma分布拟合的
风向用的是four-component the mixture of von Mises分布拟合的
如何用matlab只带的frank copula函数,生成下面左边的图形,右边是原数据的概率图
我试着用了matlab帮助里面的copulapdf例子(下面代码是帮助里的例子),他那个必须u1和u2都是同样的数,我把风速数据带进去就全为0,可能是我理解有问题。
u = linspace(0,1,10);
[u1,u2] = meshgrid(u,u);
y = copulapdf('Frank',[u1(:),u2(:)],1);
surf(u1,u2,reshape(y,10,10))
xlabel('u1')
ylabel('u2')