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点求根要用什么函数呢？
...