查看: 171|回复: 4|关注: 0

[已解决] 求助积分与画曲线问题

[复制链接]

新手

12 麦片

财富积分


050


3

主题

6

帖子

0

最佳答案

ym=135.6;rk=511.1;rd=502.7;rm=403.5;re=333.1;rb=271.0;%已知
t0=0.8;te=0.15;a=0.1;%已知
b=rd^2-rm^2+(rk^2-rd^2)*(1-t0+a*rk/(rk-rd))-2*a*(rk^3-rd^3)/(3*(rk-rd));%常数
syms r1;
g1=(rd^2-rm^2+(1-t0+a*rk/(rk-rd))*(r1^2-rd^2)-2*a*(r1^3-rd^3)/(3*(rk-rd)))/(b^2-(rd^2-rm^2+(1-t0+a*rk/(rk-rd))*(r1^2-rd^2)-2*a*(r1^3-rd^3)/(3*(rk-rd)))^2)^1/2;
y1=int(g1,r1);%求g1的积分
问题是积分结果无表达式也无错误提示。图也没法画。原公式见图片。还要求用该积分公式求出yd,yd是该积分曲线上纵坐标rd对应的横坐标值。然后画出二维曲线图(区间:rd~rk)谢谢!

回复主题 已获打赏: 0 积分

举报

论坛优秀回答者

9

主题

1629

帖子

342

最佳答案
  • 关注者: 81
发表于 2020-9-2 09:46:25 | 显示全部楼层
没有看到图片,看现有代码是要求不定积分,不定积分有表达式,在命令窗口输入y1
或将
  1. y1=int(g1,r1);%求g1的积分
复制代码

改成
  1. y1=int(g1,r1)%求g1的积分
复制代码
回复此楼 已获打赏: 0 积分

举报

新手

12 麦片

财富积分


050


3

主题

6

帖子

0

最佳答案
 楼主| 发表于 2020-9-2 10:45:46 | 显示全部楼层
好的,谢谢!表达式很怪异,也没法画曲线啊!图片大小也在4M以下,但就是没法传上去。公式 公式.png
回复此楼 已获打赏: 0 积分

举报

论坛优秀回答者

入门

203 麦片

财富积分


50500


0

主题

471

帖子

41

最佳答案
  • 关注者: 12
发表于 2020-9-7 23:18:40 | 显示全部楼层 |此回复为最佳答案
解答如下
clear all
clc
ym=135.6;rk=511.1;rd=502.7;rm=403.5;re=333.1;rb=271.0;%已知
t0=0.8;te=0.15;a=0.1;%已知
b=rd^2-rm^2+(rk^2-rd^2)*(1-t0+a*rk/(rk-rd))-2*a*(rk^3-rd^3)/(3*(rk-rd))%常数
r=linspace(rd,rk,200);
for i=1:length(r)
g1=@(r1) (rd.^2-rm.^2+(1-t0+a.*rk./(rk-rd)).*(r1.^2-rd.^2)-2.*a.*(r1.^3-rd.^3)./(3.*(rk-rd)))./(b.^2-(rd.^2-rm.^2+(1-t0+a.*rk./(rk-rd)).*(r1.^2-rd.^2)-2.*a.*(r1.^3-rd.^3)./(3.*(rk-rd))).^2).^(1/2);
y1(i)=integral(g1,r(i),rk);%求g1的积分
end
plot(r,y1)
回复此楼 已获打赏: 2 积分

举报

新手

12 麦片

财富积分


050


3

主题

6

帖子

0

最佳答案
 楼主| 发表于 4 天前 | 显示全部楼层
wjb986555360 发表于 2020-9-7 23:18
解答如下
clear all
clc

谢谢您的解答!我还有一个类似的问题,想请大牛帮助解决下,我做别的专业的,零时用到了,请帮帮忙啊!
回复此楼 已获打赏: 0 积分

举报

您需要登录后才可以回帖 登录 | 注册

本版积分规则

关闭

站长推荐上一条 /4 下一条

快速回复 返回顶部 返回列表