matlab中知道参数如何写带积分的方程

Question

假如w=[0.0529;0.3981;0.0892;0.4596]
y=[58.80;2.45;32.71;1.01]
z=[-3.15;2.78;-0.29;-0.60]
如何按下面的公式，写一个随时代入x的方程，还有一个把原公式对x不定积分后，再代入x的方程。x是一组弧度角

Answer 1

随时代入x是比较方便的

clc; clear
w=[0.0529;0.3981;0.0892;0.4596];
y=[58.80;2.45;32.71;1.01];
z=[-3.15;2.78;-0.29;-0.60];
beta(:,1) = y;
beta(:,2) = z;
beta(:,3) = w;
f = @(x)fun(x, beta);
% 这里给定x
x = [1,2,3];
f(x)

function f = fun(x, yzw)
n = numel(yzw);
m = n/3;
y = yzw(1:m);
z = yzw(m+1:2*m);
w = yzw(2*m+1:3*m);
I0 = @(y) 1/(2*pi)*integral(@(x)exp(y.*cos(x)),0,2*pi);
f = zeros(size(x));
for i = 1:m
    f = f + w(i)*1/(2*pi*I0(y(i)))*exp(y(i)*cos(x-z(i)));
end
end

不定积分做符号运算是有难度的