我matlab学习之路

离线sailor0000

Never Stop The Beat !!!

UID ：126185

注册：2018-01-31
登录：2018-06-08
发帖：68
等级：仿真一级

20楼发表于: 2018-02-10 10:49:24

顶起，工作中有同事使用matlab建模，他们是研究算法开发。matlab语法简单，画图功能强大，可以让建模人员把精力集中于模型本身而不需要再研究编程语言，这应该是matlab最强大的地方——入门建单，但是可深可浅，交由使用人员决定。

Never Stop The Beat !!!

离线守夜人

UID ：118005

注册：2015-09-16
登录：2025-04-27
发帖：92
等级：论坛版主

21楼发表于: 2018-02-11 19:02:20

回去过年了，没有网，更新会停一停

离线守夜人

UID ：118005

注册：2015-09-16
登录：2025-04-27
发帖：92
等级：论坛版主

22楼发表于: 2018-02-19 18:19:10

微带带通滤波器的特性阻抗计算

复制代码

syms w Z0 Z1 Z2
syms g0 g1 g2
syms Zed Zod Zes Zos
w=0.1 %带宽
Z0=50 %特性组抗
g0=1.000
g1=1.5963
g2=1.0967
Z1=(pi*w/(2*g0*g1))^(1/2)
Z2=(pi*w)/(2*(g1*g2)^(1/2))
Zed=Z0*[1+Z1+Z1^2] %单数偶
Zod=Z0*[1-Z1+Z1^2] %单数奇
Zes=Z0*[1+Z2+Z2^2] %双数偶
Zos=Z0*[1-Z2+Z2^2] %双数奇

离线守夜人

UID ：118005

注册：2015-09-16
登录：2025-04-27
发帖：92
等级：论坛版主

23楼发表于: 2018-02-19 21:45:39

过年放烟花了

复制代码

