[Scilab-users] filled contour plot

[Rafael Guerra]

Hi Claus,

Pls check if this solution is satisfactory:


//START OF CODE

frequencies = [200 2000 20000]; // Hz

angles = [0 30 60 90 120 150 180]; // degrees - convert to radians!

mag00 = [69.253448 92.163385 87.414209]; // on-axis freq-response

mag30 = [68.461357 90.569082 76.073549];

mag60 = [65.842262 90.745621 74.152067];

mag90 = [59.420658 80.893365 55.975304];

mag120 = [52.505173 80.795915 44.363582];

mag150 = [53.829987 59.793715 42.443964];

mag180 = [55.650341 73.457973 55.608187];



magnitudes = [mag00' mag30' mag60' mag90' mag120' mag150' mag180'];



negative_angles = -angles(2:7);

neg_angle_magnitudes = magnitudes(:,2:7);

angles = [angles negative_angles];

magnitudes = [magnitudes neg_angle_magnitudes];

[angles idx] = gsort(angles,'g','i');

magnitudes = magnitudes(:,idx);



clf;

contourplt = scf();

contourplt.color_map=jetcolormap(6); // blue = cold, red = hot :-)

nz = 40:10:90;

contour2d(frequencies,angles,magnitudes,nz,logflag="ln");

contourf(frequencies,angles,magnitudes,nz);

legends(string(nz),1:length(nz),"lr");

xlabel("frequency [Hz]");

ylabel("angle (deg)");

//END OF CODE

Regards,
Rafael


-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.scilab.org/pipermail/users/attachments/20180121/0e5fdf9d/attachment.htm>