实验5综合应用程序设计.doc

实验5. 综合应用程序设计 实验目的 1.理解命令.m文件与函数.m文件的异同； 2.掌握定义和调用matlab函数的方法。 实验环境 计算机 matlab7.1 实验说明 熟悉matlab7.1运行环境，正确操作，认真观察； 实验学时：2学时； 学会使用help； 编写程序前先写出程序的流程图框架； 自主编写调试程序直到正确为止，必要时参考相关资料，； 保存整理实验结果，提交实验报告。 实验内容 1. 设两个系统各自有5个冗余件可供选择，可靠性与成本如下表（数据来自上海大众4S店）所示。则对成本与可靠性的集成研究可化为以下两个系统中四个冗余件的优化组合。（注：尽量不用循环，用矩阵结构，保存文件名为s5_学号_01.m） 图 调角器系统与电机系统可靠性的并串联结构 表 调角器系统与电机系统可靠度及成本数据 调角器系统A 电机系统B R1i C1（=R1i×65.470） R2j C2（=R2j×119.281） 1 R11 =0.961 R21 = 0.951 2 R12 =0.920 R22 0.820 3 R13 =0.876 R23 0.788 4 R14 =0.792 R24 0.712 5 R15 =0.723 R25 0.657 目标： Min C = C1 + C2 约束条件： C1 = 65.470× (i取1,2,3,4,5中任两个数) C2 = 119.281× (j 取1,2,3,4,5中任两个数) R1 ＝1－ (i取1,2,3,4,5中任两个数) R2 ＝1－ (j取1,2,3,4,5中任两个数) R1 ≥ 0.9825 ， R2 ≥ 0.9819 实验程序及结果 （1）程序： R1i=[0.961 0.920 0.876 0.792 0.723]; R2j=[0.951 0.820 0.788 0.712 0.657]; A=[R1i(1)+R1i;R1i(2)+R1i;R1i(3)+R1i;R1i(4)+R1i;R1i(5)+R1i];%两个R1i相加 B=[R2j(1)+R2j;R2j(2)+R2j;R2j(3)+R2j;R2j(4)+R2j;R2j(5)+R2j];%两个R2j相加 C=(1-R1i)*(1-R1i);%两个(1-R1i)相乘 D=(1-R2j)*(1-R2j); %两个(1-R2j)相乘 R1=1-C;%R1的值 R2=1-D;%R2的值 C1=65.470*(A.*(R1=0.9825).*(R2=0.9819)); %在满足（R1=0.9825）（R2=0.9819）的条件下求C1，C2 C2=119.281*(B.*(R1=0.9825).*(R2=0.9819)); Cm=C1+C2; C=Cm-diag(diag(Cm)) %令Cm中的对角线上元素值为0，即排除选择的是同一个冗余件的情况。 %此处Cm就是C=C1+C2的所有可能情况 C=C(find(C~=0));%把Cm中的非零元素提取出来，组成新的矩阵Cm Min=min(C) （2）结果： C = 0 334.3957 327.6980 313.1332 302.0553 334.3957 0 0 0 0 327.6980 0 0 0 0 313.1332 0 0 0 0 302.0553 0 0 0 0 Min = 302.0553 2．定义一个函数文件,功能是采用重叠相加法完成无限长序列与有限长序列的线性卷积，然后用该函数文件计算序列x[k]=2k+1,0≤k≤18与序列h[k]={1,2,3,4;k=0,1,2,3}的线性卷积结果。 （注：可能用到的函数conv,保存文件名为s5_学号_02.m） （1）函数文件 %% 重叠相加法， 实现长短序列的线性卷积 function y=overlap_add(x,h,M) % x为较长的输入序列，h为较短的系统函数，M为分段大小 N=length(h); % 序列h(n)的长度 Lx=length(x); % 序列x(n)的长度 if NM % 算