2015-04-24 34 views
0

我尝试在GNU Octave上编写线性&非线性回归,作为我的Numerical Method类的最终项目的一部分。不过,我的代码中有一个错误,这是一个简单的解析错误,但我不知道如何解决它。我尝试修正语法或尝试改变我的编码方式仍然卡住。需要帮助解析错误GNU Octave

printf("This program is use for curve fitting\n") 
printf("Please enter the data in the [] by use one space to seperate each element.\n") 
x = input("Enter the data(x group):") 
y = input("Enter the data(y group):") 
mode = input("What kind of fitting you want to do(expo((1),power(2),least-square fitting(3)):") 
switch mode 
case{3} 
    xbar = sum(x)/length(x) 
    ybar = sum(y)/length(y) 
    xy = x.*y 
    xs = x.^2 
    a1 = ((length(x)*sum(xy))-(sum(x)*sum(y)))/((sum(xs)*length(x))-(sum(x)^2)) 
    a0 = ybar-(a1*xbar) 
    disp(x) 
    disp(y) 
    printf("mean of x = %d,mean of y = %d,a1 = %d.a0 = %d",xbar,ybar,a1,a0) 
    plot(a0+(a1*x)) 
case{2} 
    xbar = sum(x)/length(x) 
    ybar = sum(y)/length(y) 
    xy = x.*y 
    xs = x.^2 
    a1 = ((length(x)*sum(xy))-(sum(x)*sum(y)))/((sum(xs)*length(x))-(sum(x)^2)) 
    a0 = ybar-(a1*xbar) 
    alpha = 10^a0 
    disp(x) 
    disp(y) 
    printf("mean of x = %d,mean of y = %d,a1 = %d.a0 = %d,alpha = %d,beta = %d\n",xbar,ybar,a1,a0,alpha,a0) 
    plot((alpha.*(x^a1)) 
case{1} 
    xbar = sum(x)/length(x) 
    ybar = sum(y)/length(y) 
    xy= x.*y 
    xs= x.^2 
    a1 = ((length(x)*sum(xy))-(sum(x)*sum(y)))/((sum(xs)*length(x))-(sum(x)^2)) 
    a0 = ybar-(a1*xbar) 
    disp(x) 
    disp(y) 
    alpha = e^a1 
    printf("mean of x = %d,mean of y = %d,a1 = %d.a0 = %d,alpha = %d,beta = %d\n",xbar,ybar,a1,a0,alpha,a0) 
    plot(a0*(e^((a1*alpha).*x)) 
otherwise 
    break; 
end 

解析错误一直告诉发生错误的地方是{1},但我不知道错误是什么。

回答

1

这种情况下,标签的前行不匹配的括号:

plot((alpha.*(x^a1)) 
+0

感谢这会让我终于完成我的项目 – RHLS