mtufilterread.cpp

Go to the documentation of this file.
00001 #include <fstream>
00002 #include <iostream>
00003 #include <vector>
00004 #include <string>
00005 #include <complex>
00006 
00007 using namespace std;
00008 
00009 int main()
00010 {
00011         const complex<double> I(0.0,1.0);
00012         ifstream filterfile;
00013         ofstream phaseresponse;
00014         string filename;
00015         double currreal, currimag, currfreq;
00016         
00017         vector<double> freqs;
00018         vector< complex<double> > hxresponse;
00019         double dummy;
00020         char line[255];
00021         
00022         cout << "Filename: ";
00023         cin >> filename;
00024         filterfile.open(filename.c_str());
00025         filterfile.getline((char *)line,255);
00026         while (filterfile.good())
00027         {
00028                 filterfile >> currfreq;         
00029                 filterfile.ignore(256,',');
00030                 filterfile.ignore(256,',');
00031                 filterfile.ignore(256,',');
00032                 filterfile.ignore(256,',');
00033                 filterfile.ignore(256,',');
00034                 filterfile.ignore(256,',');
00035                 filterfile >> currreal;
00036                 filterfile.ignore(256,',');
00037                 filterfile >> currimag;
00038                 cout << currfreq << " " << currreal << " " << currimag << endl;
00039                 filterfile.getline((char *)line,255);
00040                 if (filterfile.good())
00041                 {
00042                         freqs.push_back(currfreq);
00043                         hxresponse.push_back(currreal + I * currimag);
00044                         
00045                 }               
00046         }
00047         filterfile.close();
00048         phaseresponse.open((filename+".resp").c_str());
00049         for (int i = 0; i < freqs.size() ; ++i)
00050                 phaseresponse << freqs.at(i) << " " << arg(hxresponse.at(i)) * 180./3.1415 << endl;
00051         phaseresponse.close();
00052 }