Question

MATLAB程序中有一组y的值作积分下限，怎么来表示

H=26.15;L=263.75;
x=0:5:130;
X=263.75/2-x;
y=2.8*(X+L*(0.015+0.035*log(H/L)))/(0.5*L-L*(0.015+0.035*log(H/L)));
syms x,y
f=1-1/sqrt(pi)*vpa(int (exp(-x^2), y,inf))

Answer 1

在MATLAB程序中有一组y的值作积分下限，可以这样来处理。即题主在原代码中调整一些代码的位置就可以解决。

具体代码如下：

H=26.15;L=263.75;

syms x

X=263.75/2-x;

y=2.8*(X+L*(0.015+0.035*log(H/L)))/(0.5*L-L*(0.015+0.035*log(H/L)));

f=1-1/sqrt(pi)*vpa(int (exp(-x^2), y,inf))

x=0:5:130;

fx=eval(f);

plot(x,fx,'r.')

xlabel('x'),ylabel('f(x)');

执行上述代码，可以得到  f(x) 的曲线图。

追问

谢谢，我现在不要图像，要得到具体的f的值应该怎么编写程序呢

追答

你运行上述代码，可以得到
f(x)=0.50000000000000000679534861169047 - 0.49999999999999999320465138830953*erf(0.018760019467255907532666527325491*x - 2.1479551344887456117407965820984)