%人字架优化问题-作图法求解 20140914作图
% 2015年3月31日晚23点添加坐标轴
clear all;
close all;
clc;
F=150000;
B=0.76;
p=7800;
T=0.0025;
%sgm=420*10^6;clear all;
close all;
E=2.1*10^11;
B=0.76;
p=7800;
T=0.0025;
%sgm=420*10^6;
% 参考程序
%figure1 = figure;
% annotation1 = annotation(figure1,'arrow',[0.131 0.131],[0.92 0.96]);
% annotation2 = annotation(figure1,'arrow',[0.88 0.96],[0.108 0.108]);
figure1 = figure;
annotation1 = annotation(figure1,'arrow',[0.13 0.13],[0.91 0.99]);%经过精细调整，该数值比网上效果好（但要承认借鉴的他们的优秀思想）。
annotation2 = annotation(figure1,'arrow',[0.88 0.96],[0.109 0.109]);%经过精细调整，该数值比网上效果好（但要承认借鉴的他们的优秀思想）。
syms D h c1;
%ezplot('sqrt(B^2+h^2)*F/(pi*T*D*h)-pi^2*E*(T^2+D^2)/(8*(B^2+h^2))')
%ezplot('sqrt(0.76^2+h^2)*150000/(pi*0.0025*D*h)-pi^2*2.1*10^11*(0.0025^2+D
%^2)/(8*(0.76^2+h^2))',[0,0.16,0,1.6]),set(findobj('Type','line'),'LineStyle','*','color','y') plot(D,h,'*')
c1=ezplot('sqrt(0.76^2+h^2)*150000/(pi*0.0025*D*h)-pi^2*2.1*10^11*(0.0025^2+D^2)/(8*(0.76^2+h^2))',[0,0.16,0,1.6]);%要注意单位
set(c1 ,'color','b','linewidth',2);
hold on
ezplot('2*pi*7800*0.0025*D*(0.76^2+h^2)^0.5-10')
hold on
ezplot('2*pi*7800*0.0025*D*(0.76^2+h^2)^0.5-8.47')
set(D ,'color','m','linewidth',2);
hold on
plot(h,D,'color',[0 1 0])
set(findobj('Type','line'),'LineStyle','--','color','r')
hold on
ezplot('150000*(0.76^2+h^2)^0.5/(pi*0.0025*D*h)-420*10^6',[0,0.16,0,1.6]),set(findobj('Type','line'),'LineStyle','--','color','k')
hold on
ezplot('150000*(0.76^2+h^2)^0.5/(pi*0.0025*D*h)-703*10^6',[0,0.16,0,1.6])
patch([D_range, fliplr(D_range)], [h_lower, fliplr(h_upper)], 'm')
alpha(0.65)
legend('Extended DynamiCS','DynamiCS')
%ezplot('150000*(0.76^2+h^2)^0.5/(pi*0.0025*D*h)-703*10^6',[0,0.16,0,1.6]),set(findobj('Type','line'),'LineStyle','--','color','g')
