Question

如何使用MATLAB得出这个积分的函数图像？

式子比较复杂，就使用手写代替了，作图的代码如何写呢？

Answer 1

如何使用MATLAB得出题主给出的积分函数图像？

第一步：将n分成若干份，即 n=0:1:100

第二步：利用for循环语句，求解当n等于某值时的y（n）

第三步：自定义被积函数，即

fun=@(q) 表达式

第四步：利用 integral数值积分，求解y（n）值

第五步：利用 plot函数，绘制实数部分，虚数部分的图形，即

plot(n1,y1,'r.-',n1,imag(y1),'k--')

完善代码运行可以得到如下结果。

追问

我运行的结果和您的不一致，不知道什么地方出现了问题？希望您给予指正，谢谢您

追答

看不清，你提供的图片，请直接复制给出。

追问

n1=0:100;
for i=1:length(n1);
n=n1(i);
fun=@(q)q.^2.*((1.6e-19).^(2-q*9.41e-10).*sinh(q*7.5e-11)./(8.16e-11.*(1-(1.6e-19).^(q*5.21e-10)).*(q+7.8538e8))).^2./(sqrt(1-q.^2/(25.12*n)));
n2=5.012*n;
y=integral(fun,0,n2);
y1(i)=y;
end

追答

fun自定义函数，与我的有出入。
fun=@(q)q.^2.*((1.6e-19).^(2-q*9.41e-10).*sinh(q*7.5e-11)./(8.16e-11*(1-1.6e-19).^(-q*5.12e-10)).*(q+7.8538e8)).^2./sqrt(1-q.^2/(25.12*n));

追问

嗯嗯，好，谢谢您，我自己再核对一下