#include<stdlib.h>
#define NUM 20
void main()
{
double sum //总和
double fez //分子
double fem //分母
int i//循环计算到NUM
int n//单项算分子分母
for( sum=0,i=1i<=NUMi++ )
{
for( fez=1,fem=1,n=1n<in++)
{
fez*=n
fem*=(2*(n+1)-1)
}
printf("%-4d %-2.6f %-2.6f\n",n,fez,fem)
sum+=fez/fem
}
printf("总和:%f\n",sum)
}
MATLAB求悬臂梁自由端挠度主要用数值分析的原理梁模型有Euler-Bernoulli, Timoshiko等
数值方法有很多,FDM,FEM,打靶等等,不同模型不同方法都有不同的程序,建议先搞清楚需要用什么模型和方法
当然也可以用带符号变量的办法来用MATLAB直接解微分方程,求的相对的一个解析解,不过这样手算部分会比较多
欢迎分享,转载请注明来源:内存溢出
评论列表(0条)