组员：xxxxxxx
matlab 程序设计题目。
20分）1） 以组为单位，自选一个或者自拟一个物理或者数学问题。
2） 利用matlab编写一个函数文件解决这个物理或者数学问题。
3） 要求：明确函数文件的输入变量与输出变量；利用简单算例对函数文件进行验证；必须用到for循环结构和if选择结构；程序命令不能少于30行。
4） 最后对该问题进行总结说明为什么要选这个题，并且谈谈通过做题的学习感受。
问题模型：已知一座桥，桥长l,桥宽l，桥厚d（厚度不变），判断一个重为g（g=f）的汽车在通过桥的整个过程中，桥的最大正应力是否满足桥的设计要求（桥的许用正应力为a）。并求出汽车在通过整个桥的过程中桥的最大正应力及此时车的位置以及最大正应力出现的位置。
桥的简化模型如下：
选题原因：切合实际，能够解决实际生活中的问题，对已知的桥进行安全检验。
做题感受：matlab在生活及工程应用中起很大作用，能够大大简化人得计算量，通过设计题目，我们发现自身的知识还需增强。
问题模型的验证（两个）及程序如下：
检验一：检验二：
程序设计：function y=fun(l,d1,d,l,f,a) %l 桥宽，d1 桥厚，d单位长度，l 桥长，f 车重 a 许用切应力。
桥厚等于桥长除以单位长度再加一
p=l./d+1;%3-5行表示桥厚向量（以d1为等厚）
q=ones(1,p);
d=q.*d1;
iz=l.*d.*d.*d./12; %iz 截面对中性轴的惯性矩。
e=d./2;
w=iz./e; %w截面抗弯矩系数。
x=[0:d:l];
for i= 1:1:(length(x)-1);%力作用在不同点。
for j= 1:1:length(x)-1;%不同点的正应力。
if j<=i %j点处的正应力。
a(j)=f.*(l-x(i)).x(j)./l.*x(j));
elsea(j)=(l-x(i)).f.*x(j)./l-f.*(x(j)-x(i)).w(j);
endend
[am1,m1]=max(a);%am1：f作用在i点处桥上的最大正应力 m1：该最大正应力的位置。
am(i)=am1;
m(i)=m1;
endamm,n]=max(am);%amm:f作用在不同点处的最大正应力，n：该最大正应力时车的位置。
mm=m(n);%mm:amm取最大值时最大正应力的位置。
y(1)=amm;%f作用在不同点处的最大正应力。
y(2)=n-1;%该最大正应力时车的位置。
y(3)=m(n)-1;%amm取最大值时最大正应力的位置。
if amm>a%判断车能否安全通过桥。
disp('the bridge is in danger')
elsedisp(' the bridge is in the security range')
endend
MATLAB程序设计基础作业
1.求下列表达式的值，观察matlab工作空间使用情况并保存全部变量。2 其中。提示 用冒号表达式生成向量，求各点函数值时用点乘运算 2.已知，求下列表达式的值 1 和 其中为单位阵 2 和 3 和。4 和 5 和。3.设有矩阵a和b 1 求它们的乘积c。2 将矩阵c的右下角3 2子矩阵赋给d。3 ...
Matlab 3 MATLAB程序设计
辽宁工程技术大学上机实验报告。1,已知函数计算。function y f x if x 1 x 0 y x 1elseif x 0 x 1 y 1else x 1 x 2 y x 2end f 1 y ans f 0.5 y ans f 1.5 ans y ans 2 用for end循环语句求 1...
matlab程序设计作业题
1 计算的值。解 输入编程 1.369 2 sin 7 10 pi sqrt 26.48 2.9 结果 3.3097 2 计算下列函数在x 2.43时的函数值并把结果表示成分数形式。解 1 输入编程 format rat x 2.43 y sqrt x sqrt x sqrt x sqrt 2 x ...