Question

所以，我正在尝试编写一个程序，它在给定的区间[a，b]中计算拉格朗日插值。我的代码工作（我在程序的不同点打印出结果），但我似乎无法得到一个好的情节。这是我的代码：

funcprot(0)
clear all

clf()

    a = -5 
    b = 5 
    j = 100 
    n = [3;9;15;36] 

function fx = func(x) 

    fx = -x/(2+x^6)

endfunction

function y = Lagrange(t, f, x) 

    n_t = size(t,'r') 
    n_x = size(x,'r')
    disp(n_t)
    disp(n_x)
    y = zeros(n_x,1)

    for k = 1:n_x 

    POL = 0

    for i = 1:n_t

    L_k = 1

    for l = 1:n_t

    if (l ~= i) then 

    L_k = L_k * (x(k)-t(l))/(t(i)-t(l))

    end 

    end

    POL =  POL + f(i)*L_k

    end 

    y(k) = POL

    end

endfunction

    x = linspace(a,b,j)' 
    fx = func(x)
    m = size(n, 'r')



    for i=1:m

    xi = linspace(a,b,n(i))'
    fxi = func(xi)

    yi1 = Lagrange(xi, fxi, x) 

    subplot(2,m,i)
    plot2d(x,[yi1,fx], rect=[a,-0.4,b, 0.4], style=[2,5])
    plot2d(xi,fxi, style=[-6])
    title(strcat(['Number of grid-points =: ',string(n(i))]))

    c = zeros(n(i),1)

    for j=1:n(i) 

    c(j) = (a+b)/2+((b-a)/2)*cos(((2*j-1)*%pi)/(2*n(i)))

    end

    fc = func(c)

    yi2 = Lagrange(c,fc,x) 
    subplot(2,m,i+m)
    plot2d(x,[yi2,fx],rect=[a,-0.4,b, 0.4], style=[2,5])
    plot2d(c,fc, style=-6)
    title(strcat(['Chebyshev = ',string(n(i))]))

    end

我总是收到以下错误消息：

plot2d: Wrong size for input argument #7: 2 < 101 expected.

这意味着，我的情节函数plot2d在某种程度上是错误的，但我无法理解为什么。我试过改变＆＃39; rect＆＃39;和＆＃39; style＆＃39;，但这没有帮助，而且，他们只确定曲线的边界和形状。我开始相信我的错误在其他地方，但我似乎无法找到。我什么时候删除＆＃39; rect＆＃39;和＆＃39; style＆＃39;，我得到了一些图表，但它们根本不正确，因为它们只是直线而不是Lagrange-Interpolation。如果有人可以指出我的错误，我将非常感激。干杯！

Answer 1

这里有一些提示让它运行：尝试匹配plot2D描述中请求的顺序和数量的参数：

plot2d(x,[yi1,fx], style=[2,5],,, rect=[a,-0.4,b, 0.4],)
plot2d(xi,fxi, style=[-6],,,)
...
plot2d(x,[yi2,fx], style=[2,5],,,rect=[a,-0.4,b, 0.4],)
plot2d(c,fc, style=-6,,,)

我不能理解为什么你必须这样做，因为它们是选择性的论点，但至少，你有你的结果（或者看起来你已经向前迈进了一步，至少！）

Lagrange-Interpolation-由于错误信息而没有显示图表，我无法修复

1 个答案: