C语言一元多项式的乘积与和(自己复习使用)

设计函数分别求两个一元多项式的乘积与和。

输入格式:

输入分2行,每行分别先给出多项式非零项的个数,再以指数递降方式输入一个多项式非零项系数和指数(绝对值均为不超过1000的整数)。数字间以空格分隔。

输出格式:

输出分2行,分别以指数递降方式输出乘积多项式以及和多项式非零项的系数和指数。数字间以空格分隔,但结尾不能有多余空格。零多项式应输出0 0

样例">输入样例:

4 3 4 -5 2  6 1  -2 0
3 5 20  -7 4  3 1

输出样例:

15 24 -25 22 30 21 -10 20 -21 8 35 6 -33 5 14 4 -15 3 18 2 -6 1
5 20 -4 4 -5 2 9 1 -2 0
#include<stdio.h>
#include<malloc.h>
typedef struct PolyNode{
	int coef;//系数 
	int expon;//指数 
	struct PolyNode *link;//指向下一个节点的指针 
}*Polynomial;
Polynomial ReadPoly();
Polynomial Mult(Polynomial P1,Polynomial P2);
Polynomial Add(Polynomial P1,Polynomial P2);
void PrintPoly(Polynomial P);
void Attach(int c,int e,Polynomial *pRear);
int Compared(int a,int b);
int main()
{
	Polynomial P1,P2,PP,PS;
	P1=ReadPoly();
	P2=ReadPoly();
    PP=Mult(P1,P2);
	PrintPoly(PP);
    PS=Add(P1,P2);
    PrintPoly(PS);
}
Polynomial ReadPoly()//创建多项式,思想:创建一个无用头结点,然后通过attach函数一个一个的把输入进来的系数c和指数e放到节点并连接起来, 
{
	Polynomial P,Rear,t;
	int c,e,N;
	scanf("%d",&N);//多项式的项数
	P=(Polynomial)malloc(sizeof(struct PolyNode));//构造一个临时不存数据的头节点
	P->link=NULL;
	Rear=P; 
	while(N--)
	{
		scanf("%d %d",&c,&e);
		Attach(c,e,&Rear);
	 } 
	 t=P;
	 P=P->link;
	 free(t);//把临时建立的头结点去掉,然后返还真正头节点的指针 
	 return P;
}
void Attach(int c,int e,Polynomial *pRear)//将多项式的每一项用链表以此链接起来 ,实现思路:把传进来的系数指数放到一个新的结点,然后将传进来的rear(即最后一个结点的指针)指向新建的结点,再把rear再次指向最后一个结点 
{
	Polynomial P;
	P=(Polynomial)malloc(sizeof(struct PolyNode));
	P->coef=c;
	P->expon=e;
	P->link=NULL;//记得要给后一个结点指向NULL 
	(*pRear)->link=P;//pRear是指针的指针,加个*就是指针rear了 
	*pRear=P;
}
Polynomial Add(Polynomial P1,Polynomial P2)//将两个多项式按次幂顺序加起来 思路:构造一个无用头结点,然后依次对比两个多项式的幂次,根据相减结果大小将项按次幂顺序插到头结点后面,如果一个插完了另一个还有多,就将剩余的全部接到后面,最后还要将头结点指向真的第一个结点,家的free掉 
{
	Polynomial front,rear,temp;
	int sum;
	rear=(Polynomial)malloc(sizeof(struct PolyNode));//先构造一个没用的头结点用于链接多项式 
	rear->link=NULL;
	front=rear;//front记录多项式的链表头节点 
	while(P1&&P2)
	     switch(Compared(P1->expon,P2->expon)){
	     	case 1:
	     		Attach(P1->coef,P1->expon,&rear);
	     	     P1=P1->link;
	     		break;
	     	case -1:
	     		Attach(P2->coef,P2->expon,&rear);
			    P2=P2->link;
			    break;
	     	case 0:
	     		sum=P1->coef+P2->coef;
			    if(sum)
			      Attach(sum,P1->expon,&rear);
			    P1=P1->link;
		        P2=P2->link;
			    break; 
	     		}
	for(;P1;P1=P1->link)Attach(P1->coef,P1->expon,&rear);
	for(;P2;P2=P2->link)Attach(P2->coef,P2->expon,&rear);
		rear->link=NULL;
		temp=front;
		front=front->link;//令front指向结果多项式的第一个非0项 
		free(temp);//释放临时的空表头节点 
		return front;
	} //讨论了三种情况的对应方法 
