贷款内部收益率IRR计算表格、程序和代码

// 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说明找到符合要求的IRR if (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; }