版块导航: 正在加载中...

登录注册

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

24小时热门版块排行榜

返回列表

616840392

金虫 (小有名气)

应助: 0 (幼儿园)
金币: 1181.1
散金: 280
红花: 2
帖子: 118
在线: 150.3小时
虫号: 2256776
注册: 2013-01-23
性别: GG
专业: 力学

[求助] 求助使用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°）。
请会做这方面的朋友帮忙解决一下，谢谢！

表达式.jpg

回复此楼

» 猜你喜欢

2026年循环经济功能材料国际会议（ICFMCE 2026）已经有0人回复
2026年第五届电气、电子与信息工程国际会议（ISEEIE 2026）已经有0人回复
物理学I论文润色/翻译怎么收费? 已经有76人回复
哈尔滨理工大学物理系招收物理学考研调剂已经有5人回复
0702一志愿吉大B区求调剂已经有5人回复
求调剂已经有0人回复
0702一志愿吉大B区求调剂有论文已经有0人回复
请问还有没有用Latex写文章的小伙伴们？已经有0人回复
光学工程学硕调剂信息已经有29人回复
欢迎加入课题组已经有0人回复
散金币，求好运，祝面上顺利！已经有53人回复

» 本主题相关价值贴推荐，对您同样有帮助:

matlab高手进！发现一个奇怪的现象…… 已经有4人回复
大神！这个求定积分不等式的证明怎么做呢？已经有3人回复
用matlab求解一个隐三角函数方程已经有10人回复
求问，关于定积分和不定积分软件计算已经有3人回复
1/x与正切函数tanx的不定积分求解已经有7人回复
求高手帮忙， MATLAB求积分上限。其中A,E,R为可输入的变值，求hc 已经有5人回复
matlab编程约束条件问题已经有10人回复
数学定积分无法求解，求高手已经有9人回复
MATLAB中，没有数学表达式的话，能不能做积分？已经有3人回复
求助一个关于计算磁阻的定积分已经有4人回复
新手matlab求助，谢谢各位已经有10人回复
求助高手！！！！！！！！！！！matlab最小约束的多变量寻优函数求解方程已经有5人回复
matlab不定积分的问题已经有10人回复
特殊函数求解高手请进，这个类似于Gamma函数的定积分能否求解？已经有12人回复
Matlab数值积分结果是虚数已经有6人回复
求高手解答Matlab求解微分方程的初始值和边界值问题已经有5人回复
10次代数方程怎么找解析解最好详细给出matlab程序谢谢已经有5人回复
求助如何用Matlab求解复杂的不规则积分区间的二重定积分的数值解？已经有3人回复
锂离子电池matlab仿真的问题已经有6人回复
matlab 积分计算算不出结果，求帮助！！！已经有10人回复
matlab 气泡边缘提取求助高手！！已经有5人回复
【求助】【请教】高物试题，请教高手！！邮箱hgaoliang@163.com 已经有6人回复
【求助】请教Matlab高手求定积分近似值问题已经有4人回复
【求助】请教Matlab高手根据积分画图问题已经有4人回复

零五零六零七又零八，长亭折柳，届届应有毕业时；或春或夏或秋亦或冬，故地重游，时时可为归宁日。

1楼 2014-12-17 16:03:28

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

专家经验: +1059
计算强帖: 8
应助: 1712 (讲师)
贵宾: 8.888
金币: 68127.7
散金: 1938
红花: 443
沙发: 4
帖子: 4373
在线: 3291.4小时
虫号: 1122189
注册: 2010-10-14
专业: 宇宙学
管辖: 计算模拟区

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +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')

附图1.png

赞一下(1人)

回复此楼

MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。

2楼2014-12-17 17:09:20