应该这两句有问题，sum=sum+(f1-freal)，RMES=sqrt((sum^2)/(N-1));；可以改成sum=sum+(f1-freal)^2，RMES=sqrt((sum)/(N-2));运行结果为RMES = 0.8385
