matlab 图片倒影_matlab-画个拱桥和倒影？

matlab可用于各行各业的应用中，现在我们就简单画一个拱桥试一试。

r=input('请输入半径: ');

z=input('请输入弧度: ');

figure()

%r=4;

%画拱洞，用李萨如图形法

x=15;

y=0;

nseg=1000;

theta = 0 : (2 * pi / nseg) : (2 * pi);

pline_x = r * cos(theta) + x;

pline_y = r * sin(theta) + y;

plot(pline_x, pline_y);

axis([0 30 -8 8]); %限定坐标范围

hold on;

%画水面

x1=0:0.1:30;

y1=zeros(size(x1));

plot(x1,y1)

hold on;

%画测量数据

xg=0:2:30;

yg=[4,4.2,4.6,4.8,5,5.3,6.4,6.9,6.7,6.1,5.6,5.1,4.7,4.4,4,3.6];

plot(xg,yg,'o')

hold on

%画拱桥

xxg=-100:0.1:100;

p=polyfit(xg,yg,2); % 拟合出的二次函数的系数

yyg=polyval(p,xxg); % 拟合曲线的纵坐标

xxg=xxg(yyg>=0);

yyg=yyg(yyg>=0);

%画倒影

plot(xxg,yyg)

hold on

plot(xxg,-yyg)

hold on

fill(xxg,-yyg,'b')

fill(pline_x, pline_y,'w')

%计算面积

xxg2=0:0.1:30;

yyg2=yyg(find(xxg==0):find(xxg==30));

%polyarea(xxg2,yyg)

Aera=trapz(xxg2,yyg2)-pi/2*r^2;

disp(['面试是： ',num2str(Aera)])

%进行安全度分析

if(z>0.1)

disp('经过安全性分析，该拱桥在应用中可能存在安全隐患。')

else

disp('经过安全性分析，该拱桥的安全性良好，可以使用。')

end

那么最后的结果是