1个回答
展开全部
如果知道A,w的大概范围,可以这样:
function y=da59(z) % 保存为 da59.m 文件
A=z(1);w=z(2);
syms x
y=A*w*sin(w*x);
f=sqrt(1+y^2);
s=int(f,x,0,pi/(2*w));
y=double([s-(pi*A/2);w*A-15/sqrt(A^2-225)]);
----------------------------------------------------------
运行以下:
Aw0=[~,~] % 初值
[Aw,fval,exitflag]=fsolve('da59',Aw0)
function y=da59(z) % 保存为 da59.m 文件
A=z(1);w=z(2);
syms x
y=A*w*sin(w*x);
f=sqrt(1+y^2);
s=int(f,x,0,pi/(2*w));
y=double([s-(pi*A/2);w*A-15/sqrt(A^2-225)]);
----------------------------------------------------------
运行以下:
Aw0=[~,~] % 初值
[Aw,fval,exitflag]=fsolve('da59',Aw0)
追问
兄弟 懂了 !谢谢了,最后一句不懂意思,能否解释一下?
追答
% Aw 为所求的解
%fval为解代入方程的验证结果
% exitflag=1表示函数收敛于解;反之,=0
推荐律师服务:
若未解决您的问题,请您详细描述您的问题,通过百度律临进行免费专业咨询