资源简介
DANDAOFANGZHEN.m
代码片段和文件信息
syms y0 y1 y2 y3 y4 y5 x x0 x1 x2 x3 x4 x5 a Ma
t0=(x-x1)*(x-x2)*(x-x3)*(x-x4)*y0/((x0-x1)*(x0-x2)*(x0-x3)*(x0-x4)) %对于阻力因数的线性差值,采用4阶拉格朗日插值法
t1=(x-x0)*(x-x2)*(x-x3)*(x-x4)*y1/((x1-x0)*(x1-x2)*(x1-x3)*(x1-x4))
t2=(x-x1)*(x-x0)*(x-x3)*(x-x4)*y2/((x2-x1)*(x2-x0)*(x2-x3)*(x2-x4))
t3=(x-x1)*(x-x2)*(x-x0)*(x-x4)*y3/((x3-x1)*(x3-x2)*(x3-x0)*(x3-x4))
t4=(x-x1)*(x-x2)*(x-x3)*(x-x0)*y4/((x4-x1)*(x4-x2)*(x4-x3)*(x4-x0))
Cx0=t0+t1+t2+t3+t4
Cx15=subs(Cx0{xx0x1x2x3x4y0y1y2y3y4}{a2468100.0430.05110.06510.08470.112})
Cx21=subs(Cx0{xx0x1x2x3x4y0y1y2y3y4}{a2468100.0360.04360.05580.07360.0973})
Cx27=subs(Cx0{xx0x1x2x3x4y0y1y2y3y4}{a2468100.03080.03720.04810.06410.0849})
Cx33=subs(Cx0{xx0x1x2x3x4y0y1y2y3y4}{a2468100.02650.03230.04190.0560.0746})
Cx40=subs(Cx0{xx0x1x2x3x4y0y1y2y3y4}{a2468100.02220.02720.03560.04780.0644})
Cx=subs(Cx0{xx0x1x2x3x4y0y1y2y3y4}{Ma1.52.12.73.34.0Cx15Cx21Cx27Cx33Cx40})
f0=(x-x1)*(x-x2)*(x-x3)*(x-x4)*(x-x5)*y0/((x0-x1)*(x0-x2)*(x0-x3)*(x0-x4)*(x0-x5)) %对于升力因数的线性差值,采用5阶拉格朗日插值法
f1=(x-x0)*(x-x2)*(x-x3)*(x-x4)*(x-x5)*y1/((x1-x0)*(x1-x2)*(x1-x3)*(x1-x4)*(x1-x5))
f2=(x-x1)*(x-x0)*(x-x3)*(x-x4)*(x-x5)*y2/((x2-x1)*(x2-x0)*(x2-x3)*(x2-x4)*(x2-x5))
f3=(x-x1)*(x-x2)*(x-x0)*(x-x4)*(x-x5)*y3/((x3-x1)*(x3-x2)*(x3-x0)*(x3-x4)*(x3-x5))
f4=(x-x1)*(x-x2)*(x-x3)*(x-x0)*(x-x5)*y4/((x4-x1)*(x4-x2)*(x4-x3)*(x4-x0)*(x4-x5))
f5=(x-x1)*(x-x2)*(x-x3)*(x-x4)*(x-x0)*y5/((x5-x1)*(x5-x2)*(x5-x3)*(x5-x4)*(x5-x0))
Cya0=f0+f1+f2+f3+f4+f5
Cya15=subs(Cya0{xx0x1x2x3x4x5y0y1y2y3y4y5}{a12468100.03020.03040.03060.03090.03110.0313})
Cya20=subs(Cya0{xx0x1x2x3x4x5y0y1y2y3y4y5}{a12468100.02790.02800.02840.02860.02880.029})
Cya25=subs(Cya0{xx0x1x2x3x4x5y0y1y2y3y4y5}{a12468100.02610.02640.02670.02690.02720.0274})
Cya30=subs(Cya0{xx0x1x2x3x4x5y0y1y2y3y4y5}{a12468100.02470.02480.02510.02540.02570.0259})
Cya35=subs(Cya0{xx0x1x2x3x4x5y0y1y2y3y4y5}{a12468100.02260.02270.02310.02330.02360.0238})
Cya40=subs(Cya0{xx0x1x2x3x4x5y0y1y2y3y4y5}{a12468100.02090.0210.02130.02160.02190.0221})
Cya=subs(Cya0{xx0x1x2x3x4x5y0y1y2y3y4y5}{Ma1.522.533.54Cya15Cya20Cya25Cya30Cya35Cya40})
mu=zeros(1100); %为迭代量预置空间,
Y=zeros(1100);
theta
- 上一篇:图像处理识别图中棒材个数
- 下一篇:基于matlab,读取图像文件并计算其图像对比度
相关资源
- 《MATLAB仿真在现代通信中的应用》徐
- KECA1.m
- zw_支持向量机与朴素贝叶斯算法matl
- test1_mat.m
- Finnal6.m
- odefun1.m
-
飞控simuli
nk仿真.zip -
光伏发电最大功率点跟踪simuli
nk仿真 - MATLAB全球磁场绘图程序.zip
- MATLAB人脸识别考勤系统摄像头,记录
- MATLAB多方法车牌识别系统bp+模板+GUI
- 四旋翼MATLAB仿真文件实验报告.zip
-
4.1基于Simuli
nk的图像、视频处理.z -
Simuli
nk的低通数字滤波器的仿真分析 - matlabR2017b迅雷种子文件及安装教程.
- 最优化方法及其matlab程序实现.zip
- 5d5680a96365f0e9a0c7114890f97460.mobi
- 《全面详解LTE:MATLAB建模、仿真与实
- XXXXXXRobustControlDesignwithMATLAB.pdf
- matlab图像处理包括数据反演以及其他
- Matlab常用32种算法.zip
- 传感器网络节点定位算法matlab代码.
- 基于DCT变换和DFT变换的数字图像压缩
- 10900257雷达系统设计MATLAB仿真.pdf
- 神经网络、遗传算法、支持向量机、
- 《先进PID控制MATLAB仿真第3版》仿真程
- 雷达系统设计MATLAB仿真.pdf196612
- 滑模变结构控制MATLAB仿真(第1版)(
- [Matlab作品]matlab实现视频中动态目标跟
- zw_qq_15023225-7501331-S函数的BP神经网络
评论
共有 条评论