交通规划福尼斯法计算出行分布交通量

最新推荐文章于 2024-06-24 17:58:39 发布

原创最新推荐文章于 2024-06-24 17:58:39 发布 · 1.2k 阅读

8 ·

CC 4.0 BY-SA版权

文章标签：

#c语言 #矩阵

本文探讨了如何运用C语言编程技术和矩阵运算来实施福尼斯法，以计算城市交通规划中的出行分布流量。通过理解交通网络的结构，建立矩阵模型，实现了对交通流量的精确计算和分配。

交通规划福尼斯法计算出行分布交通量

做交通规划作业福尼斯法时实在不想手算了，写了个c语言代码来完成。
代码不难，主要是为了记录自己的学习过程。

#include <stdio.h>
#include <math.h>
double OD[3][3]={
    {17.0, 7.0, 4.0},
    {7.0, 38.0, 10.0},
    {5.0, 5.0, 17.0}
    };//初始值
double e = 0.03;//误差常数
double u[3]={41.7, 88.8, 38.0};
double v[3]={41.3, 89.4, 37.8};//
double T;//总量
double O[3];
double D[3];
double fo[3];//发生增长系数
double fd[3];//吸引增长系数
void GetODT(){
    for (int i = 0; i < 3; i++){
        for (int j = 0; j < 3; j++)
        {
            O[i]=0;
            D[j]=0;
        }
        T=0;
    }
    for (int i = 0; i < 3; i++){
        for (int j = 0; j < 3; j++)
        {
            O[i]+=OD[i][j];
            D[j]+=OD[i][j];
            T+=OD[i][j];
        }
    }
}//求和计算O、D、T的值
void PrintODT(){
    for (int i = 0; i < 3; i++){
        for (int j = 0; j < 3; j++){
            printf("OD[%d][%d] = %.4f\t", i,j, OD[i][j] );
            }
        printf("O[%d] = %.4f\t", i, O[i]);
        printf("\n");
    }
    for (int j = 0; j < 3; j++){
    printf("  D[%d] = %.4f\t", j, D[j]);
    }
    printf("  T = %.4f\t", T);
    printf("\n");
    printf("\n");
}//打印OD表
int main(){
int flag=1;
int m=1;
int f=1;
while(1){
    m=0;
    GetODT();
    PrintODT();
    for (int i=0; i<3; i++){
        fo[i] = u[i]/O[i];
        printf("fo[%d]=%.4f\t", i, fo[i]);
        for (int j=0; j < 3; j++){
            OD[i][j]*=fo[i];
            if(fabs(fo[i]-1)>=e){
                m=1;
            }
        }
    }
    printf("\n");
    if(m==0&&f==0){
        flag=0;
        break;
    }
    GetODT();
    PrintODT();
    f=0;
    for (int j=0; j<3; j++){
        fd[j] = v[j]/D[j];
        printf("fd[%d]=%.4f\t", j, fd[j]);
        for (int i=0; i < 3; i++){
            OD[i][j]*=fd[j];
            if(fabs(fd[j]-1)>=e){
                f=1;
            }
        }
    }
    printf("\n");
    if(m==0&&f==0){
        flag=0;
        break;
    }
}
return 0;
}