PAT(乙级) 1054 求平均值

本题的基本要求非常简单:给定 N 个实数,计算它们的平均值。但复杂的是有些输入数据可能是非法的。一个“合法”的输入是 [−1000,1000] 区间内的实数,并且最多精确到小数点后 2 位。当你计算平均值的时候,不能把那些非法的数据算在内。

输入格式:

输入第一行给出正整数 N(≤100)。随后一行给出 N 个实数,数字间以一个空格分隔。

输出格式:

对每个非法输入,在一行中输出 ERROR: X is not a legal number,其中 X 是输入。最后在一行中输出结果:The average of K numbers is Y,其中 K 是合法输入的个数,Y 是它们的平均值,精确到小数点后 2 位。如果平均值无法计算,则用 Undefined 替换 Y。如果 K 为 1,则输出 The average of 1 number is Y

输入样例 1:

7
5 -3.2 aaa 9999 2.3.4 7.123 2.35

输出样例 1:

ERROR: aaa is not a legal number
ERROR: 9999 is not a legal number
ERROR: 2.3.4 is not a legal number
ERROR: 7.123 is not a legal number
The average of 3 numbers is 1.38

输入样例 2:

2
aaa -9999

输出样例 2:

ERROR: aaa is not a legal number
ERROR: -9999 is not a legal number
The average of 0 numbers is Undefined

 

        这题的通过率很低,大概是因为情况太复杂了,这题是比较复杂的字符串处理,要想到各种情况,稍微不注意就出错了,而且这题也有不严谨的地方,末尾有小数点在本题中是被认为合法的,这导致我错了第4个测试点。一些人用了库函数来轻松解决问题,但我还是考虑了所有情况,做了复杂的判断和处理。

        首先说我是用字符串来接收一个实数的,先做判断,再将其转换为数字,判断和转换分别写了两个函数。在判断函数中我们要考虑以下几种情况:

        1、遇到了小数点。如果这个小数点在第一位,或者之前出现过小数点,或者这个小数点离末尾的距离大于3,那么可以说这个数是不合法的。

        2、遇到了负号。如果这个负号不是在第一位,或者出现了多次,或者这个字符串中只包含这一个负号,或者负号后面直接跟了小数点,那么可以说这个数是不合法的。

        3、遇到了除数字、小数点、负号之外的其他字符,那么这个数是不合法的。

        以上这些情况不一定全部有用,可能有一些判断是多余的,但既然想到了这些就都写进程序里了。判断后则将其转换为数字,这个数有可能是带小数部分的,而且整数和小数的转换方法不同,所以要分开处理最后加在一起。得到数之后判断是否在-1000到1000的范围内,然后求和求平均值输出。代码如下,详见注释:

#include<bits/stdc++.h>
using namespace std;
bool judge(string s)
{
	bool flag1=false,flag2=false;//分别标记小数点和负号的出现次数 
	for(int i=0;i<s.length();i++)
	{
		if(s[i]=='.')
		{
			if(flag1==true || i==0)//出现多次和出现在首位的情况 
				return false;
			if(s.length()-i>3)//小数部分超过2位的情况 
				return false;
			flag1=true;
		}
		else if(s[i]=='-')
		{
			if(flag2==true || i!=0 || s.length()==1 || s[i+1]=='.') 
				return false;
			flag2=true;
		}
		else if(s[i]<'0' || s[i]>'9')
			return false;
	}
	return true;
}
float tran(string s)
{
	bool flag=false;
	if(s[0]=='-')//处理负数的情况,先将负号去掉 
	{
		flag=true;
		s=s.substr(1,s.length());
	}
	int pos=s.find(".");
	string s1=s.substr(0,pos),s2="";
	if(pos!=-1)//将小数部分单独截取出来转换 
		s2=s.substr(pos+1,s.length());
	float res=0,res1=0;
	for(int i=0;i<s1.length();i++)//整数部分的转换 
		res=res*10+s1[i]-'0';
	for(int i=s2.length()-1;i>=0;i--)//小数部分的转换 
	{
		res1+=s2[i]-'0';
		res1/=10;
	}
	res+=res1;
	if(flag)
		res=-res;
	return res;
}
int main()
{
	int n,count=0;
	float sum=0;
	cin>>n;
	string s;
	for(int i=0;i<n;i++)
	{
		cin>>s;
		if(!judge(s) || tran(s)>1000 || tran(s)<-1000)
		{
			cout<<"ERROR: "<<s<<" is not a legal number"<<endl;
			continue;
		}
		count++;
		sum+=tran(s);
	}
	if(count==1)
	{
		cout<<"The average of 1 number is ";
		printf("%.2f",sum);
	}
	else if(count>1)
	{
		cout<<"The average of "<<count<<" numbers is ";
		printf("%.2f",sum/(count*1.0));
	}
	else
		cout<<"The average of 0 numbers is Undefined";
	return 0;
}

 

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 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、付费专栏及课程。

余额充值