matlab中if语句使用方法和应用实例

最新推荐文章于 2025-03-26 20:59:13 发布

九灵猴

最新推荐文章于 2025-03-26 20:59:13 发布

阅读量10w+

点赞数 23

文章标签： matlab

1、我们首先需要了解if语句是一种选择判断语句，可以和for语句结合使用，也可以单独使用，这里创建了3个变量，作出了一个简单判断如果a小于b，b等于a，需要注意的是if语句后面要有end语句结束。运行脚本之后，在命令行窗口可以看到b已经等于a了。如下图所示：

2、介绍一下if嵌套语句，这里作出的判断是如果a小于b,b等于a，如果继续b还大于c，b等于c，如下图所示：

3、介绍if语句中的else语句，else是否则的意思，这里作出的判断是如果a小于b,b等于a，如果满足a小于b条件下，b大于c，b就等于c，否则b=c-a,如下图所示：

博客等级

码龄6年

1
原创

102
点赞

457
收藏

30
粉丝

关注

私信

热门文章

分类专栏

matlab 6篇

最新评论

采用Newmark-β法求解振动方程——MATLAB
站在井底憧憬星空: 代码好多地方有问题，需要修正一下 clc clear m=[1,2,3,4]; m=diag(m); k= [800 -800 0 0; -800 2400 -1600 0; 0 -1600 4800 -3200; 0 0 -3200 8000]; c=0.05*m+0.02*k; f0=100; t1=5; nt=2000; dt=0.01; gama=0.5; beta=0.25; a0=1/ beta/dt/dt; a1=gama/ beta /dt; a2=1/beta/dt; a3=1/2/ beta -1; a4=gama/beta-1; a5=dt/2*(gama/ beta-2); a6=dt*(1-gama); a7=dt*gama; d=zeros(4,nt+1); v=zeros(4,nt+1); a=zeros(4,nt+1); tt=zeros(1,nt+1); for i=2:(nt+1) tt(1,1)=1*dt; tt(1,i)=i*dt; t=(i-1)*dt; if (t<t1) f=[f0*sin(4*pi*t/t1);0;0;0]; else f=[0;0;0;0]; end ke=k+a0*m+a1*c; fe=f+m*(a0*d(:,i-1)+a2*v(:,i-1)+a3*a(:,i-1))+c*(a1*d(:,i-1)+a4*v(:,i-1)+a5*a(:,i-1)); d(:,i)=inv(ke)*fe; a(:,i)=a0*(d(:,i)-d(:,i-1))-a2*v(:,i-1)-a3*a(:,i-1); v(:,i)=v(:,i-1)+a6*a(:,i-1)+a7*a(:,i); end figure; plot(tt(1,:),d(1,:)) ylabel('位移/[m]','fontsize',20); grid on; xlabel('时间/[s]','fontsize',20); title('顶层楼板位移响应时程','fontsize',25);
采用Newmark-β法求解振动方程——MATLAB
站在井底憧憬星空: clc clear m=[1,2,3,4]; m=diag(m); k= [800 -800 0 0; -800 2400 -1600 0; 0 -1600 4800 -3200; 0 0 -3200 8000]; c=0.05*m+0.02*k; f0=100; t1=5; nt=2000; dt=0.01; gama=0.5; beta=0.25; a0=1/ beta/dt/dt; a1=gama/ beta /dt; a2=1/beta/dt; a3=1/2/ beta -1; a4=gama/beta-1; a5=dt/2*(gama/ beta-2); a6=dt*(1-gama); a7=dt*gama; d=zeros(4,nt+1); v=zeros(4,nt+1); a=zeros(4,nt+1); tt=zeros(1,nt+1); for i=2:(nt+1) tt(1,1)=1*dt; tt(1,i)=i*dt; t=(i-1)*dt; if (t<t1) f=[f0*sin(4*pi*t/t1);0;0;0]; else f=[0;0;0;0]; end ke=k+a0*m+a1*c; fe=f+m*(a0*d(:,i-1)+a2*v(:,i-1)+a3*a(:,i-1))+c*(a1*d(:,i-1)+a4*v(:,i-1)+a5*a(:,i-1)); d(:,i)=inv(ke)*fe; a(:,i)=a0*(d(:,i)-d(:,i-1))-a2*v(:,i-1)-a3*a(:,i-1); v(:,i)=v(:,i-1)+a6*a(:,i-1)+a7*a(:,i); end figure; plot(tt(1,:),d(1,:)) ylabel('位移/[m]','fontsize',20); grid on; xlabel('时间/[s]','fontsize',20); title('顶层楼板位移响应时程','fontsize',25);
采用Newmark-β法求解振动方程——MATLAB
m0_74002118: 请问阻尼矩阵c是怎么确定的呢，具体地说是其系数0.05和0.02是怎么算出来的？
采用Newmark-β法求解振动方程——MATLAB
m0_74002118: 请问阻尼矩阵c是怎么确定的呢
matlab中for语句使用方法和应用实例
tangguomanman: i=1：5和i=1：1：5也不会显示中间的结果，会覆盖，需要自己创建一个矩阵存中间数据

大家在看

【教程4＞第7章＞第1节】信道编译码概述

最新文章

目录

评论

被折叠的条评论为什么被折叠?

到【灌水乐园】发言

查看更多评论

添加红包

成就一亿技术人!

hope_wisdom

发出的红包

实付元

使用余额支付

点击重新获取

扫码支付

钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。