应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 求助使用Matlab求解定积分的问题，望高手帮忙解答。
21/1返回列表
查看: 533  |  回复: 1
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

616840392

金虫 (小有名气)

[求助] 求助使用Matlab求解定积分的问题，望高手帮忙解答。 已有1人参与

所求积分的表达式如图所示：
主要目的是求出表达式右边的积分，其中alpha1和alpha是未知数。想通过积分求得关于alpha1的表达式然后代入alpha1的取值范围画出相应的图形。
目前写的求解关于alpha1的表达式的程序如下：
（其中t代表alpha，x代表alpha1）
clear;clc;
u = 0.3;
E = 2.068e11;
L = 0.016;
alpha2 = 6.2819*pi/180;
syms t x
f=1.2.*(1+u).*(alpha2-t).*cos(t).*(cos(x)).^2/(E*L*(sin(t)+(alpha2-t)*cos(t)));
f_s=int(f,t,-x,alpha2)
上述程序无法得到想要的表达式，最后的结果如下：
f_s =int(-(cos(t)*cos(x)^2*((39*t)/25 - 308114935332307701/1801439850948198400))/(3308800000*sin(t) - 3308800000*cos(t)*(t - 7900382957238659/72057594037927936)), t, -x, 7900382957238659/72057594037927936)
由此说明上述方法存在问题，请高手帮忙解决一下。另外希望可以告诉一下如何根据求得的关于alpha1的表达式f_s画出图形（以f_s为纵坐标，alpha1为横坐标，一直alpha1的取值范围为：（18.9965°，24.0475°）。
请会做这方面的朋友帮忙解决一下，谢谢！

求助使用Matlab求解定积分的问题，望高手帮忙解答。
表达式.jpg
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
零五零六零七又零八，长亭折柳，届届应有毕业时；或春或夏或秋亦或冬，故地重游，时时可为归宁日。
1楼 2014-12-17 16:03:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
616840392: 金币+40, ★★★★★最佳答案, 非常感谢！以后相关问题可能还需要多多向你请教！谢谢！ 2014-12-18 10:16:33
CODE:
clear all;clc;
u = 0.3;
E = 2.068e11;
L = 0.016;
alpha2 = 6.2819*pi/180;
alpha1=[18.9965:0.1:24.0475]/180*pi;
for i=1:length(alpha1)
    x=alpha1(i);
    f(i)=quadgk(@(t)1.2.*(1+u).*(alpha2-t).*cos(t).*(cos(x)).^2./(E*L*(sin(t)+(alpha2-t).*cos(t))),-x,alpha2  );
end

figure(1)
plot(alpha1,f,'b-o')

求助使用Matlab求解定积分的问题，望高手帮忙解答。-1
附图1.png
赞一下(1人) 回复此楼
MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。
2楼2014-12-17 17:09:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 616840392 的主题更新
21/1返回列表
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