zhendongjisuan.m

内燃机转子仿真

clear all;
clc;
close all;
M=[];
K=[];
C=[];
f0=;
force=input('the number of the force:')
switch force
    case 1;
        f=f0*eye(n);%阶跃
    case 2;
        f=[];%等倾上升阶跃
    case 3; 
        t1=;
        ft1=[];%指数渐进阶跃
    case 4;
        f=[];%矩形脉冲
    case 5;  
        f=[];%梯形脉冲
    case 6;
        f=[];%三角形脉冲
    case 7;  
        f=[];%后峰锯齿形脉冲
    case 8;
        f=[];%前峰锯齿形脉冲
    case 9;  
        f=[];%半圆脉冲
    case 10;
        f=[];%倒争弦脉冲
    case 11;  
        f=[];%1/4余弦脉冲
    case 12;
        f=[];%指数衰减
    case 13;  
        f=[];%一周正弦
    case 14;
        f=[];%一周余弦
    case 15;  
        f=[];%一周矩形波
    case 16;
        f=[];%直线上升
    case 17;  
        f=[];%抛物线上升
n=;
[m,c,k]=arryf(M,C,K,n);