马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?注册
x
问题如上,本人曾用int函数求解,但最终结果与应得结果相差很大,应得结果0.37,本人用int函数求解结果14.44,希望有大神可以帮助求解下!
本人编写程序如下:
clc
close all
clear all
syms w
k=4.9873;
l(1)=0.3066;
l(2)=0.3160;
r1(1)=sqrt(0.0681+(0.9365-w)^2);
r(1)=sqrt(0.0681+(0.6299-w)^2);
r2(1)=sqrt(0.0681+(0.3233-w)^2);
A=30/(sin(k*l(1))*sin(k*l(2)));
f1=sin(k*(l(2)-abs(w)))*sin(k*r1(1))/r1(1);
y(1)=int(f1,'w',-l(2),l(2));
z(1)=vpa(y(1));
f2=sin(k*(l(2)-abs(w)))*sin(k*r2(1))/r2(1);
y(2)=int(f2,'w',-l(2),l(2));
z(2)=vpa(y(2));
f3=(-2)*cos(k*l(1))*sin(k*(l(2)-abs(w)))*sin(k*r(1))/r(1);
y(3)=int(f3,'w',-l(2),l(2));
z(3)=vpa(y(3));
R=A*(z(1)+z(2)+z(3)) |