1.#QNAN0,gcc vc6.0 的区别

hlyces 2010-01-21 11:34:24
gcc-4.4.2-32

下:[2240] term:return value:1.#QNAN0

vc 6.0

下:[2244] term:return value:1.000000
/* Copyright(C)
* For free
* All right reserved
*
*/
/**
* @file main.cpp
* @brief calculator
* @author hlyces, wsklrw@163.com
* @version 1.0
* @date 2010-01-16
*/
#include "std_lib_facilities.h"

#include <windows.h>

class Token{
public:
char kind;
double value;
string name;
Token(char ch)
:kind(ch),value(0){}
Token(char ch, double val)
:kind(ch),value(val){}
Token( char ch, string n): kind( ch), name( n){}
};

class Token_stream
{
public:
Token_stream( ):full( false),buffer( 0){

}
Token get( );
void putback( Token t);
void ignore( char c);
private:
bool full;
Token buffer;
};

class Variable{
public:
string name;
double value;
Variable( string n, double v):name( n),value( v){}
};

const string prompt = ">";
const string result = "=";
const string declkey= "let";
const char let = 'L';
const char name = 'a';
const char print = ';';
const char quit = 'q';
const char number = '8';
vector <Variable> var_table;
Token_stream ts;



double expression();

double term();

double primary();

double get_value( string s);

void set_value( string s, double d);

double statement( );

bool is_declared( string var);

double define_name( string var, double val);

double declaration( );

void Token_stream::putback( Token t)
{
OutputDebugString("putback..");
if( full) error( "putback( ) into a full buffer");
buffer=t;
full=true;
}

Token Token_stream::get( )
{
if( full)
{
full=false;
return buffer;
}
char ch;
cin >> ch;
switch( ch){
case quit:
case print:
case '(':
case ')':
case '+':
case '-':
case '*':
case '/':
case '%':
OutputDebugString("get operator");
return Token( ch);
case '.':
case '0': case '1': case '2': case '3': case '4':
case '5': case '6': case '7': case '8': case '9':
{
cin.putback( ch);
double val;
cin >> val;
OutputDebugString("get number");
return Token( number,val);
}
case '=':
OutputDebugString("get '='");
return Token(ch);
default :
if( isalpha( ch))
{
string s;
s+=ch;
while( cin.get( ch) &&( isalpha( ch)||isdigit( ch)))
s+=ch;
if( s == declkey)
{
OutputDebugString("get 'let'");
return Token( let);
}
cin.putback(ch);
return Token( name, s);
}
error( "Bad token");
}
}

void Token_stream::ignore( char c)
{
if( full && c==buffer.kind)
{
full = false;
return;
}
full=false;

char ch=0;
while( cin>>ch)
if( ch==c) return;
}





double get_value( string s)
{
for( int i=0; i<var_table.size( ); ++i)
if( var_table[i].name==s)
{
return var_table[i].value;
}
error( "get : undefined Variable",s);
}

void set_value( string s, double d)
{
for( int i=0; i<var_table.size( );++i)
if( var_table[i].name==s)
{
var_table[i].value=d;
return;
}
error( "set : undefined Variable",s);

}

bool is_declared( string var)
{
for( int i=0; i<var_table.size( ); ++i)
if( var_table[i].name==var) return true;
return false;
}

double define_name( string var, double val)
{
if( is_declared( var)) error( var," declared twice");
var_table.push_back( Variable( var,val));
return val;
}

double declaration( )
{
Token t=ts.get( );
if( t.kind!=name) error( "name expected in declaration");
string var_name=t.name;

Token t2=ts.get( );

if( t2.kind!='=') error( "missing '=' in declaration of ",var_name);

double d = expression( );
define_name( var_name, d);
return d;
}

double statement( )
{
Token t=ts.get( );
switch( t.kind)
{
case let:
return declaration( );
default:
ts.putback( t);
return expression( );
}
}





double expression()
{
double left=term();
Token t=ts.get();
while( true)
{
switch(t.kind)
{
case '+' :
left+=term( );
t=ts.get( );
break;
case '-' :
left-=term( );
t=ts.get( );
break;
default:
ts.putback( t);
char buffer[30];
sprintf(buffer,"expression:return value:%f",left);
OutputDebugString(buffer);
return left;
}
}
}
double term()
{
double left=primary();
Token t=ts.get();
while( true)
{
switch(t.kind)
{
case '*' :
left*=primary( );
t=ts.get( );
break;
case '/' :
{

double d=primary( );
if( d == 0) error( "divide by zero");
left/=d;
t=ts.get( );
break;
}
case '%' :
{
int i1=narrow_cast<int>( left);
int i2=narrow_cast<int>(term());
if( i2==0) error( "%:divide by zero");
left=i1%i2;
t= ts.get( );
break;
}
default:
ts.putback( t);
char buffer[30];
sprintf(buffer,"term:return value:%f",left);
OutputDebugString(buffer);
return left;
}
}
}

