多项式操作的c++实现

#include"define.h"

typedef struct
{
 float coef;//系数
 int expn;//指数
}term,ElemType;
typedef struct LNode
{
    ElemType data;
    struct LNode *next;
}*Link;
typedef struct
{
    Link head,tail;
    int len;
}LinkList;

typedef LinkList polynomial;

status MakeNode(Link &p,ElemType e)
{
    //分配由p指向的值为e的节点
    p=(Link)malloc(sizeof(LNode));
    if(p)
    {
        p->data.coef=e.coef;
  p->data.expn=e.expn;
        p->next=NULL;
        return OK;
    }
    else
        return 0;
}

Link GetHead(LinkList L)
{
    return L.head;
}

status SetCurElem(Link &p,ElemType e)
{
    //用e更新p所指节点
 if(!p)
 {
  printf("p is NULL!");
  return ERROR;
 }
    p->data.coef=e.coef;
 p->data.expn=e.expn;
 return OK;
}

Link LocateElem(LinkList L,ElemType e,status(* compare)(ElemType,ElemType))
{
    //返回线性链表中第一个与e满足compare()判定关系的位置;没有则返回NULL
    Link p;
    p=L.head->next;
    while(p!=NULL&&!(*compare)(p->data,e))
        p=p->next;
    if(p!=NULL)
        return p;
    else
        return NULL;
}

status InitList(LinkList &L)
{
    L.head=L.tail=(Link)malloc(sizeof(LNode));
    if(!L.head)
        exit(-2);
    L.len=0;//长度不算头结点
    L.head->next=NULL;
    return OK;
}

status ecqual(ElemType a,ElemType b)
{
 if(a.expn==b.expn)
  return OK;
 return ERROR;
}

void InsFirst(Link h,Link s)
{
    //插入节点(头插)
    s->next=h->next;
    h->next=s;
}

void GetCurElem(Link p,ElemType &a)
{
    //p所直节点的值赋给a

  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
#ifndef FUNCTION_H_ #define FUNCTION_H_ #include #include #include "polyfit.h" #include using namespace std; dxs::dxs() { ifstream fin("多项式拟合.txt"); fin>>n; x=new float[n]; y=new float[n]; for(int i=0;i>x[i]; } for(i=0;i>y[i]; } cout<>nn; m=nn+1; u=new float*[m]; for(i=0;i<m;i++) { u[i]=new float[m+1]; }//创建m行,m+1列数组 } void dxs::dfine() { for(int i=0;i<m;i++) { for(int j=0;j<m+1;j++) { u[i][j]=0; } } for(i=0;i<m;i++) { for(int j=0;j<m;j++) { for(int k=0;k<n;k++) { u[i][j]=u[i][j]+pow(x[k],j+i); } } } for(i=0;i<m;i++) { for(int k=0;k<n;k++) { u[i][m]=u[i][m]+pow(x[k],i)*y[k]; } } } void dxs::show() { for(int i=0;i<m;i++) { for(int j=0;j<m+1;j++) { cout<<u[i][j]<<" ";//<<endl; } cout<<endl; } ////显示具有m行m+1列u数组的各元素值 } void dxs::select_main(int k,float **p,int m) { double d; d=*(*(p+k)+k); //cout<<d; int l=k; int i=k+1; for(;i fabs(d)) { d=*(*(p+i)+k); l=i; } else continue; } if(d==0) cout<<"错误"; else { if(k!=l) { for(int j=k;j<m+1;j++) { double t; t=*(*(p+l)+j); *(*(p+l)+j)=*(*(p+k)+j); *(*(p+k)+j)=t; } } } } void dxs::gaosi() { for(int k=0;k<m;k++) { select_main(k,u,m);//调用列主元函数 for(int i=1+k;i<m;i++) { // *(*(p+i)+k)=(float) *(*(p+i)+k) / *(*(p+k)+k); u[i][k]=(float) u[i][k] / u[k][k]; } for(i=k+1;i<m;i++) { for(int j=k+1;j=0;i--) { float a=0; for(int j=i+1;j<m;j++) { //a=a + (*(*(p+i)+j) * *(*(p+j)+m)); a=a+u[i][j] * u[j][m]; } //*(*(p+i)+n-1)= (*(*(p+i)+n-1) - a) / *(*(p+i)+i); u[i][m]= (u[i][m] -a) / u[i][i]; } cout<<"方程组的解为:"<<endl; for(i=0;i<m;i++) { cout<<"a"<<i+1<<"="; cout<<u[i][m]<<endl; // l[i]=*(*(p+i)+n-1); } cout<<"y="<<u[0][m]; for(i=1;i<m;i++) { cout<<showpos<<u[i][m]<<"x"; if(i!=1)cout<<"^"<<noshowpos<<i; } cout<<endl; } dxs::~dxs() { delete[]x,y; delete []*u; } #endif
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值