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},但我不知道错误是什么。
感谢这会让我终于完成我的项目 – RHLS