// testirr.cpp : 定义控制台应用程序的入口点。
//

#include "stdafx.h"
#include <math.h>//const double zero = 1e-5;//int n;/*
//代码是在csdn上找到的 我拿下来改了一下
//https://blog.csdn.net/dinghaoseu/article/details/50322117
//这是原来的代码
double quickpow(double a, int b)
{double ans = 1;while(b){if(b & 1) ans = ans * a;a = a * a;b >>= 1;}return ans;
}double makeans(double irr)
{double ans = 0;for(int i = 1; i <= n; i++)ans += ( a[i] / (quickpow(1 + irr, i)) );return ans;
}
*/
//计算流出Npv
double getNpvOut(double irr, const double arrOutMoney[], int n)
{double npv = 0;double denominator = 1 + irr; //分母double multiplier = denominator; //乘数for (int i = 1; i <= n; ++i){npv += arrOutMoney[i] / denominator; //即：arrOutMoney[i] / (quickpow(1 + irr, i)denominator *= multiplier;}return npv;
}//最小最大IRR（内部收益率） 用于折算年化率计算
//最小值必须大于-1 最大值1其实就足够了，折算成年化收益率是120%（国家规定利率不能超过30%）
//数字越小计算速度越快，为保险计这里填10
#define IRR_MIN -1.0f
#define IRR_MAX 10.0f//二分迭代寻找合适的irr值（如果存在多个irr，取第一次找到的值，不保证大小顺序）
//arrInOutMoney是分期现金流 nArrLen是arrInOutMoney元素个数
//arrInOutMoney[0]必须是负数，代表总分期金额（负数），后面是每期还款金额（正数）
//返回[IRR_MIN, IRR_MAX]之间的数字代表符合要求的IRR值，<IRR_MIN代表找不到合适的IRR
double binarySearchGetIrr(const double arrInOutMoney[], int nArrLen)
{double l = IRR_MIN, r = IRR_MAX; //irr取值-1~10之间 int n = nArrLen - 1;//int nCnt = 0;while (l < r){//每次从最大和最小期望IRR中间取一个数值进行npv测算double mid = (l + r) / 2;//现在是用除法求Npv，其实可以改成用乘法求，效率会高一点//因为计算机处理乘法速度比较快double npvOut = getNpvOut(mid, arrInOutMoney, n);//++nCnt;//如果结果等于0说明找到符合要求的IRRif (fabs(npvOut + arrInOutMoney[0]) <= 1e-5) //double类型不能直接与0比较判断是否相等
        {    //printf("nCnt = %d\n", nCnt); //经测试，一般分12期迭代次数在30~60之间return mid;}//irr越大，npvOut越小，故npvOut太大时irr就应该落在mid和r之间，反之则反之else if (npvOut > -arrInOutMoney[0]) {l = mid;}else{r = mid;}}//printf("nCnt = %d\n", nCnt);//找不到返回比IRR_MIN还小的值return IRR_MIN - 1;
}//输入按月分期现金流，输出对应irr和折算年化收益率
//arrInOutMoney是分期现金流 nArrLen是arrInOutMoney元素个数（一般是7或者13）
//arrInOutMoney[0]必须是负数，代表总分期金额，后面是每期还款金额（正数）
//nCheckFlag = 0，代表直接计算irr，否则会先对数据合法性做检查
//返回0代表成功 其它代表失败 失败原因存放在errBuf（调用者需要保证至少有256个字节空间）
int GetIrrAndAnnualizedRate(const double arrInOutMoney[], int nArrLen,OUT double *pIrr, OUT double *pAnnualizedRate, int nCheckFlag, OUT char *errBuf)
{double irr = 0;if (nCheckFlag != 0){double inMoney, outMoney;int i;if (arrInOutMoney == NULL || nArrLen < 2){strcpy(errBuf, "arrInOutMoney需要非空并且元素个数大于2个");return -10;}inMoney = arrInOutMoney[0];outMoney = 0;for (i = 1; i < nArrLen; ++i){if (arrInOutMoney[i] < 0){//不支持多次现金流入（因为没这个需求，不要浪费计算力）outMoney = -1;break;}outMoney += arrInOutMoney[i];}if (inMoney >= 0 || (-inMoney > outMoney)){strcpy(errBuf, "第一个元素必须是负现金流，之后每个元素均是正现金流，""并且正现金流之和要大于负现金流");return -20;}}irr = binarySearchGetIrr(arrInOutMoney, nArrLen);if (irr < IRR_MIN){sprintf(errBuf, "%.5f(%.2f%%)~%.5f(%.2f%%)之间无法找到合适的irr，请检查现金流是否输入异常",IRR_MIN, IRR_MIN * 12 * 100, IRR_MAX, IRR_MAX * 12 * 100);return -30;}*pIrr = irr;*pAnnualizedRate = irr * 12 * 100;return 0;
}int getIrrDemo()
{double irr, annualizedRate;double a[100 * 12 + 1];int n, nRet;char errBuf[256];printf("**************如果要退出，请在还款期数填0**************\n");while ((printf("input 还款期数 n（0代表退出）:")) && ~scanf("%d", &n) && n){printf("n = %d\n", n);if (n >= sizeof(a) / sizeof(a[0])){printf("n值太大，不支持\n");continue;}printf("输入分期金额（负数）:");if (scanf("%lf", &a[0]) != 1) {printf("输入的金额不能包含非数字和小数点\n");getchar();continue;}printf("输入%d期还款金额（正数），每输入一期按一次回车:", n);for (int i = 1; i <= n; i++){scanf("%lf", &a[i]);}nRet = GetIrrAndAnnualizedRate(a, n + 1, &irr, &annualizedRate, 1, errBuf);if (nRet != 0){printf("error:[%s]\n", errBuf);continue;}//计算irr常用的方法是迭代计算，即不断尝试可能值，根据尝试结果缩小范围，直到找到符合要求的值//网上能找到的迭代算法有二分迭代，牛顿迭代，加速迭代，其中二分迭代最好理解，最容易开发irr = binarySearchGetIrr(a, n + 1);if (irr < IRR_MIN){printf("找不到合适的irr\n");}else{printf("irr = %.6f  年化收益率(12 * irr) = %.4f%%\n", irr, annualizedRate);}}return 0;
}int main(int argc, char *argv[])
{getIrrDemo();//system("pause");return 0;
}