请大家给我看下这个程序

xindao

某线性系统（一种速度反馈液压放大器的简化），微分方程为：

1. dy=5*((u[t]+0.5)*sqrt(8-y[t])-sqrt(y[t]+2))/(sqrt(8-y[t])-2);
   
2. for(i=0;i<100;i++)u[i];
   
3. for(j=0;j<100;j++)
   
4. {
   
5.     for(i=0;i<j+1;i++)
   
6.         u[i]=0.1;
   
7.     ApartTwoVolterra(u,y2);//各阶输出隔离
   
8.     for(i=0;i<100-j;i++)
   
9.       {
   
10.         b=y2[i+1+j]/0.01/pow(0.01,2);
    
11.         for(k=I;k<i+j+1;k++)
    
12.         for(h=k;h<i+j+1;h++)
    
13.             b=b-jie2[k*100+h];
    
14.           jie2[i*100+i+j]=jie[(i+j)*100+i]=b/2;
    
15.        }
    
16.   }
    
17.   Void ApartTwoVolterra(double*u,double*y)
    
18.    {//省略变量定义
    
19.          for(j=0;j<10;j++)
    
20.           {
    
21.             for(i=0;i<Flong;i++)u[i]=u[i]*1.2
    
22.             for(i=0.y[0]=0;i<Flong;i++)
    
23.             {  
    
24.               dy=5*((u[t]+0.5)*sqrt(8-y[t])-sqrt(y[t]+2))/(sqrt(8-y[t])-2);
    
25.               y[i+1]=y[i]+dy*0.01;//简化的仿真过程
    
26.               y1[j*Flong+i]=y[i];
    
27.              }
    
28.          }
    
29.        for(i=o;i<10;i++)
    
30.        for(j=0;j<6;j++)
    
31.           a[i*6+j]=pow(pow(1.2,i),j+1);
    
32.        for(i=0;i<Flong;++)
    
33.        {
    
34.           for(j=0;j<10;j++)
    
35.               b[j]=y1[j*Flong+i];
    
36.           equation(a,b,x,10,6);//调用10X6的线性方程组的最小二乘解函数
    
37.           y[i]=x[1];
    
38.         }
    
39.    }

复制代码

这是Volterra函数二阶核的简化求解程序。我看不是很懂。希望大侠们帮忙给我做个完整一点的。感激不尽！

[ 本帖最后由 风花雪月 于 2008-8-21 10:26 编辑 ]

xindao

这里是一些算法思想

express

按LSD意见办理

风花雪月

这是Volterra函数二阶核的简化求解程序。我看不是很懂。希望大侠们帮忙给我做个完整一点的。感激不尽！

没明白你的意思