⛄一、获取代码方式
获取代码方式1:
完整代码已上传我的资源:【气动学】基于matlab涡格法气动力计算【含Matlab源码 3038期】
点击上面蓝色字体,直接付费下载,即可。
获取代码方式2:
付费专栏Matlab物理应用(初级版)
备注:
点击上面蓝色字体付费专栏Matlab物理应用(初级版),扫描上面二维码,付费29.9元订阅海神之光博客付费专栏Matlab物理应用(初级版),凭支付凭证,私信博主,可免费获得1份本博客上传CSDN资源代码(有效期为订阅日起,三天内有效);
点击CSDN资源下载链接:1份本博客上传CSDN资源代码
⛄二、部分源代码
%wogefa zhaoxian%
clear
clc
alpha=2;
l0=input(‘请输入翼根弦长’)
ll=input(‘请输入翼梢弦长’)
L=input(‘请输入展长’)
ka=input(‘请输入前缘后掠角(以弧度计量)’)
ka=ka100000;
ka=round(ka);
ka=ka/100000;
S=(ll+l0)0.5L;
bA=2/S(L/6*(l0-ll)2+L/2*ll*(l0-ll)+L/2*ll2);
u=50;
tic
[Q,XZ,X1,X2,Z1,Z2,W]=fenkuai(l0,ll,L,ka,u);
%列方程组%
[T,qi]=qiujie(XZ,X1,X2,Z1,Z2,u);
%以下部分求解方程组%
v=100sin(alpha/180pi);
V(1:2uu)=v4pi;
H=inv(T)V’;
P(u,2u)=0;
for i=1:1:u
for j=1:1:2u
P(i,j)=H(2u*(i-1)+j);
end
end
%以下部分计算答案%
zi(1:u,1:2u)=L/2/u;
F=1.112100sum(sum(P.zi));
Cy=F/(0.51.112S100^2);
M=1.112100sum(sum(P.zi.W));
Mz=1.112100sum(sum(P.zi.W))/(0.51.112SbA100^2);
fprintf(‘升力系数为’)
Cy
fprintf(‘力矩系数为’)
mz=M/(0.51.112SbA*100^2)
fprintf(‘焦点位置为’)
XXX=Mz/Cy
fprintf(‘奇点数量为’)
qi
toc
⛄三、运行结果
⛄四、matlab版本及参考文献
1 matlab版本
2014a
2 参考文献
[1] 门云阁.MATLAB物理计算与可视化[M].清华大学出版社,2013.
3 备注
简介此部分摘自互联网,仅供参考,若侵权,联系删除