MATLAB数据可视化学习笔记

地形地貌图形绘制

%%地形地貌图形绘制
[x,y]=meshgrid(1:10);%构造测量网格
h=[0,0.02,-0.12,0,-2.09,0,-0.58,-0.08,0,0;
    0.02,0,0,-2.38,0,-4.96,0,0,0,-0.1;
    0,0.1,1,0,-3.04,0,-0.53,0,0.1,0;
    0,0,0,3.52,0,0,0,0,0,0;
    -0.43,-1.98,0,0,0,0.77,0,2.17,0,0;
    0,0,-2.29,0,0.69,0,2.59,0,0.3,0;
    -0.09,-0.31,0,0,0,4.27,0,0,0,-0.01;
    0,0,0,5.13,7.4,0,1.89,0,0.04,0;
    0.1,0,0.68,0,0,1.75,0,-0.11,0,0;
    0,-0.01,0,0,0.3,0,0,0,0,0.01];%测量数据点
[xi,yi]=meshgrid(1:0.1:10);%构造插值网格
hi=interp2(x,y,h,xi,yi,'spline');%二维插值命令
surf(hi);%画出地貌图
xlabel('x');ylabel('y');zlabel('h');

车灯光源投影区域的绘制

CUMCM 2002 A题 车灯线光源的优化设计

安装在汽车头部的车灯的形状为一旋转抛物面，车灯的对称轴水平地指向正前方, 其开口半径36毫米，深度21.6毫米。经过车灯的焦点，在与对称轴相垂直的水平方向，对称地放置一定长度的均匀分布的线光源。要求在某一设计规范标准下确定线光源的长度。
该设计规范在简化后可描述如下。在焦点F正前方25米处的A点放置一测试屏，屏与FA垂直，用以测试车灯的反射光。在屏上过A点引出一条与地面相平行的直线，在该直线A点的同侧取B点和C点，使AC=2AB=2.6米。要求C点的光强度不小于某一额定值（可取为1个单位），B点的光强度不小于该额定值的两倍（只须考虑一次反射）。
请解决下列问题：
（1）在满足该设计规范的条件下，计算线光源长度，使线光源的功率最小。
（2）对得到的线光源长度，在有标尺的坐标系中画出测试屏上反射光的亮区。
（3）讨论该设计规范的合理性。

%%车灯光源投影区域的绘制
p=0.03;x=25.0216;
for y1=-0.002:0.0004:0.002
    y0=(-0.036:0.001:0.036)'*ones(1,73);
    z0=ones(73,1)*(-0.036:0.001:0.036);
    x0=(y0.^2+z0.^2)/(2*p);
    xn=(p^3+4*x0*2*p.*x0+p*(-4*y1*y0+3*2*p*x0))./(2*(p^2+2*p*x0));
    yn=(2*p*x0.*y0+p^2*(-y1+y0)+y1*(y0.^2-z0.^2))./(p^2+2*p*x0);
    zn=(p^2+2*p*x0+2*y1*y0).*z0./(p^2+2*p*x0);
    y=y0+(yn-y0).*(x-x0)./(xn-x0);
    z=z0+(zn-z0).*(x-x0)./(xn-x0);
    plot(y,z,'b.')
    xlabel('y')
    ylabel('z')
    hold on
end

学习来源：MATLAB在数学建模中的应用