用python写一个小程序，输入坐标求线性回归

你好：

上面的程序，请看如下代码：

# -*- coding: cp936 -*-

end=input("是否结束(y/n)：")

while end=="n":

    print "Number of coordinates：2"

    xx=input("x's：")

    yy=input("y's：")

    a=float(list(xx)[0])

    b=float(list(xx)[1])

    c=float(list(yy)[0])

    d=float(list(yy)[1])

    print "第一个点是：("+str(a)+","+str(c)+")"

    print "第一个点是：("+str(b)+","+str(d)+")"

    x0=c-a

    y0=float(d-b)

    print "直线方程为：",

    if x0==0:

        print "x=",a

    else:

        print "y=%r(x-%r)+%r"%(y0/x0,a,c)

    print

#include<stdio.h>

#include<stdlib.h>

#include<iostream.h>

typedef struct data

{

float x

float y

}Data//变量x和函数值y的结构

Data d[20]//最多二十组数据

float f(int s,int t)//牛顿插值法，用以返回插商

{

if(t==s+1)

return (d[t].y-d[s].y)/(d[t].x-d[s].x)

else

return (f(s+1,t)-f(s,t-1))/(d[t].x-d[s].x)

}

float Newton(float x,int count)

{

int n

while(1)

{

cout<<"请输入n值(即n次插值):"//获得插值次数

cin>>n

if(n<=count-1)// 插值次数不得大于count-1次

break

else

system("cls")

}

//初始化t，y，yt。

float t=1.0

float y=d[0].y

float yt=0.0

//计算y值

for(int j=1j<=nj++)

{

t=(x-d[j-1].x)*t

yt=f(0,j)*t

//cout<<f(0,j)<<endl

y=y+yt

}

return y

}

float lagrange(float x,int count)

{

float y=0.0

for(int k=0k<countk++)//这儿默认为count-1次插值

{

float p=1.0//初始化p

for(int j=0j<countj++)

{//计算p的值

if(k==j)continue//判断是否为同一个数

p=p*(x-d[j].x)/(d[k].x-d[j].x)

}

y=y+p*d[k].y//求和

}

return y//返回y的值

}

void main()

{

float x,y

int count

while(1)

{

cout<<"请输入x[i],y[i]的组数，不得超过20组:"//要求用户输入数据组数

cin>>count

if(count<=20)

break//检查输入的是否合法

system("cls")

}

//获得各组数据

for(int i=0i<counti++)

{

cout<<"请输入第"<<i+1<<"组x的值:"

cin>>d[i].x

cout<<"请输入第"<<i+1<<"组y的值:"

cin>>d[i].y

system("cls")

}

cout<<"请输入x的值:"//获得变量x的值

cin>>x

while(1)

{

int choice=3

cout<<"请您选择使用哪种插值法计算:"<<endl

cout<<" (0):退出"<<endl

cout<<" (1):Lagrange"<<endl

cout<<" (2):Newton"<<endl

cout<<"输入你的选择:"

cin>>choice//取得用户的选择项

if(choice==2)

{

cout<<"你选择了牛顿插值计算方法，其结果为:"

y=Newton(x,count)break//调用相应的处理函数

}

if(choice==1)

{

cout<<"你选择了拉格朗日插值计算方法，其结果为:"

y=lagrange(x,count)break//调用相应的处理函数

}

if(choice==0)

break

system("cls")

cout<<"输入错误!!!!"<<endl

}

cout<<x<<" , "<<y<<endl//输出最终结果

}

假设你的实际实在A2单元格，则可在B2单元格输入公式：

=X%*(60+MAX(0,MIN(40,(100-60)*(A2-80)/(150-80))))

化简后公式为：

=X%*(60+MAX(0,MIN(40,(A2-80)*4/7)))

---