double primary()
{
Token t=ts.get();
switch(t.kind)
{
case '(' :
{
double d=expression( );
t=ts.get( );
if( t.kind!=')') error("')' expected") ;
return d;
}
case number :
char buffer[30];
sprintf(buffer,"primary:return value:%f",t.value);
OutputDebugString(buffer);
return t.value;
case '-':
return - primary( );
case '+':
return primary( );
default:
error( "primary expected");
}

}



void clean_up_mess( )
{
ts.ignore( print);
//system("cls");
}

void calculate( )
{
while(cin)
{
try{
cout << prompt;
Token t = ts.get();
while( t.kind == print) t=ts.get( );
if( t.kind == quit)
{
return ;
}
ts.putback( t);
cout << result << statement( ) << endl;
}
catch( exception & e)
{
cerr << e.what( ) << endl;
clean_up_mess( );
}
}
}



int main()
{
try
{
calculate( );
keep_window_open( );
return 0;
}
catch( exception & e)
{
cerr << e.what( ) << endl;
keep_window_open( "~~");
return 1;
}
catch(...)
{
cerr << "exception \n" ;
keep_window_open( "~~");
return 2;
}
}



...全文
422 6 打赏 收藏 转发到动态 举报
写回复
用AI写文章
6 条回复
切换为时间正序
请发表友善的回复…
发表回复
cocat 2010-01-21
  • 打赏
  • 举报
回复
好长噢~帮顶
老吴笔记 2010-01-21
  • 打赏
  • 举报
回复
哪里越界嘞
不晓得,你试了下喽
hlyces 2010-01-21
  • 打赏
  • 举报
回复
...

我在VC下能得到正解


gcc下却说越界??


是不是double 要写成long double?
老吴笔记 2010-01-21
  • 打赏
  • 举报
回复
百度一下,你就知道
Google一下,你就知道的太多了
hlyces 2010-01-21
  • 打赏
  • 举报
回复

/*
simple "Programming: Principles and Practice using C++" course header to
be used for the first few weeks.
It provides the most common standard headers (in the global namespace)
and minimal exception/error support.

Students: please don't try to understand the details of headers just yet:
All will be explained. This header is primarily used so that you don't have
to understand every concept all at once.
*/

#ifndef H112
#define H112 201001L

#include<iostream>
#include<fstream>
#include<sstream>
#include<cmath>
#include<cstdlib>
#include<string>
#include<vector>
#include<algorithm>
#include<stdexcept>
using namespace std;

const string s="";
template<class T> string to_string(const T& t)
{
ostringstream os;
os << t;
return os.str();
}

struct Range_error : out_of_range { // enhanced vector range error reporting
int index;
Range_error(int i) :out_of_range("Range error: "+to_string(i)), index(i) { }
};




struct Exit : runtime_error {
Exit(): runtime_error("Exit") {}
};

// error() simply disguises throws:
inline void error(const string& s)
{
throw runtime_error(s);
}

inline void error(const string& s, const string& s2)
{
error(s+s2);
}

inline void error(const string& s, int i)
{
ostringstream os;
os << s <<": " << i;
error(os.str());
}



template<class T> char* as_bytes(T& i) // needed for binary I/O
{
void* addr = &i; // get the address of the first byte
// of memory used to store the object
return static_cast<char*>(addr); // treat that memory as bytes
}


inline void keep_window_open()
{
// cin.clear();
// cout << "Please enter a charcter to exit\n";
system("pause");
// char ch;
// cin >> ch;
return;
}

inline void keep_window_open(string s)
{
if (s=="") return;
cin.clear();
cin.ignore(120,'\n');
for (;;) {
cout << "Please enter " << s << " to exit\n";
string ss;
while (cin >> ss && ss!=s)
cout << "Please enter " << s << " to exit\n";
return;
}
}

// make std::min() and std::max() accessible:
#undef min
#undef max

#include<iomanip>
inline ios_base& general(ios_base& b) // to augment fixed and scientific
{
b.setf(ios_base::fmtflags(0),ios_base::floatfield);
return b;
}