function firework
% 烟花烟花满天飞
OldHandle = findobj( 'Type', 'figure', 'Tag', 'FireWork' ) ;
if ishandle(OldHandle)
close(OldHandle) ;
end
% 图形窗口初始化
fig = figure('units','normalized','position',[0.1 0.1 0.8 0.8],...
'menubar','none','name','烟花满天飞',...
'numbertitle','off','color',[0 0 0],'tag','FireWork');
% 烟花炸开前的初始位置
h0 = line(0,0,0,'marker','.',...
'markersize',18,...
'MarkerEdgeColor',[1 1 1],...
'MarkerFaceColor',[1 1 1],...
'EraseMode' , 'xor',...
'color',[0 0 0]);
% 设置坐标系显示属性
axis equal
axis([-50 50 -50 50 0 100])
axis off
view(-42,22)
% 设置参数
rate = 1:-0.01:0; % 颜色衰减率
v0 = 45; % 烟花头初始速度
g = -9.8; % 重力加速度
usedtime = -v0/g; % 烟花头爆炸前所经历时间
zs = v0*usedtime+0.5*g*usedtime^2; %烟花头爆炸前达到的最高高度
theta = rand(250,1)*2*pi; % 各粒子速度的方位角
phi = rand(250,1)*2*pi-pi; % 各粒子速度的仰角
age = 20; % 粒子生存期
% 常用颜色矩阵
colormat = [1.0 0.5 0.5
1.0 0.75 0.5
1.0 1.0 0.5
0.75 1.0 0.5
0.5 1.0 0.5
0.5 1.0 0.75
0.5 1.0 1.0
0.5 0.75 1.0
0.5 0.5 1.0
0.75 0.5 1.0
1.0 0.5 1.0
1.0 0.5 0.75];
% 随机产生各粒子对应的颜色序号
colorid = randsample(12,250,true);
% 粒子对应的颜色矩阵
colormat = colormat(colorid,:);
% 粒子颜色与背景色（夜色）的距离
colordist = sqrt(sum(colormat.^2,2));
v1 = 20; % 粒子的初始速度
k = 1; % 颜色衰减率初始序号
timerA = timer('TimerFcn',@TimerFcnA,...
'executionmode','fixedspacing','Period',0.001);
start(timerA);
h = getappdata(gcf,'HandleParticle');
timerB = timer('TimerFcn',{@TimerFcnB},...
'executionmode','fixedspacing','period',0.001);
%--------------------------------------------------------------------------
% 定时器回调函数（烟花头）
%--------------------------------------------------------------------------
function TimerFcnA(timerA,event)
ta = get(timerA,'TasksExecuted')*0.1;
if ta <= usedtime
z = v0*ta+0.5*g*ta^2;
set(h0,'zdata',z,...
'MarkerEdgeColor',[1 1 1],...
'MarkerFaceColor',[1 1 1])
drawnow
%pause(0.01)
else
delete(h0)
stop(timerA);
x0 = zeros(2,250);
y0 = zeros(2,250);
z0 = zs*ones(2,250);
h = line(x0,y0,z0,'marker','h',...
'markersize',12,...
'MarkerEdgeColor',[1 1 1],...
'MarkerFaceColor',[1 1 1],...
'EraseMode' , 'xor',...
'color',[0 0 0]);
setappdata(gcf,'HandleParticle',h);
start(timerB);
end
end
%--------------------------------------------------------------------------
% 定时器回调函数（粒子）
%--------------------------------------------------------------------------
function TimerFcnB(timerB,event)
tb = get(timerB,'TasksExecuted')*0.15;
if age>0 && any(colordist>=0.05)
colormat = colormat*rate(k);
colordist = sqrt(sum(colormat.^2,2));
for i = 1:250
xi = v1*cos(phi(i))*cos(theta(i))*tb;
yi = v1*cos(phi(i))*sin(theta(i))*tb;
zi = zs+v1*sin(phi(i))*tb+0.5*g*tb^2;
set(h(i),'XData',xi,'YData',yi,'ZData',zi,...
'MarkerEdgeColor',colormat(i,:),...
'MarkerFaceColor',colormat(i,:))
end
drawnow
age = age-0.1;
k = k+1;
else
stop(timerB);
delete(h)
end
end
%-------------------
end

离线守夜人

UID ：118005

注册：2015-09-16
登录：2025-04-27
发帖：92
等级：论坛版主

24楼发表于: 2018-02-20 20:39:24

下雪了

复制代码

A=imread('bg.jpg');%背景图像
h=figure;
imshow(A); %显示背景图像
axes
axis([0 1 0 1]); %设置坐标范围，并隐藏坐标
axis off;
set(gcf,'color','k'); %设置背景颜色为黑色
%设置显示的雪花数目
N=20;
handles=zeros(1,N);
x=rand(2,N);
% fontsize(大) step 步长
fontsize=floor(rand(1,N)*17+32); %设置字号，32-48之间
new_handles_N=0;
%% 绘图部分
for i=1:N
handles(i)=text(x(1,i),x(2,i),'*','fontsize',fontsize(i),'color','w'); %handles 存储每一个“雪花”的句柄
end
while 1
if ~ishandle(h)
return
end
for i=1:N
temp=get(handles(i),'position');
step=get(handles(i),'fontsize')/48*0.05; % 不同大小的雪花，速度不同。
if temp(2)<0 %判断是否飘出坐标范围
new_handles_N=new_handles_N+1;
new_handles(new_handles_N)= copyobj(handles(i),gca);
if new_handles_N==500 %扫雪了~\(≧▽≦)/~啦啦啦
delete(new_handles);
new_handles_N=0;
end
temp(1)=rand(1);
temp(2)=1;
else
temp(1)=temp(1)+rand(1)*0.1-0.05;
temp(2)=temp(2)-step; %速度
end
set(handles(i),'position',temp,'rotation',get(handles(i),'rotation')+5);
end
pause(.2)
end