function jie_fangcheng
clc;
clear all;
close all;
format long
xx1=[linspace(0.1048,0.1256,1000)];
xx2=[linspace(-0.1256,-0.1048,1000)];
yy1=obj(xx1);
a1=0.1*(max(xx1)-min(xx1));
b1=0.1*(max(yy1)-min(yy1));

xx2=[linspace(-0.1256,-0.1048,1000)];
yy2=obj(xx2);
a2=0.1*(max(xx2)-min(xx2));
b2=0.1*(max(yy2)-min(yy2));

figure(1)
subplot(2,1,1)
plot([min(xx1)-a1 max(xx1)+a1],[0 0],'r-',xx1,yy1,'b-');axis([min(xx1)-a1 max(xx1)+a1 min(yy1)-b1 max(yy1)+b1]);


subplot(2,1,2)
plot([min(xx2)-a2 max(xx2)+a2],[0 0],'r-',xx2,yy2,'b-');axis([min(xx2)-a2 max(xx2)+a2 min(yy2)-b2 max(yy2)+b2]);


[x fval]=fsolve(@obj,[ -0.105 -0.12 0.105 0.12  ])   

function f=obj(x)

cl=1500;cs=3000;bs=1800;
pl=1000;ps=2500;H=50;h=25;f=30;w=2*pi*f;
kal=w/cl;
kas=w/cs;
kbs=w/bs;

sln=sqrt(kal^2-x.^2);
ssn=sqrt(x.^2-kas^2);
lsn=sqrt(x.^2-kbs^2);
f=(ps/pl)*((x/kbs).^4).*(sln./ssn).*(4*ssn.*ssn./(x.^2)-((x.^2+lsn)./x.^2).^2)-tan(sln*H);

2Â¥: Originally posted by Mr__Right at 2016-05-03 20:32:38
1. xÈ«¶¼ÊÇÒÔ2´ÎÐÎʽ³öÏֵģ¬¼ò»¯µÄ»°£¬²»ÈçÏÈÇóx^2;
2. ·ûºÅ½â¿ÉÄÜÊÇÇó²»³öÀ´µÄ£¬Çó¹ØÓÚ t=x^2 µÄ·ÇÏßÐÔ·½³Ì F(t)=0µÄ¸ùµÄʱºò£¬ÏÈ»­Í¼³öÀ´¿ÉÒÔ¿´µ½ÔÚ(0.01,0.016)Ö®¼ä tÓжà¸ù

3Â¥: Originally posted by ÉÏÌ«¿Õ˵ at 2016-05-03 23:21:54
¶àлÄú°¡£¬¶ÔÓÚµÚ2µãÇó¸ùÒªÓÃʲôº¯ÊýÄØ£¿
...