void PrintPoly(Polynomial P)//输出多项式 
{
	int flag=0;//辅助调整输出格式用 ,两项之间的空格就通过flag添加 
	if(!P)
	{
		printf("0 0\n");
		return;
	}
	while(P)
	{
		if(!flag)
		  flag=1;
		else
		  printf(" ");
		printf("%d %d",P->coef,P->expon);
		P=P->link;
	}
	printf("\n");
	return;
}
int Compared(int a,int b)//在add里使用的一个小函数 
{
	if(a>b)
	return 1;
    if(a<b)
    return -1;
    if(a==b)
    return 0;
}
Polynomial Mult(Polynomial P1,Polynomial P2)//思路:先通过P1第一项乘P2所有项得到一个有序的多项式,再通过两个循环扫描把剩余的相乘会得到的项数以此插入,保证有序 
{
	Polynomial t1,t2,P,Rear,t;
	int c,e;
	if(!P1||!P2)  return NULL;//只要有一个为空则返回空 
	t1=P1;t2=P2;
	P=(Polynomial)malloc(sizeof(struct PolyNode));P->link=NULL;
	Rear=P;
	while(t2){ //先用P1的第一项乘以P2,得到P ,即一个有序的多项式链 
		Attach(t1->coef*t2->coef,t1->expon+t2->expon,&Rear);
		t2=t2->link;
	}
	t1=t1->link;
	while(t1)//接下来就是在一个已知有序的多项式中插入其他的各个项,并且保证有序 
	{
		t2=P2;Rear=P;
		while(t2)
		{
			e=t1->expon+t2->expon;
			c=t1->coef*t2->coef;
			while(Rear->link&&Rear->link->expon>e)//扫描,如果扫描发现一个的项的幂次比自己大,那就rear往后指一个继续扫描,直到出现后一个比要插入的小或等  
			   Rear=Rear->link;
			if(Rear->link&&Rear->link->expon==e)
			{
				if(Rear->link->coef+c)//如果系数相加不等于0那就合并
				Rear->link->coef+=c;
				else{//如果系数相加等于0那就要把这个节点删掉
				t=Rear->link; 
				Rear->link=t->link;
				free(t);	
				} 
			 } 
			 else{
			 	t=(Polynomial)malloc(sizeof(struct PolyNode));//创建一个新的节点,把他插入 
			 	t->coef=c;t->expon=e;
			 	t->link=Rear->link;
			 	Rear->link=t;
			 	Rear=Rear->link;//Rear往后挪 
			 }
			t2=t2->link ;
		}
		t1=t1->link;
	}
	t2=P;P=P->link;free(t2);
	return P;
}

	
  • 1
    点赞
  • 3
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
在信号处理领域,DOA(Direction of Arrival)估计是一项关键技术,主要用于确定多个信号源到达接收阵列的方向。本文将详细探讨三种ESPRIT(Estimation of Signal Parameters via Rotational Invariance Techniques)算法在DOA估计中的实现,以及它们在MATLAB环境中的具体应用。 ESPRIT算法是由Paul Kailath等人于1986年提出的,其核心思想是利用阵列数据的旋转不变性来估计信号源的角度。这种算法相比传统的 MUSIC(Multiple Signal Classification)算法具有较低的计算复杂度,且无需进特征值解,因此在实际应用中颇具优势。 1. 普通ESPRIT算法 普通ESPRIT算法两个主要步骤:构造等效旋转不变系统和估计角度。通过空间平移(如延时)构建两个子阵列,使得它们之间的关系具有旋转不变性。然后,通过对子阵列数据进最小二乘拟合,可以得到信号源的角频率估计,进一步转换为DOA估计。 2. 常规ESPRIT算法实现 在描述中提到的`common_esprit_method1.m`和`common_esprit_method2.m`是两种不同的普通ESPRIT算法实现。它们可能在实现细节上略有差异,比如选择子阵列的方式、参数估计的策略等。MATLAB代码通常会包含预处理步骤(如数据归一化)、子阵列构造、旋转不变性矩阵的建立、最小二乘估计等部。通过运两个文件,可以比较它们在估计精度和计算效率上的异同。 3. TLS_ESPRIT算法 TLS(Total Least Squares)ESPRIT是对普通ESPRIT的优化,它考虑了数据噪声的影响,提高了估计的稳健性。在TLS_ESPRIT算法中,不假设数据噪声是高斯白噪声,而是采用总最小二乘准则来拟合数据。这使得算法在噪声环境下表现更优。`TLS_esprit.m`文件应该包含了TLS_ESPRIT算法的完整实现,包括TLS估计的步骤和旋转不变性矩阵的改进处理。 在实际应用中,选择合适的ESPRIT变体取决于系统条件,例如噪声水平、信号质量以及计算资源。通过MATLAB实现,研究者和工程师可以方便地比较不同算法的效果,并根据需要进调整和优化。同时,这些代码也为教学和学习DOA估计提供了一个直观的平台,有助于深入理解ESPRIT算法的工作原理。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值