这是一道老师出的matlab实验题，可以麻烦您看看告诉我怎么做吗？

求解方法：

1、自定义你的积分方程函数。主要内容

func=@(t,u)(1-(1-qfunc(u))^(M-1))exp(-(u-sqrt(t))^2/2)exp(-t/x)/(xsqrt(2pi));

y=integral2(func,0,+inf,-inf,+inf);

2、利用循环语句求解XdB与Y的关系值。

3、利用semilogy（）函数绘出M=2、4、8、16情况下的XdB与Y的曲线图

semilogy(XdB,YdB(1,:),XdB,YdB(2,:),XdB,YdB(3,:),XdB,YdB(4,:))

Function FuncMin(arr() As Integer) As Integer

FuncMin = arr(LBound(arr))

For i = LBound(arr) + 1 To UBound(arr)

If arr(i) < FuncMin Then

FuncMin = arr(i)

End If

Next

End Function

test调用：

Private Sub Command1_Click()

Dim arr(1 To 3) As Integer

arr(1) = 2

arr(2) = 1

arr(3) = 4

Print FuncMin(arr)

End Sub

您好，un为目标函数，它可用前面的方法定义；

x0为初始值；

A、b满足线性不等式约束 ，若没有不等式约束，则取A=[ ]，b=[ ]；

Aeq、beq满足等式约束 ，若没有，则取Aeq=[ ]，beq=[ ]；

lb、ub满足 ，若没有界，可设lb=[ ]，ub=[ ]；

nonlcon的作用是通过接受的向量x来计算非线性不等约束 和等式约束 分别在x处的估计C和Ceq，通过指定函数柄来使用，如：>>x = fmincon(@myfun,x0,A,b,Aeq,beq,lb,ub,@mycon)，先建立非线性约束函数，并保存为myconm：function [C,Ceq] = mycon(x)

C = …

% 计算x处的非线性不等约束 的函数值。

Ceq = …

% 计算x处的非线性等式约束 的函数值。

lambda是Lagrange乘子，它体现哪一个约束有效。

output输出优化信息；

grad表示目标函数在x处的梯度；

hessian表示目标函数在x处的Hessiab值。

注意：

1 fmincon 函数提供了大型优化算法和中型优化算法。默认时，若在 fun 函数中提供了梯度（options 参数的 GeadObj 设置为 'on'），并且只有上下界存在或只有等式约束，fmincon 函数将选择大型算法。 当既有等式约束又有梯度约束时，使用中型算法。

2 fmincon 函数的中型算法一般是使用序列二次规划。在每一步迭代中求解二次规划子问题，并用 BFGS 法更新 Lagrangian 乘子和 Hessian 矩阵。

3 fmincon 函数的大型算法采用了subspace trust region 优化算法。这种算法是把目标函数在点x的邻域泰勒展开（x可以认为是人为提供的初始猜测），这个展开的邻域就是所谓的trust region，泰勒展开进行到二阶项为止。

4 fmincon 函数可能会给出局部最优解，这与初始值的选取有关。