// run-time checked narrowing cast (type conversion):
template<class R, class A> R narrow_cast(const A& a)
{
R r = a;
if (A(r)!=a) error(string("info loss"));
return r;
}


inline int randint(int max) { return rand()%max; }

inline int randint(int min, int max) { return randint(max-min)+min; }

#endif
hlyces 2010-01-21
  • 打赏
  • 举报
回复
[Quote=引用 4 楼 laowu_ 的回复:]
哪里越界嘞
不晓得,你试了下喽
[/Quote]
在gcc下
你有没试??
输入1+2;
会得到=nan
而在VC下是
=3
内容概要:本文围绕“非线性流量的数据驱动Koopman模型预测控制研究”展开,提出一种基于数据驱动的Koopman算子理论方法,用于构建非线性系统的线性化状态空间模型,并结合模型预测控制(MPC)实现对复杂非线性系统的高效控制。研究通过引入扩展动态模态分解(EDMD)等观测函数,将非线性动力学映射至高维特征空间,在该空间中实现近似线性化表征,进而融合线性MPC框架进行优化求解。全文系统阐述了Koopman算子的数学基础、隐式线性化机制及在非线性流量控制中的建模流程,并通过Matlab代码完成了算法实现与仿真实验,验证了该方法在处理无精确物理模型、强非线性、时变动态系统中的有效性与鲁棒性,尤其适用于工业流程控制、能源系统调度等实际工程场景。; 适合人群:具备自动控制理论、非线性系统分析基础,熟悉Matlab编程,从事控制工程、系统辨识、智能优化、能源系统建模等方向的研究生、科研人员及工程技术人员。; 使用场景及目标:①应用于难以建立精确数学模型的复杂非线性系统(如流体动力系统、电力电子系统、机器人动力学等)的建模与实时控制;②实现数据驱动下的模型预测控制,提升系统响应速度与控制精度;③为先进控制策略(如MPC)提供一种可行的线性化建模范式,推动现代控制理论与数据科学、机器学习的深度融合。; 阅读建议:建议读者结合提供的Matlab代码深入理解Koopman方法的具体实现过程,重点关注观测函数构造、核函数选择、矩阵逼近、降维处理及MPC控制器设计等关键技术环节,并尝试将其迁移至其他非线性系统中进行复现实验与性能对比,以全面掌握其适用范围与局限性。
内容概要:本文详细介绍了一种基于Simulink的光伏储能单相逆变器并网仿真模型,系统涵盖了光伏阵列、储能单元、DC-AC单相逆变器及并网接口的完整结构,重点实现了储能环节的能量管理与逆变器并网控制策略的建模仿真。通过Simulink平台构建系统模型,验证了逆变器输出电能质量、并网稳定性以及控制系统的动态响应性能,采用SPWM调制、PI闭环控制等关键技术,确保并网电流与电网电压同频同相,满足并网电能质量要求。该模型不仅可用于分布式能源系统的仿真研究,还可作为新能源并网技术的教学与工程实践工具。; 适合人群:电气工程、自动化、新能源科学与工程等相关专业的高校本科生、研究生、科研人员,以及从事光伏发电系统设计、储能控制与并网技术研发的工程技术人员。; 使用场景及目标:①深入理解光伏储能系统中能量转换、存储与并网控制的整体工作原理;②支持课程设计、毕业设计或科研项目中对单相逆变器控制策略(如SPWM、PI调节、锁相技术等)的仿真验证与参数优化;③为后续研究更复杂的控制算法(如MPPT、低电压穿越、谐波抑制等)提供可扩展的仿真基础平台。; 阅读建议:建议结合MATLAB/Simulink环境动手搭建与调试模型,逐步理解各模块(如光伏建模、储能充放电控制、逆变器驱动、锁相环、PI调节器等)的功能与交互关系,重点关注控制系统的设计逻辑与参数整定过程,并可通过修改负载条件或电网参数测试系统鲁棒性,为进一步拓展至三相系统或多机并网场景奠定基础。

65,212

社区成员

发帖
与我相关
我的任务
社区描述
C++ 语言相关问题讨论,技术干货分享,前沿动态等
c++ 技术论坛(原bbs)
社区管理员
  • C++ 语言社区
  • encoderlee
  • paschen
加入社区
  • 近7日
  • 近30日
  • 至今
社区公告
  1. 请不要发布与C++技术无关的贴子
  2. 请不要发布与技术无关的招聘、广告的帖子
  3. 请尽可能的描述清楚你的问题,如果涉及到代码请尽可能的格式化一下

试试用AI创作助手写篇文章吧