using-MATLAB.zip_人字架

%人字架优化问题-作图法求解 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')