Calibracion de sensor mq131

Hola a todos!

Estoy intentando hacer funcionar el sensor mq131 y he leído en muchos foros me sigue quedando un vacio muy grande, ya que de este sensor no hay nada o casi nada. No tengo claro como hay que calibrarlo, además yo estoy probando con resistencia de 100kilo hasta 200kilo como todo el mundo pone en los foros. Si pongo la de 100k me da unos valores analógicos comprendidos entre 500 hasta 600 mas o menos y si pongo la de 200k me de unos valores de 700 hasta 800. Y si expongo mi sensor a gas ozono esos valores decrecen siempre. Luego yo estoy metiendo este programa para calibrarlo pero siempre me da 0ppb y eso es porque no tengo ni idea de como coger la tabla de calibración y me gustaría haber si alguien puede ayudarme y decirme que mal tengo en mi programa que pongo a continuación:

************************Hardware Related Macros************************************/
#define         MQ_PIN                       (0)     //define which analog input channel you are going to use
#define         RL_VALUE                     (100)    //define the load resistance on the board, in kilo ohms
#define         RO_CLEAN_AIR_FACTOR          (9.83)  //RO_CLEAR_AIR_FACTOR=(Sensor resistance in clean air)/RO,
                                                    //which is derived from the chart in datasheet

/***********************Software Related Macros************************************/
#define         CALIBARAION_SAMPLE_TIMES     (50)    //define how many samples you are going to take in the calibration phase
#define         CALIBRATION_SAMPLE_INTERVAL  (500)   //define the time interal(in milisecond) between each samples in the
                                                    //cablibration phase
#define         READ_SAMPLE_INTERVAL         (50)    //define how many samples you are going to take in normal operation
#define         READ_SAMPLE_TIMES            (5)     //define the time interal(in milisecond) between each samples in 
                                                    //normal operation

/**********************Application Related Macros**********************************/
#define         GAS_O3                      (0)

/*****************************Globals***********************************************/
float           O3Curve[3]  =  {0.69, 0.69,-0.76};    //two points are taken from the curve in datasheet. 
                                                    //with these two points, a line is formed which is "approximately equivalent" 
                                                    //to the original curve. 
                                                    //data format:{ x, y, slope}; point1: (lg200, lg8.5), point2: (lg10000, lg0.03) 

float           Ro           =  10;                  //Ro is initialized to 10 kilo ohms

void setup()
{
 Serial.begin(9600);                                //UART setup, baudrate = 9600bps
 Serial.print("Calibrating...\n");                
 Ro = MQCalibration(MQ_PIN);                        //Calibrating the sensor. Please make sure the sensor is in clean air 
                                                    //when you perform the calibration                    
 Serial.print("Calibration is done...\n"); 
 Serial.print("Ro=");
 Serial.print(Ro);
 Serial.print("kohm");
 Serial.print("\n");
}

void loop()
{
  Serial.print("O3:"); 
  Serial.print(MQGetGasPercentage(MQRead(MQ_PIN)/Ro,GAS_O3) );
  Serial.print( "ppb" );
  Serial.print("\n");
  delay(200);
}

/****************** MQResistanceCalculation ****************************************
Input:   raw_adc - raw value read from adc, which represents the voltage
Output:  the calculated sensor resistance
Remarks: The sensor and the load resistor forms a voltage divider. Given the voltage
        across the load resistor and its resistance, the resistance of the sensor
        could be derived.
************************************************************************************/ 
float MQResistanceCalculation(int raw_adc)
{
 return ( ((float)RL_VALUE*(1023-raw_adc)/raw_adc));
}

/***************************** MQCalibration ****************************************
Input:   mq_pin - analog channel
Output:  Ro of the sensor
Remarks: This function assumes that the sensor is in clean air. It use  
        MQResistanceCalculation to calculates the sensor resistance in clean air 
        and then divides it with RO_CLEAN_AIR_FACTOR. RO_CLEAN_AIR_FACTOR is about 
        10, which differs slightly between different sensors.
************************************************************************************/ 
float MQCalibration(int mq_pin)
{
 int i;
 float val=0;

 for (i=0;i<CALIBARAION_SAMPLE_TIMES;i++) {            //take multiple samples
   val += MQResistanceCalculation(analogRead(mq_pin));
   delay(CALIBRATION_SAMPLE_INTERVAL);
 }
 val = val/CALIBARAION_SAMPLE_TIMES;                   //calculate the average value

 val = val/RO_CLEAN_AIR_FACTOR;                        //divided by RO_CLEAN_AIR_FACTOR yields the Ro 
                                                       //according to the chart in the datasheet 

 return val; 
}
/*****************************  MQRead *********************************************
Input:   mq_pin - analog channel
Output:  Rs of the sensor
Remarks: This function use MQResistanceCalculation to caculate the sensor resistenc (Rs).
        The Rs changes as the sensor is in the different consentration of the target
        gas. The sample times and the time interval between samples could be configured
        by changing the definition of the macros.
************************************************************************************/ 
float MQRead(int mq_pin)
{
 int i;
 float rs=0;

 for (i=0;i<READ_SAMPLE_TIMES;i++) {
   rs += MQResistanceCalculation(analogRead(mq_pin));
   delay(READ_SAMPLE_INTERVAL);
 }

 rs = rs/READ_SAMPLE_TIMES;

 return rs;  
}

/*****************************  MQGetGasPercentage **********************************
Input:   rs_ro_ratio - Rs divided by Ro
        gas_id      - target gas type
Output:  ppm of the target gas
Remarks: This function passes different curves to the MQGetPercentage function which 
        calculates the ppm (parts per million) of the target gas.
************************************************************************************/ 
int MQGetGasPercentage(float rs_ro_ratio, int gas_id)
{
 if ( gas_id == GAS_O3) {
    return MQGetPercentage(rs_ro_ratio,O3Curve);
 }  
 return 0;
}

/*****************************  MQGetPercentage **********************************
Input:   rs_ro_ratio - Rs divided by Ro
        pcurve      - pointer to the curve of the target gas
Output:  ppm of the target gas
Remarks: By using the slope and a point of the line. The x(logarithmic value of ppm) 
        of the line could be derived if y(rs_ro_ratio) is provided. As it is a 
        logarithmic coordinate, power of 10 is used to convert the result to non-logarithmic 
        value.
************************************************************************************/ 
int  MQGetPercentage(float rs_ro_ratio, float *pcurve)
{
 return (pow(10,( ((log(rs_ro_ratio)-pcurve[1])/pcurve[2]) + pcurve[0])));
}

Espero que alguien pueda ayudarme.
Un saludo