Posts on 极客Vinkey

Cisco AnyConnect优雅脚本与openconnect进阶教程

Fri, 07 Aug 2026 14:12:58 +0000

https://github.com/ma-wenqian/vpn-playbook

概述

一般情况下，能用 Cisco AnyConnect（现称 Cisco Secure Client）的场景，基本都可以用 OpenConnect 替代。从便捷到进阶，本仓库提供了对应的脚本，辅助自动登录学校（或类似配置的）VPN。

目前 Cisco 基本都启用了双重认证：先输入 VPN 地址，然后输入账号、密码，接着输入二次验证码——也就是常见存储在 Google Authenticator、iCloud 钥匙串、Microsoft Authenticator 中的那个 6 位验证码。这个验证码本质上是基于一个固定的 Base32 密钥和当前时间戳计算得到的（TOTP），并没有什么"黑盒"逻辑。因此，只要在绑定 2FA 时提取到这个固定密钥，就可以自己算出验证码，从而实现全程自动登录：

Unix 平台：使用 oathtool 工具，OTP=$(oathtool --totp -b "$SECRET")
Windows 平台：使用 totp-go 命令行工具，totp-go "$SECRET"

⚠️ 请像保管密码一样保管这个密钥！ 任何拿到这个密钥的人都可以无限期地为你的账号生成有效的二次验证码——它不像单次验证码那样会失效。切勿把填入了真实用户名/密码/密钥的配置文件提交到（哪怕是私有的）代码仓库。本仓库中的 config.bat / config.sh 都只是占位模板，请在本地填入真实信息后不要把这部分改动提交到 git。

仓库结构

.
├── Cisco/
│   ├── unix/
│   │   ├── cisco-vpn.sh          # 基于 oathtool + `vpn -s`（脚本化 stdin）实现 connect/disconnect
│   │   ├── config.sh             # 你的 VPN 地址/账号/密码/TOTP 密钥（模板）
│   │   └── install.sh            # 一键安装脚本（curl | bash）
│   └── win/
│       ├── cisco-vpn.bat        # 双击运行的入口脚本，用 totp-go 生成验证码并直接写入 vpncli.exe 的脚本化输入
│       ├── config.bat           # 你的 VPN 地址/账号/密码/TOTP 密钥（模板，变量名与 unix 端一致）
│       ├── icon/                # 桌面快捷方式用到的 connect.ico / disconnect.ico
│       └── install.ps1          # 一键安装脚本，顺带安装 totp-go，还能自动创建桌面快捷方式（irm | iex）
└── Openconnect/
    ├── openconnect-setup.sh       # OpenConnect + microsocks 的 Docker 部署（更轻量）
    └── openconnect-gost-setup.sh  # OpenConnect + gost 的 Docker 部署（功能更全）

移动端

iOS：无法做到完全自动化，顶多是借助自动填充功能填入验证码。
Android：本仓库暂未研究。

Cisco（桌面端，仅命令行）

不要启动 GUI 界面，全部使用命令行实现自动化操作。

Windows快速切换简体繁体拼音输入-Ctrl+Shift+F

Mon, 27 Apr 2026 17:37:23 +0000

新发现：Ctrl+Shift+F

今天不知道电脑怎么突然就变成繁体中文输入法了，我想着也没有下载啥软件、修改系统设置，在网上搜了一下，发现应该是误触了快捷键，那么相当于也知道了一个快速切换的方法：Ctrl+Shift+F

对应系统设置：时间和语言→语言和地区→Microsoft拼音→按键最下方，如果按了没反应，大概率是和其他软件的快捷键冲突了，回到桌面内就可以了。

以后就可以快速通过Ctrl+Shift+F切换简体和繁体中文了~

关于微软输入法变成繁体问题

VSCode+Latex Workshop环境配置

Thu, 04 Dec 2025 17:52:14 +0000

本文总结了在windows 和wsl/Ubuntu平台部署latex环境，在vscode中写文件和编译的过程。

在线版本的overleaf可以免去配置环境，方便进行多人协作。但是在前期写作阶段，我比较倾向于本地环境，方便修改图片，不用反复上传。可以在初稿完成后，上传至overleaf协作修改。

一共两个步骤：

配置对应系统的tex
vscode安装Latex Workshop插件

即可实现一键编译。编译速度上应该是mac>linux>wsl>windows，所以在windows中可以额外配置wsl的latex环境。

安装Tex

这个环节可能要很长时间，可以选择从国内镜像站下载。安装的过程中可以同步进行下一步操作。

Windows

安装texlive，下载install-tl-windows.exe，一键安装。安装过程需要一段时间，可以放在后台等。同样的网络下，wsl中安装速度比windows快了很多。

Ubuntu

根据texlive网页说明，在terminal中运行，

cd /tmp # working directory of your choice
Download: wget https://mirror.ctan.org/systems/texlive/tlnet/install-tl-unx.tar.gz
           or: curl -L -o install-tl-unx.tar.gz https://mirror.ctan.org/systems/texlive/tlnet/install-tl-unx.tar.gz

zcat < install-tl-unx.tar.gz | tar xf - # note final - on that command line
cd install-tl-2*
perl ./install-tl --no-interaction # as root or with writable destination
# may take several hours to run

安装完成后需要添加路径/usr/local/texlive/YYYY/bin/PLATFORM 到 PATH，以实际屏幕显示的路径进行设置。

我的科研绘图调色板

Thu, 27 Nov 2025 18:46:45 +0000

本文收集了一些科研图配色，并提供颜色效果预览和hex、rgb数值。

🌈 我的科研绘图调色板

双色

色块预览	线条预览	HEX 值	RGB 值	颜色名称
		`#b22222`	`(178, 34, 34)`
		`#000080`	`(0, 0, 128)`

色块预览	线条预览	HEX 值	RGB 值	颜色名称
		`#f5a889`	`(245, 168, 137)`
		`#acd6ec`	`(172, 214, 236)`

色块预览	线条预览	HEX 值	RGB 值	颜色名称
		`#148f28`	`(20, 143, 40)`
		`#ea71ae`	`(234, 113, 174)`

三色

马卡龙配色

端口转发

Tue, 04 Nov 2025 17:24:27 +0000

本文总结windows、Ubuntu平台端口转发和端口占用查询方法，mac的可以查看网站其他帖子。

Windows

使用netsh interface portproxy命令进行管理（添加，删除）即可，快捷方便。

创建端口转发

v4tov4 代表用的是ipv4地址 listenaddress 监听地址 listenport 监听端口号 connectaddress 转发的目的地地址 connectport 转发目的地端口号

netsh interface portproxy add v4tov4 listenaddress=127.0.0.1 listenport=1085 connectaddress=100.122.123.54 connectport=1085

查看端口使用

netsh interface portproxy show all

Ubuntu

如何将Cisco Anyconnect共享给局域网使用

Wed, 29 Oct 2025 15:18:25 +0000

cisco默认split tunnel功能是关闭的，所以开启后默认拦截所有路由且包括本地局域网，导致无法访问路由器管理页面、打印机、部署在本地的其他其他服务等，这就很不方便。

ios端的anyconnect支持connect on demand可以手动把局域网或者特定ip/域名加入白名单，其中在mac端可以下载ipad上的anyconnect同样可以实现对应功能，不要用mac端的cisco（这个也是拦截所有的）。所以一种曲线救国的方式就是在mac平台开启一个本地代理server连接让局域网使用，这个就比较简单了，但是本人现在的主力设备是window，经过探索发现下面的方式。

所需软件

WSL
3proxy
frp

WSL配置

wsl，并设置为mirror网络模式，这样wsl中可以访问局域网，跳过cisco阻断所有网络，同时也可以实现在wsl中访问局域网设备，比如安装chrom后，用google-chrome启动浏览器，可以访问本地网络。

WSL中安装3Proxy并开启Socks代理

使用3proxy进行安装，直接在局域网启动就行，安全性没问题。默认开启1080 sock端口配置文件在`/etc/3proxy/conf/3proxy.cfg

nscache 65536
nserver 8.8.8.8
nserver 8.8.4.4

config /conf/3proxy.cfg
monitor /conf/3proxy.cfg

log /logs/3proxy-%y%m%d.log D
rotate 60
counter /count/3proxy.3cf


include /conf/counters
include /conf/bandlimiters

auth none
allow *
socks -i127.0.0.1

admin -p8080 -i127.0.0.1

利用FRP将端口穿透至局域网

frp 进行内网穿透，将上述代理端口映射出去。理论上可以映射给一个公网服务器，这样相当于有了随时可以连接的代理服务器，但不建议，因为存在安全隐患。

frps 服务端

在openwrt/istoreos中安装，当前测试的版本是0.65.0

frpc 客户端

在wsl中安装，我用的是GitHub自动安装脚本https://github.com/stilleshan/frpc 然后根据提示信息修改端口

serverAddr = "10.0.0.1"  
serverPort = 7000  
auth.method = "token"  
auth.token = ""  
  
  
[[proxies|proxies]]  
name = "cisco_vpn"  
type = "tcp"  
localIP = "127.0.0.1"  
localPort = 1080  
remotePort = 1080

然后sudo systemctl restart frpc重启服务

记录一下脾胃虚弱中药方

Fri, 28 Apr 2023 11:53:30 +0000

症状表现

胃口差，爱吃酸、甜的食物，嗜睡，身体很疲惫，晚上失眠，脸色暗沉……

不好的习惯：在脾胃不好的前提下，早上起的晚，早餐经常吃酸奶（有点凉）+面板+鸡蛋，爱吃水果、鸭肉

生活调理

首先在身体恢复前要控制饮食，生冷酸辣甜的不吃、水果不吃，不能喝冷水，只能喝热水/温水，晚上早点睡觉，缓解压力

中医调理

一共花了¥1,159.36

3月17日，我第一次去的是佛山一个诊所，这边给的是中药颗粒需要自己用开水泡开，因为没有做胃镜所以给了一个治疗胃溃疡的药，共花费¥367。
当时发了一个朋友圈，一位老师建议我去省中医的脾胃科看一下，4月4日去了省中医琶洲分院，挂号费¥10，药费¥318.79, 医生：庄壮锋
- 西药（担心有胃溃疡）：（国）雷贝拉唑钠肠溶胶囊/20mg*14粒/盒，¥27.2
- 中药/7天用量，¥272.34，代煎费2.75*7=¥19.25
4月14日，拿着药方去了校医院开了3天的中药，报销后¥114.83
4月20日，去了省中医大德路总院，我说有效果，现在感觉好多了，医生说这次吃完可以不用来了，药有一点点更新，挂号费¥10，药费¥338.74
- 西药（担心有胃溃疡）：（国）雷贝拉唑钠肠溶胶囊/20mg*14粒/盒，¥27.2
- 中药/7天用量，¥292.24，代煎费2.75*7=¥19.25

药方

药方04/04

药品名称	总量	金额（元）
木香（广木香）	10.00g*7	¥4.55
郁金	15.00g*7	¥12.21
延胡索(醋延胡索/玄胡索/元胡/制延胡索）	15.00g*7	¥20.35
香附（醋/制香附）	10.00g*7	¥4.55
合欢皮	20.00g*7	¥5.78
党参	15.00g*7	¥36.75
广陈皮	5.00g*7	¥55.48
法半夏	15.00g*7	¥34.13
炙甘草	5.00g*7	¥3.55
茯神	20.00g*7	¥17.50
白术	15.00g*7	¥16.67
佛手	15.00g*7	¥24.94
海螵蛸（鱼古）	20.00g*7	¥21.00
浙贝母	10.00g*7	¥14.88

药方04/20

药品名称	总量	金额（元）
木香(广木香)	10.00g*7	¥4.50
郁金	15.00g*7	¥12.21
延胡索(醋延胡索/玄胡索/元胡/制延胡索）	15.00g*7	¥20.35
黄芪(北芪)	20.00g*7	¥24.50
合欢皮	20.00g*7	¥5.78
党参	15.00g*7	¥36.75
广陈皮	5.00g*7	¥55.48
法半夏	15.00g*7	¥34.13
炙甘草	5.00g*7	¥3.55
茯神	20.00g*7	¥17.50
白术	15.00g*7	¥16.67
佛手	15.00g*7	¥24.94
海螵蛸(鱼古)	20.00g*7	¥21.00
浙贝母	10.00g*7	¥14.88

痊愈

在2023年底，我去校医院再次号脉，医生说现在已经没问题了。同时我也在B在更新了一条视频：

Mac查看端口占用和释放端口的教程

Tue, 25 Apr 2023 20:08:57 +0000

FATAL Port 4000 has been used. Try other port instead.

1.查看端口被哪个程序占用

sudo lsof -i tcp:port

示例：

sudo lsof -i tcp:8080 查看8080端口是否被占用

2.看到进程的PID，可以将进程杀死。

sudo kill -9 PID

示例：

sudo kill -9 3210

MacOS教程| Mac上查看端口占用和释放端口的教程

清除iPhone短信app右上角的小红点

Sun, 23 Apr 2023 23:15:02 +0000

最新的的ios系统已经可实现一键全部已读，去掉小红点了。

来源

去年换了手机后，不知道从哪天开始，短信app右上角的数字没办法清零，我很难受。

最后在知乎上找到了一个好办法，归根到底是ios系统bug¹，很古老的未读短信自动隐藏（未加载出来），但又会在右上角显示。

曲折

Apple社区只会给出垃圾解决办法

短信上未读提醒数无法消除 - Apple 社区 https://discussionschinese.apple.com/thread/253437710?answerId=256489744322#256489744322 根本不管用！

正解

第一步：先在“未知发件人”里一直往下拉，拉到最底部，也就是这台手机未知发件人里最早的那条短信。

第二步：点击右上角三个点，“选择信息”，点击左上角“全选”，再点击左下角“未读短信”

第三步：回到“所有短信”，一样一直往下拉，拉到最底部。 ps：直接在“未读短信”里这样操作也行，但有时候会失败，所以最保险的做法还是回到“所有短信”操作

第四步：点击右上角三个点，“选择短信”，点击左下角“全部已读”。然后会出现好像卡死一样，不着急，停留在界面等，一直等到处理完成，就大功告成了！！ ps：等待处理的过程中（点击完“全部已读”后）不要返回主界面，什么都不要点，就等着它处理完，不然可能会导致失败，当然失败的话按着步骤重新来一遍就行。

Iphone短信全部已读了但是还是有未读数量怎么办啊? - 知乎作者：做梦去了趟比邻星 ↩︎

问卷调查数据采用Excel处理生成数据看板

Fri, 14 Apr 2023 10:52:57 +0000

月总结-2023

Thu, 30 Mar 2023 09:21:04 +0000

3月很ok

3月收获了好消息，而且现在科研道路越走越坦然和自信了，表现在grasshopper建模速度突飞猛进，一个下午的工作量堪比研一一周的，我真的佩服我自己哈哈哈哈！

面对实验的意外、失败这些，我也没有那么消极了，抱着早发现早治疗、早发现早换赛道、小插曲也是工作量、非必要不加班、多运动的科研宗旨。而且新学期我也搬了工位，记得胡老师说过定期换换位置，有助于提升科研积极性，我把实验室的一个绿植放在旁边，感觉很舒服。

又用我的兴趣搭建了实验室wiki网页（基于Docusaurus）、在线预约系统（基于MRBS）、内网网盘（Alist），买了一个10年域名（要报销哈哈哈），估计够用了。

有空，游泳健身了解一下！

4月是养精蓄锐

这个月放慢了实验速度，更多时间在分析实验逻辑和模拟上来，我觉得这样是对的，埋头做实验久了需要跳出来，针对前一阶段的内容提炼分析，每一次的实验都应对后续实验产生影响，这里指实验细节和逻辑设置，而不是数据影响。生活上，这个月一直在吃中药调理肠胃，现在感觉胃口好多了，在此基础上我健身更规律了，完成了健身房的打卡任务，可以领83块钱了，而且最近肌肉有明显增加，看起来更壮了。 5月主线任务：科研+雅思+健身

Matlab绘图代码共享-硕士阶段

Sun, 12 Mar 2023 14:46:54 +0000

拉压弯曲线图及柱状图2022/12

更新于2023.03.13

{% folding child:codeblock open:flase color:yellow 点击展开代码详情 %}

%% 基础信息
% 反复修改的代码简直屎山！
clc;clear
close all


% color1=[77,133,189]/255;
% color2=[247,144,61]/255;
% color3=[89,169,90]/255;

% 马卡龙配色
color1=[147 213 220]/255;
color2=[242 202 201]/255;
color3=[198 223 200]/255;

% color1=[255 230 215]/255;
% color2=[225 240 219]/255;
% color3=[222 234 248]/255;

% color1=[240   181   125]/255;
% color2=[211   225   174]/255;
% color3=[113   171   182]/255;
% 
% color1=[   76   139   192]/255;
% color2=[  232    49    51]/255;
% color3=[   94   183    91]/255;

% 大论文使用的配色
% color1=[  110 158 206]/255;
% color2=[  118 186 128]/255;
% color3=[  230 146 143]/255;

FontSize=13;
color3s={color1;color2;color3};
colors={color1;color1;color1;color2;color2;color2;color3;color3;color3};
colors9=[  
   199    92   100
  240   181   125
  211   225   174
  113   171   182
   75    90   161
   76   139   192
  232    49    51
  163    96   173
   94   183    91]/255;


ki=[ 1	0	0
    0.75	0.25	0
    0.5	0.5	0
    0.25	0.75	0
    0	1	0
    0	0.75	0.25
    0	0.5	0.5
    0	0.25	0.75
    0	0	1];
bar_colors=ki*[color1;color3;color2];

legends=["v=240mm/s  T=210°C  h=0.16mm"
    "v=250mm/s  T=215°C  h=0.20mm"
    "v=260mm/s  T=220°C  h=0.24mm"
    "v=240mm/s  T=215°C  h=0.24mm"
    "v=250mm/s  T=220°C  h=0.16mm"
    "v=260mm/s  T=210°C  h=0.20mm"
    "v=240mm/s  T=220°C  h=0.20mm"
    "v=250mm/s  T=210°C  h=0.24mm"
    "v=260mm/s  T=215°C  h=0.16mm"]';
legends9=[legends;repmat("",5,9)];
%%

hexc=[
'32aeec'
'E3863c'
'E9262e'];

r=hex2dec(hexc(:,[1,2]));
g=hex2dec(hexc(:,[3 4]));
b=hex2dec(hexc(:,[5 6]));
rgb=[r,g,b];
%% 拉伸


Tensile_x=cell(9,6);       % read strain to Tensile_x
Tensile_y=cell(9,6);       % read stress to Tensile_y

% 读取原始数据并储存

% 设置导入选项并导入数据
opts = delimitedTextImportOptions("NumVariables", 8);
% 指定范围和分隔符
opts.DataLines = [8, Inf];
opts.Delimiter = "\t";
% 指定列名称和类型
opts.VariableNames = ["Load", "LoadY", "Position", "Displacement", "Extension", "ExtensionY", "StepIndex", "Time"];
opts.VariableTypes = ["double", "double", "double", "double", "double", "double", "double", "double"];
% 指定文件级属性
opts.ExtraColumnsRule = "ignore";
opts.EmptyLineRule = "read";
T_L0=50;              % 标距L0
T_S0=4*13;     % 横截面积S0
for Gi=1:9
    for ti=1:6
        filesname=sprintf("/Users/mawenqian/Documents/硕士/VSL/科研/小论文1/基本试件打印参数正交试验/G%d/N%d.txt",Gi,ti);
        readdata = readtable(filesname, opts);
        % 转换为输出类型
        readdata = table2array(readdata);
        Tensile_x{Gi,ti}=readdata(:,5)/T_L0*100;   % 应变 6sh/L^2
        Tensile_y{Gi,ti}=readdata(:,1)/T_S0;        % 应力 3Fl/2bh^2
        clear readdata filesname
    end
end
clear opts
save Tensile_raw_data Tensile_x Tensile_y
%%

% 通过下方原始图像人工判断获得，亦可绘制后展示。
Tensile_start_line=[
    133	137	137	144	145	139
    151	151	147	157	160	155
    176	178	177	179	191	174
168	205	143	194	162	192
157	162	163	154	222	170
196	242	175	206	246	174
    143	141	145	130	138	146
    133	129	135	126	130	127
145	150	163	145	138	142];
Tensile_end_line=[
    2022	1764	3568	2073	2596	3338
    2222	1912	2363	3740	1927	1832
    2613	1821	2390	2353	2041	2040
4000	3000	2449	2667	2557	1937
2000	1943	2000	1776	2000	2000
1719	1719	1719	1719	1719	1719
  2312	1975	2537	2161	1730	1979
1407	1238	1345	1260	1593	1684
1656	2419	2000	2000	2000	2000];
Tensile_maxy_line=[
740	709	744	730	755	716
1066	1103	1064	1167	1058	1157
834	709	735	779	779	723
1377	1441	1622	1512	1698	1399
1043	1074	1069	1140	1072	1036
1214	1142	1164	1214	1142	1164
1178	1099	1213	1252	1191	1182
1097	1091	1081	1035	1100	1063
853	903	920	896	891	922];

% ttt=1;
% for Gi=[ 7 8]
% for ti=1:6
%     subplot(4,6,ttt)
% 
% 
% x=Tensile_x{Gi,ti};
% y=Tensile_y{Gi,ti};
% [~,Tensile_maxy_line(Gi,ti)]=max(y);
% plot(y,MarkerIndices=Tensile_end_line(Gi,ti),MarkerSize=6,Marker="*",MarkerFaceColor='r')
% title(sprintf("%d-%d",Gi,ti))
%     ttt=ttt+1;
% end
% 
% end

%%
% % 绘制原始数据图像总图3x3
% figure()
% set(gcf,'unit','centimeters','Position',[5 5 30 30]);
% for Gi=1:9
%     subplot(3,3,Gi)
%     hold on
%     for ti=1:6
%         x=Tensile_x{Gi,ti}(Tensile_start_line(Gi,ti):Tensile_end_line(Gi,ti));x=x-x(1);
%         y=Tensile_y{Gi,ti}(Tensile_start_line(Gi,ti):Tensile_end_line(Gi,ti));y=y-y(1);
% Tensile_x{Gi,ti}=x;
% Tensile_y{Gi,ti}=y;
% plot(Tensile_y{Gi,ti},LineWidth=1.5)
% %         plot(Tensile_x{Gi,ti},Tensile_y{Gi,ti},LineWidth=1.5,Marker='o',MarkerEdgeColor='red',MarkerIndices=[Tensile_start_line(Gi,ti)])
%         %         plot(Tensile_y{Gi,ti},LineWidth=1.5)
%     end
%     xlabel("Strain (%)")
%     ylabel("Stress (MPa)")
% %     xlim([0 0.1])
%     % ylim([0,70])
%     % legend(Location="northwest",Box="off");
%     box('on');
%     title(['Tensile ',num2str(Gi)])
%     set(gca,'FontName','Arial','FontSize',FontSize)
% end
%%
%
% % 分组绘制原始图像9张单独图（应力）
% for Gi=[1 2 3 7 8]
%     figure()
%     set(gcf,'unit','centimeters','Position',[5 5 40 28]);
%     for ti=1:6
%         subplot(2,3,ti)
% %         [~,maxy]=max(Tensile_y{Gi,ti}(1:Tensile_end_line(Gi,ti)));
% %         Tensil_maxy_line(Gi,ti)=maxy;
%         plot(Tensile_y{Gi,ti},LineWidth=1.5,Marker='o',MarkerEdgeColor='red',MarkerIndices=[Tensile_start_line(Gi,ti),Tensil_maxy_line(Gi,ti),Tensile_end_line(Gi,ti)])
%         xlabel("N")
%         ylabel("Stress (MPa)")
%         % xlim([0 10])
%     % ylim([0,70])
% %     legend(Location="northwest",Box="off");
%     box('on');
%
%     title(['Tensile ',num2str(Gi), '-',num2str(ti)])
%     set(gca,'FontName','Arial','FontSize',FontSize)
%     end
% end

% 绘制受拉总图
%%
Tensile_maxy=zeros(9,6);

clear fh1
fh1=figure();
set(gcf,'unit','centimeters','Position',[15 15 15 15]);
axes('parent',fh1,'unit','centimeters','position',[1.5 1.5 12.5 12.5]);
hold on
LineStyles=repmat(["-",":","--"],1,3);
for Gi=1:9
    for ti=1:6
        x=Tensile_x{Gi,ti}(Tensile_start_line(Gi,ti):Tensile_end_line(Gi,ti));x=x-x(1);
        y=Tensile_y{Gi,ti}(Tensile_start_line(Gi,ti):Tensile_end_line(Gi,ti));y=y-y(1);
        Tensile_maxy(Gi,ti)=max(y);
        plot(x,y,LineWidth=1.5,Color=colors{Gi},LineStyle=LineStyles(Gi))
        clear x y
    end
    
    xlabel("Strain (%)")
    ylabel("Stress (MPa)")
    xlim([0 8])
    ylim([0,60])
 legend(legends9,Location="northeast",Box="off");
    box('on');
 
    set(gca,'FontName','Arial','FontSize',FontSize)  
end



%% 压缩1屈服强度
compressive1_x=cell(9,6);       % read strain to compressive1_x
compressive1_y=cell(9,6);       % read stress to compressive1_y
Com1_L0=25.4;                   % 标距L0
Com1_S0=12.7*12.7*pi*0.25;      % 横截面积S0


% 读取原始数据并储存
% 设置导入选项并导入数据
opts = spreadsheetImportOptions("NumVariables", 2);
% 指定工作表和范围
opts.Sheet = "Sheet1";
opts.DataRange = "A2:B2000";
% 指定列名称和类型
opts.VariableNames = ["LoadValue", "PositionValue", "PlayTime", "ExtendValue"];
opts.VariableTypes = ["double", "double", "double", "double"];
for Gi=19:27
    for ti=1:6
        filesname=sprintf("/Users/mawenqian/Documents/硕士/VSL/科研/小论文1/基本试件打印参数正交试验/G%d/N%d.xls",Gi,ti);
        readdata = readtable(filesname, opts);
        % 转换为输出类型
        readdata = table2array(readdata);
        compressive1_x{Gi-18,ti}=readdata(:,2)/Com1_L0*100;  % 应变
        compressive1_y{Gi-18,ti}=readdata(:,1)/Com1_S0;      % 应力
        clear readdata filesname
    end
end
clear opts
% 通过下方原始图像人工判断获得，亦可绘制后展示。
compressive1_start_line=[
    73	78	74	77	12	1
    77	87	79	78	75	77
    82	82	80	82	86	83
    369	70	64	72	72	74
    71	67	66	63	64	66
    70	70	68	67	69	69
    62	6	6	5	4	4
    3	3	3	181	1	7
    65	12	5	4	5	64];
compressive1_end_line=[
    471	413	432	390	303	322
    422	432	413	380	380	378
    471	442	480	390	303	322
    827	529	493	497	503	482
    511	526	453	475	496	496
    513	509	528	486	473	494
    528	516	475	414	463	470
    463	454	519	605	463	457
    619	497	496	365	420	478];
compressive1_maxy_line=[310	307	318	314	234	219
    354	377	344	374	352	345
    321	325	333	324	303	322
    728	412	381	403	404	404
    357	350	340	344	345	347
    349	330	332	341	329	339
    331	295	293	279	286	314
    306	285	300	499	313	293
    436	274	300	289	298	323];
% % 绘制原始数据图像总图3x3
% figure()
% set(gcf,'unit','centimeters','Position',[5 5 30 30]);
% for Gi=19:27
%     subplot(3,3,Gi-18)
%     hold on
%     for ti=1:6
%         [~,maxy]=max(compressive1_y{Gi-18,ti}(1:compressive1_end_line(Gi-18,ti)));
%         plot(compressive1_x{Gi-18,ti},compressive1_y{Gi-18,ti},LineWidth=1.5,Marker='o',MarkerEdgeColor='red',MarkerIndices=[maxy])
%         %         plot(compressive1_y{Gi-18,ti},LineWidth=1.5)
%     end
%     xlabel("Strain (%)")
%     ylabel("Stress (MPa)")
%     % xlim([0 10])
%     % ylim([0,70])
%     %     legend(Location="northwest",Box="off");
%     box('on');
%     title(['compressive1 ',num2str(Gi-18)])
%     set(gca,'FontName','Arial','FontSize',FontSize)
% end
% % 分组绘制原始图像9张单独图（应力）
% for Gi=27:-1:19
%     figure()
%     set(gcf,'unit','centimeters','Position',[5 5 40 28]);
%     for ti=1:6
%         subplot(2,3,ti)
%         [~,maxy]=max(compressive1_y{Gi-18,ti}(1:compressive1_end_line(Gi-18,ti)));
% compressive1_maxy_line(Gi,ti)=maxy;
%         plot(compressive1_y{Gi-18,ti},LineWidth=1.5,Marker='o',MarkerEdgeColor='red',MarkerIndices=[compressive1_start_line(Gi-18,ti),compressive1_maxy_line(Gi,ti),compressive1_end_line(Gi-18,ti)])
%         xlabel("N")
%         ylabel("Stress (MPa)")
%         % xlim([0 10])
%         % ylim([0,70])
%         %     legend(Location="northwest",Box="off");
%         box('on');
%
%         title(['compressive1 ',num2str(Gi-18), '-',num2str(ti)])
%         set(gca,'FontName','Arial','FontSize',FontSize)
%     end
% end



% 绘制受压总图
compressive1_maxy=zeros(9,6);

% figure()
% set(gcf,'unit','centimeters','Position',[5 5 30 30]);
% for Gi=19:27
%     Gi=Gi-18;
%     subplot(3,3,Gi)
% 
%     hold on
%     for ti=1:6
%         compressive1_maxy(Gi,ti)=compressive1_y{Gi,ti}(compressive1_maxy_line(Gi,ti))-compressive1_x{Gi,ti}(compressive1_start_line(Gi,ti));
%         x=compressive1_x{Gi,ti}(compressive1_start_line(Gi,ti):compressive1_end_line(Gi,ti));x=x-x(1);
%         y=compressive1_y{Gi,ti}(compressive1_start_line(Gi,ti):compressive1_end_line(Gi,ti));y=y-y(1);
%         compressive1_maxy(Gi,ti)=max(y);
%         plot(x,y,LineWidth=1.5,Color=colors{Gi})
%     end
%     xlabel("Strain (%)")
%     ylabel("Stress (MPa)")
%     xlim([0 8])
%     ylim([0,60])
% 
%     legend(legends(Gi),Location="north",Box="off");
%     box('on');
%     title(['Compressive ',num2str(Gi)])
%     set(gca,'FontName','Arial','FontSize',FontSize)
% end

clear fh1
fh1=figure();
set(gcf,'unit','centimeters','Position',[15 15 15 15]);
axes('parent',fh1,'unit','centimeters','position',[1.5 1.5 12.5 12.5]);
hold on
LineStyles=repmat(["-",":","--"],1,3);
for Gi=19:27
    Gi=Gi-18;
    hold on
    for ti=1:6
        x=compressive1_x{Gi,ti}(compressive1_start_line(Gi,ti):compressive1_end_line(Gi,ti));x=x-x(1);
        y=compressive1_y{Gi,ti}(compressive1_start_line(Gi,ti):compressive1_end_line(Gi,ti));y=y-y(1);
        compressive1_maxy(Gi,ti)=max(y);
        plot(x,y,LineWidth=1.5,Color=colors{Gi},LineStyle=LineStyles(Gi))
    end
    xlabel("Strain (%)")
    ylabel("Stress (MPa)")
    xlim([0 8])
    ylim([0,60])
 legend(legends9,Location="southeast",Box="off");
    box('on');
%     title('Compressive')
    set(gca,'FontName','Arial','FontSize',FontSize)
end
%% 弯曲


Flexural_x=cell(9,6);       % read strain to Flexural_x
Flexural_y=cell(9,6);       % read stress to Flexural_y

% 读取原始数据并储存

% 设置导入选项并导入数据
opts = delimitedTextImportOptions("NumVariables", 8);
% 指定范围和分隔符
opts.DataLines = [8, Inf];
opts.Delimiter = "\t";
% 指定列名称和类型
opts.VariableNames = ["Load", "LoadY", "Position", "Displacement", "Extension", "ExtensionY", "StepIndex", "Time"];
opts.VariableTypes = ["double", "double", "double", "double", "double", "double", "double", "double"];
% 指定文件级属性
opts.ExtraColumnsRule = "ignore";
opts.EmptyLineRule = "read";
L0=6*3.2/52^2;              % 标距L0
S0=3*52/(2*12.7*3.2^2);     % 横截面积S0
for Gi=10:18
    for ti=1:6
        filesname=sprintf("/Users/mawenqian/Documents/硕士/VSL/科研/小论文1/基本试件打印参数正交试验/G%d/N%d.txt",Gi,ti);
        readdata = readtable(filesname, opts);
        % 转换为输出类型
        readdata = table2array(readdata);
        Flexural_x{Gi-9,ti}=readdata(:,4).*L0*100;   % 应变 6sh/L^2
        Flexural_y{Gi-9,ti}=readdata(:,1).*S0;        % 应力 3Fl/2bh^2
        clear readdata filesname
    end
end
clear opts
% 通过下方原始图像人工判断获得，亦可绘制后展示。
Flexural_start_line=[
    467	518	522	779	368	484
    519	337	270	390	415	252
    517	357	717	515	553	654
    770	553	453	344	498	318
    498	521	569	667	510	417
    360	431	407	364	472	358
    308	493	677	617	648	471
    379	529	641	781	915	439
    693	701	745	696	734	744];
Flexural_maxy_line=[
    4554	4543	4534	4729	4459	4457
    4391	4252	4080	4298	4264	4114
    4496	4306	4691	4533	4453	4458
    4431	4205	4134	3991	4242	4116
    4529	4633	4710	4794	4624	4509
    4078	4153	4095	4148	4172	4140
    4059	4267	4356	4334	4289	4146
    4165	4235	4366	4562	4611	4207
    4725	4725	4725	4725	4725	4725];

% % 绘制原始数据图像总图3x3
% figure()
% set(gcf,'unit','centimeters','Position',[5 5 30 30]);
% for Gi=10:18
%     subplot(3,3,Gi-9)
%     hold on
%     for ti=1:6
%         plot(Flexural_x{Gi-9,ti},Flexural_y{Gi-9,ti},LineWidth=1.5,Marker='o',MarkerEdgeColor='red',MarkerIndices=[Flexural_start_line(Gi-9,ti)])
%         %         plot(Flexural_y{Gi-9,ti},LineWidth=1.5)
%     end
%     xlabel("Strain (%)")
%     ylabel("Stress (MPa)")
%     % xlim([0 10])
%     % ylim([0,70])
%     % legend(Location="northwest",Box="off");
%     box('on');
%     title(['Flexural ',num2str(Gi)])
%     set(gca,'FontName','Arial','FontSize',FontSize)
% end
%
% % 分组绘制原始图像9张单独图（应力）
% for Gi=18:-1:10
%     figure()
%     set(gcf,'unit','centimeters','Position',[5 5 40 28]);
%     for ti=1:6
%         subplot(2,3,ti)
%         [~,maxy]=max(Flexural_y{Gi-9,ti});
%         Flexural_maxy_line(Gi-9,ti)=maxy;
%         plot(Flexural_y{Gi-9,ti},LineWidth=1.5,Marker='o',MarkerEdgeColor='red',MarkerIndices=[Flexural_start_line(Gi-9,ti),Flexural_maxy_line(Gi-9,ti)])
%         xlabel("N")
%         ylabel("Stress (MPa)")
%         % xlim([0 10])
%     % ylim([0,70])
% %     legend(Location="northwest",Box="off");
%     box('on');
%
%     title(['Flexural ',num2str(Gi-9), '-',num2str(ti)])
%     set(gca,'FontName','Arial','FontSize',FontSize)
%     end
% end

% 绘制受弯总图
Flexural_maxy=zeros(9,6);

% figure()
% set(gcf,'unit','centimeters','Position',[5 5 30 30]);
% for Gi=10:18
%     Gi=Gi-9;
%  subplot(3,3,Gi)
% 
%     hold on
%     for ti=1:6
%         Flexural_maxy(Gi,ti)=Flexural_y{Gi,ti}(Flexural_maxy_line(Gi,ti))-Flexural_x{Gi,ti}(Flexural_start_line(Gi,ti));
%         x=Flexural_x{Gi,ti}(Flexural_start_line(Gi,ti):end);x=x-x(1);
%         y=Flexural_y{Gi,ti}(Flexural_start_line(Gi,ti):end);y=y-y(1);
%         Flexural_maxy(Gi,ti)=max(y);
%         plot(x,y,LineWidth=1.5,Color=colors{Gi})
%     end
%     xlabel("Strain (%)")
%     ylabel("Stress (MPa)")
%     xlim([0 8])
%     ylim([0,70])
% 
%     legend(legends(Gi),Location="north",Box="off");
%     box('on');
%     title(['Flexural ',num2str(Gi)])
%     set(gca,'FontName','Arial','FontSize',FontSize)
% end
% 
% 
clear fh1
fh1=figure();
set(gcf,'unit','centimeters','Position',[15 15 15 15]);
axes('parent',fh1,'unit','centimeters','position',[1.5 1.5 12.5 12.5]);
hold on
LineStyles=repmat(["-",":","--"],1,3);
for Gi=10:18
    Gi=Gi-9;

    hold on
    for ti=1:6
        x=Flexural_x{Gi,ti}(Flexural_start_line(Gi,ti):end);x=x-x(1);
        y=Flexural_y{Gi,ti}(Flexural_start_line(Gi,ti):end);y=y-y(1);
        Flexural_maxy(Gi,ti)=max(y);
        plot(x,y,LineWidth=1.5,Color=colors{Gi},LineStyle=LineStyles(Gi))
    end

end

xlabel("Strain (%)")
    ylabel("Stress (MPa)")
    xlim([0 8])
    ylim([0,70])

    legend(legends9,Location="southeast",Box="off");
    box('on');
%     title(['Flexural'])
    set(gca,'FontName','Arial','FontSize',FontSize)

%% 汇总柱状图
% close all
com_bar=figure();
set(gcf,'unit','centimeters','Position',[5 5 28 28]);
axes('parent',com_bar,'unit','centimeters','position',[1.5 1.5 25 25],'fontsize',20);

x=1:9;
all_mean=reshape(mean([Tensile_maxy;compressive1_maxy;Flexural_maxy],2),9,3)';
all_mean_std=reshape(std([Tensile_maxy;compressive1_maxy;Flexural_maxy],0,2),9,3)';

bar_fig= bar(x,all_mean,FaceColor="flat",BarWidth=0.8);
hold on


for typei=1:3

    set(bar_fig(typei),'FaceColor',color3s{typei});

    errorbar(x+0.225*(typei-2),all_mean(typei,:),all_mean_std(typei,:),'LineStyle','none',LineWidth=1.5,Color=[0.5 0.5 0.5]);

    xtips = bar_fig(typei).XEndPoints;
    ytips = bar_fig(typei).YEndPoints+all_mean_std(typei,:);
    label = string(vpa(bar_fig(typei).YData,3));
    text(xtips,ytips,label,'HorizontalAlignment','center','VerticalAlignment','bottom','FontSize',12)
end
legend(["Tensile Yield Strength" "Compressive Yield Strength" "Flexural Yield Strength" ],Location="northwest",box="on")

%     xlabel("Type")
ylabel("Stress (MPa)")
set(gca,'FontName','Arial','FontSize',14)
%% 单独柱状图
figure();
ti=1;
set(gcf,'unit','centimeters','Position',[3 3 15 15]);
axes('unit','centimeters','position',[1.5 1.5 12.5 12.5]);
bar_fig_i= bar(x,all_mean(ti,:),BarWidth=0.8,FaceColor=color2);
hold on
errorbar(x,all_mean(ti,:),all_mean_std(ti,:),'LineStyle','none',LineWidth=1.5,Color=[0.5 0.5 0.5]);

xtips = bar_fig_i.XEndPoints;
ytips = bar_fig_i.YEndPoints+all_mean_std(ti,:);
label = string(vpa(bar_fig_i.YData,3));
text(xtips,ytips,label,'HorizontalAlignment','center','VerticalAlignment','bottom')


ylabel("Stress (MPa)")
set(gca,'FontName','Arial','FontSize',FontSize)
%% 保存计算获得的应力矩阵
all_sigma=[Tensile_maxy(:),compressive1_maxy(:),Flexural_maxy(:)];

save sigma_mat all_sigma all_mean all_mean_std  Tensile_maxy compressive1_maxy Flexural_maxy


%% 加载应力矩阵
load('sigma_mat.mat')
A=repmat([1	1	1	2	2	2	3	3	3]',6,1);
B=repmat([1	2	3	1	2	3	1	2	3]',6,1);
C=repmat([1	2	3	2	3	1	3	1	2]',6,1);
D=repmat([1	2	3	3	1	2	2	3	1]',6,1);
save taguchi9_mat A B C D % 保存相关变量
%% SN方差分析
sn=-10*log10(1./all_sigma.^2);
anova_soure_tbl= [{'A','B','C','D','Tensile','Compressive','Flexural'};[num2cell([A B C D]), num2cell(sn)]];
all_anova=cell(3,1);
for ni=1:3
    [~,tbli]=anovan(sn(:,ni),{A,B,C,D},"varnames",{'A','B','C','D'});
    tbli{1,1}='Factor';
    percent_source=cell2mat(tbli(2:6,5));
    percent=[{'Contribution'};num2cell(round(100*[percent_source;sum(percent_source)]./sum(percent_source),2))];


    all_anova{ni}=[tbli(:,[1 2 3 5 6 7]),percent];

end
save anovan_mat  all_anova anova_soure_tbl

%% 信噪比分析（采用越大越好计算公式）
all_sigma2square=1./all_sigma.^2;
all_sn=cell(3,1);
% all_ki=cell(3,1);
for ni=1:3
    [tbl_A,grp_A]=grpstats(all_sigma2square(:,ni),A,["mean","gname"]);
    [tbl_B,grp_B]=grpstats(all_sigma2square(:,ni),B,["mean","gname"]);
    [tbl_C,grp_C]=grpstats(all_sigma2square(:,ni),C,["mean","gname"]);
    [tbl_D,grp_D]=grpstats(all_sigma2square(:,ni),D,["mean","gname"]);

    sni=-10*log10([tbl_A,tbl_B,tbl_C,tbl_D]);
    sni_delta=max(sni)-min(sni);
    [~,sni_sort]=sort(sni_delta,'descend'); % 降序排列
    [~,sni_sort]=sort(sni_sort);           % 获得排位
    all_sn{ni}=cell(size(sni)+[3 1]);
    all_sn{ni}(1,:)={'Level' 'A','B','C','D'};
    all_sn{ni}(2:end,1)={'1' '2' '3' 'Delta' 'Rank'};
    all_sn{ni}(2:end,2:end)=[num2cell([sni;sni_delta]);num2cell(sni_sort)];

    % 直观分析计算过程（不采用）
    % tbl_A=grpstats(all_sigma(:,ni),A,@(x)sum(x,1));
    % tbl_B=grpstats(all_sigma(:,ni),B,@(x)sum(x,1));
    % tbl_C=grpstats(all_sigma(:,ni),C,@(x)sum(x,1));
    % tbl_D=grpstats(all_sigma(:,ni),D,@(x)sum(x,1));
    % all_ki{ni}=[tbl_A,tbl_B,tbl_C,tbl_D];

end
save all_sn_mat all_sn
%% 绘制信噪比SN曲线
close all
which_sigma=1; % 1:Tensile 2:Compressive1 3:Flexural
y1=cell2mat(all_sn{which_sigma}(2:4,2));
y2=cell2mat(all_sn{which_sigma}(2:4,3));
y3=cell2mat(all_sn{which_sigma}(2:4,4));
y4=cell2mat(all_sn{which_sigma}(2:4,5));
% y1=[0 0 0 ];
% y2=[0 0 0 ];
% y3=[0 0 0 ];
% y4=[0 0 0 ];
plot([y1;y2;y3;y4])
ylim_range=[18 32]; % set ylim range from the upper figure


FontSize=13;
fh1=figure();
set(gcf,'unit','centimeters','Position',[2 2 23 11]);

axes('parent',fh1,'unit','centimeters','position',[2 2 5 8],'fontsize',10);
plot([1 2 3],y1,'linewidth',3,'Marker','.','MarkerSize',30);
xlim([0.5 3.5])
ylim(ylim_range)
set(gca,'XTick',[1 2 3],'yTickLabel',num2str(get(gca,'yTick')','%.1f'));
xlabel("Filament","FontSize",FontSize,'FontName','Arial')
ylabel("Mean of SN Radio","FontSize",FontSize,'FontName','Arial')
box off

axes('parent',fh1,'unit','centimeters','position',[7 2 5 8],'fontsize',10);
plot([240 250 260],y2,'linewidth',3,'Marker','.','MarkerSize',30);
xlim([235 265])
ylim(ylim_range)
set(gca,'XTick',[240 250 260],'YTickLabel','')
xlabel("Speed(mm/s)","FontSize",FontSize,'FontName','Arial')
box off

axes('parent',fh1,'unit','centimeters','position',[12 2 5 8],'fontsize',10);
plot([210 215 220],y3,'linewidth',3,'Marker','.','MarkerSize',30)
xlim([205 225])
ylim(ylim_range)
set(gca,'XTick',[210 215 220],'YTickLabel','')
xlabel("Temperature(\circC)","FontSize",FontSize,'FontName','Arial')
box off

axes('parent',fh1,'unit','centimeters','position',[17 2 5 8],'fontsize',10);
plot([0.16 0.20 0.24],y4,'linewidth',3,'Marker','.','MarkerSize',30)
xlim([0.14 0.26])
ylim(ylim_range)
set(gca,'XTick',[0.16 0.20 0.24],'YTickLabel','')
xlabel("Layer thickness(mm)","FontSize",FontSize,'FontName','Arial')
box off

ax2 = axes('parent',fh1,'unit','centimeters','position',[2 2 20 8],...  
           'XAxisLocation','top',...  
           'YAxisLocation','right',...  
           'Color','none',...  
           'XColor','k','YColor','k');  
set(ax2,'YTick', []);  
set(ax2,'XTick', []); 
%% 线性回归预测
X=[ones(size(A)) A B C D];
all_regress=cell(3,1);
for ni=1:3
    [b,~,~,~,stats] = regress(all_sigma(:,ni),X);
    all_regress{ni}=[{'' 'A' 'B' 'C' 'D' 'R^2' 'p-value'}',num2cell([b;stats([1 3])'])];
end

save all_regress_mat all_regress

{% endfolding %}

home assistant安装方法

Thu, 10 Nov 2022 22:02:11 +0000

home assistant安装

现在实验室虽然在2楼还是很潮，打印用的材料总是受潮，所以老师让买了一个小米除湿器，这个设备可以直接加入米家远程控制，但我用的是iPhone想体验一下homekit，冲浪之后发现home assistant可以实现米家接入homekit，网上的教程有些代码过期了，因此整理一下安装过程。

docker安装

Mac下的docker

安装docker

部署镜像

#拉取hass的最新版镜像(注意：根据自己需求选择对应版本，并不是越新越好)
docker pull homeassistant/home-assistant:latest
#创建容器并运行
docker run -d --name="hass" -v /Users/wenqian/home-assistants:/config -p 8123:8123 -p homeassistant/home-assistant:latest
解释:
-d：表示在后台运行
–name：给容器设置一个名称，用于识别，这里起名叫做hass
-v：配置目录映射（容器内的指定目录直接映射到本地主机的对应目录，比如上面那句命令，意思就是将容器里hass的 config 目录映射到 /Users/wenqian/home-assistants 这个目录里面）
-p：映射端口（将容器内的8123端口直接映射到物理机的8123端口，这样服务器的IP地址+端口号即可访问hass页面）

访问hass

浏览器访问localhost:8123/ip:8123，设置用户名并登录，

安装hacs商店

官网链接 hacs

# 进入docker容器内
docker exec -it hass bash

# 下载
wget -O - https://get.hacs.xyz | bash -

hacs内安装xiaomi miot
hass内安装xiaomi miot
hass内安装apple homekit并设置桥接

docker内容器的ip不是桥接，因此homekit桥接失败，我还搞明白如何修改，看了网上代码好麻烦，直接放弃了

虚拟机安装

Mac虚拟机

我直接下载了官网推荐的VirtualBox虚拟机软件，但是无法安装64位虚拟机，遂放弃……

windows虚拟机hyper

从home assistant官网下载windows下hyper对应的虚拟机文件，参考官网设置，等待系统启动即可使用。

我们学校内有线网可以使用ipv6，访问github会快一点，所以虚拟交换机为与有线网卡的桥接，这样hass虚拟机直接获得了校内ipv4和校内ipv6，理论上我在校园网下手机homekit可以直接控制设备，不用额外搭配家庭中枢（HomePod、Apple TV）【我试了食堂的网络不行哈哈哈哈】

具体操作有人问了我再补充吧

搭建自己的随机图API

Tue, 04 Oct 2022 22:53:26 +0000

因为一直没有为博客的随机封面找到好的api，所以就在网上搜索如何自己制作，发现很简单，只需要：有图、写个php。

生成图片链接

首先选好自己的目标图片，上传至internet，可以选择自己服务器、云存储空间、GitHub图床等，获取所有图片外链即可，建议图片名字选择纯数字。

PHP生成随机链接

制作图片api的方式有很多，本文参考如下代码，通过php生成随机数字，拼接得到随机访问图片的链接，然后跳转对应链接。

<?php
//初始化随机数生成器种子，这行代码也可以删除
$seed = time();
//获取随机数
$num = rand(1,42);
//拼接图片地址
$picpath = "http://ritn3ta2k.hn-bkt.clouddn.com/api-img/".$num.".jpg";
//重定位到图片
die(header("Location: $picpath"));
?>

参考资料

Mac下Nginx的安装与配置

Sun, 02 Oct 2022 16:48:18 +0000

环境：MacBook Air M1, macOS Monterey 12.6

安装

建议选择homebrew安装，安装完成后会显示基本信息，后续如果要重新查看可通过命令实现。
brew install nginx 安装
brew info nginx 查看信息


Docroot is: /opt/homebrew/var/www

The default port has been set in /opt/homebrew/etc/nginx/nginx.conf to 8080 so that
nginx can run without sudo.

nginx will load all files in /opt/homebrew/etc/nginx/servers/.

To restart nginx after an upgrade:
  brew services restart nginx
Or, if you don't want/need a background service you can just run:
  /opt/homebrew/opt/nginx/bin/nginx -g daemon off;

使用

brew services start nginx 启动 brew services stop nginx 停止 nginx -t 测试配置文件

Mac个人使用技巧随时更

Sat, 24 Sep 2022 14:23:45 +0000

2021年12月31日拥有了我的第一台苹果电脑MacBook Air Apple M1 芯片，之前没有用过Mac系统因此会遇到一些操作上的差异，当然这个过程中也有一些好用的功能，这个文章就记录一些我认为有趣的地方。

快速定位目标文件/文件夹

在访达界面中，使用 Shift-Command-G ，打开“前往文件夹”窗口，输入路径即可 - 在其他软件中选择本地文件时，使用**/**，打开“前往文件夹”窗口，输入路径即可

搬运官网 Mac 键盘快捷键

Sat, 24 Sep 2022 14:23:41 +0000

{callout color="#f0ad4e"} 本文来源于Apple官网，原文发布日期： 2022 年 07 月 27 日。 {/callout}

您可以按下某些组合键来实现通常需要鼠标、触控板或其他输入设备才能完成的操作。

要使用键盘快捷键，请按住一个或多个修饰键，然后按快捷键的最后一个键。例如，要使用 Command-C（拷贝），请按住 Command 键并按 C 键，然后同时松开这两个键。Mac 菜单和键盘通常使用符号来表示某些按键，其中包括以下修饰键：

Command（或 Cmd）⌘
Shift ⇧
Option（或 Alt）⌥
Control（或 Ctrl）⌃
Caps Lock ⇪
Fn

{alert type=“info”} 在 Windows PC 专用键盘上，请用 Alt 键代替 Option 键，用 Windows 标志键代替 Command 键。 Apple 键盘上的某些按键具有特殊符号和功能，例如显示屏亮度、键盘亮度、调度中心等。如果您的键盘上没有这些功能，您也许可以通过创建自己的键盘快捷键来实现其中的一些功能。要将这些键用作 F1、F2、F3 或其他标准功能键，请将它们与 Fn 键组合使用。 {/alert}

剪切、拷贝、粘贴和其他常用快捷键

Command-X ：剪切所选项并拷贝到剪贴板。 Command-C ：将所选项拷贝到剪贴板。这同样适用于“访达”中的文件。 Command-V ：将剪贴板的内容粘贴到当前文稿或 App 中。这同样适用于“访达”中的文件。 Command-Z ：撤销上一个命令。随后您可以按 Shift-Command-Z 来重做，从而反向执行撤销命令。在某些 App 中，您可以撤销和重做多个命令。 Command-A ：全选各项。 Command-F ：查找文稿中的项目或打开“查找”窗口。 Command-G：再次查找：查找之前所找到项目出现的下一个位置。要查找出现的上一个位置，请按 Shift-Command-G。 Command-H：隐藏最前面的 App 的窗口。要查看最前面的 App 但隐藏所有其他 App，请按 Option-Command-H。 Command-M：将最前面的窗口最小化至“程序坞”。要最小化最前面的 App 的所有窗口，请按 Option-Command-M。 Command-O：打开所选项，或打开一个对话框以选择要打开的文件。 Command-P：打印当前文稿。 Command-S：存储当前文稿。 Command-T：打开新标签页。 Command-W：关闭最前面的窗口。要关闭 App 的所有窗口，请按下 Option-Command-W。 Option-Command-Esc：强制退出 App。 Command-空格键：显示或隐藏“聚焦”搜索栏。要从“访达”窗口执行“聚焦”搜索，请按 Command–Option–空格键。（如果您使用多个输入源以便用不同的语言键入内容，这些快捷键会更改输入源而非显示“聚焦”。了解如何更改冲突的键盘快捷键。） Control-Command–空格键：显示字符检视器，您可以从中选择表情符号和其他符号。 Control-Command-F：全屏使用 App（如果 App 支持）。空格键：使用“快速查看”来预览所选项。 Command-Tab：在打开的 App 中切换到下一个最近使用的 App。 Shift-Command-5：在 macOS Mojave 或更高版本中，拍摄截屏或录制屏幕。也可以使用 Shift-Command-3 或 Shift-Command-4 来拍摄截屏。进一步了解截屏。 Shift-Command-N：在“访达”中创建一个新文件夹。 Command-逗号 (,)：打开最前面的 App 的偏好设置。

Matlab 命令记录

Fri, 27 May 2022 14:46:54 +0000

cell

定义空cell

定义一个2✖️3的元胞结构数组

>> x=cell(2,3)

x =

  2×3 cell 数组

    {0×0 double}    {0×0 double}    {0×0 double}
    {0×0 double}    {0×0 double}    {0×0 double}

cell读取时{}()的区别

()读取得到的是元包结构，即cell数组
{}读取得到的是数据本身的内容，例如字符串

% 这个例子
>> title={{'pitch-RAO'} {'T'}    {'heave-波能俘获效率'} }

title =

  1×3 cell 数组

    {1×1 cell}    {1×1 cell}    {1×1 cell}

>> title(1)

ans =

  1×1 cell 数组

    {1×1 cell}

>> title{1}

ans =

  1×1 cell 数组

    {'pitch-RAO'}

>> title(1)(1)
错误:  无效的数组索引。
 
>> title(1){1}
错误:  无效的数组索引。
 
>> title{1}{1}

ans =

    'pitch-RAO'

>> title{1}(1)

ans =

  1×1 cell 数组

    {'pitch-RAO'}

坐标轴科学记数法表示与轴范围的影响

xlim([xmin xmax]) x轴显示范围
ylim([ymin ymax]) y轴显示范围
xticklabel xtick x坐标轴标签
yticklabel ytick x坐标轴标签

Excel-vba创建表及粘贴公式结果

Thu, 26 May 2022 22:03:46 +0000

VBA创建新表汇总子sheet结果

源代码

V1.0


Public Sub abs_and_average()
    On Error Resume Next  '这样不会报错。下标越界
    If IsEmpty(Worksheets("data")) Then
            Sheets.Add after:=Sheets(Sheets.Count)
          Sheets(Sheets.Count).Name = "data"
      Else
           Worksheets("data").UsedRange.ClearContents
      End If
  
  
  insertline=0
  For Each usedsheet In Worksheets:
  insertline=insertline+1
      If usedsheet.Range("b12") = "Cd3" Then
          usedsheet.Range("w14").FormulaR1C1 = "齿轮箱输入转速(r/min)的绝对值"
          usedsheet.Range("p10").FormulaR1C1 = "输入轴转速"
      
          '单元格输入公式，相对引用（可以通过录制宏进行修改替换）
          usedsheet.Range("w15").FormulaR1C1 = "=abs(RC[-1])"
          
      
          '判断需要填充多长，这里从第30行开始向后循环判断非空单元格
          i = 30
          While usedsheet.Cells(i + 1, 10) <> ""
              i = i + 1
          Wend
          cc = "w15:v" & i
      
          '使用range的方法，向下填充，还有其他功能，可以通过录制宏查看
          usedsheet.Range(cc).FillDown
          usedsheet.Range("q10").Formula = "=AVERAGE(U3964:U5014)"


        '把平均值输出到汇总表
        Worksheets("data").cells(insertline,1)=split(usedsheet.name,"=")(0)
        Worksheets("data").cells(insertline,2)=split(usedsheet.name,"=")(1)
        Worksheets("data").cells(insertline,3)=usedsheet.Range("Q10").value

      Else
  
      End If
  Next
  ' 显示适当列宽
 Worksheets("data").range("A:c").EntireColumn.AutoFit
 Worksheets("data").Sort.SortFields.Clear
 Worksheets("data").Sort.SortFields.Add2 Key:=Range("B1:B20"), _
      SortOn:=xlSortOnValues, Order:=xlAscending, DataOption:=xlSortNormal
  With Worksheets("data").Sort
      .SetRange Range("A1:C20")
      .Header = xlGuess
      .MatchCase = False
      .Orientation = xlTopToBottom
      .SortMethod = xlPinYin
      .Apply
  End With
  End Sub

V2.0

Public Sub abs_and_average0528()
    On Error Resume Next  '这样不会报错。下标越界
    If IsEmpty(Worksheets("data")) Then
            Sheets.Add after:=Sheets(Sheets.Count)
          Sheets(Sheets.Count).Name = "data"
      Else
           Worksheets("data").UsedRange.ClearContents
      End If
  
  
  insertline=0
  For Each usedsheet In Worksheets:
  insertline=insertline+1
      If usedsheet.Range("b10") = "Cd3" Then
        '判断需要目标行数，这里从第30行开始向后循环判断非空单元格
        i = 30
        While Cells(i+1, 10) <> ""
            i = i + 1
        Wend
        cc = "u12:x" & i
        usedsheet.range(cc).clear
        '单元格输入文字
        usedsheet.Range("U12").FormulaR1C1 = "发电机转速（r/min）"
        usedsheet.Range("v12").FormulaR1C1 = "齿轮箱输入转速(r/min)"
        usedsheet.Range("w12").FormulaR1C1 = "【绝对值】发电机转速（r/min）"
        usedsheet.Range("x12").FormulaR1C1 = "【绝对值】齿轮箱输入转速(r/min)"


        '单元格输入公式，相对引用
        usedsheet.Range("U13").FormulaR1C1 = "=RC[-1]*60/2/3.1415926"
        usedsheet.Range("V13").FormulaR1C1 = "=RC[-19]*60/2/3.1415926"
        '单元格输入公式，直接输入内容
        usedsheet.Range("w13").Formula = "=abs(u13)"
        usedsheet.Range("x13").Formula = "=abs(v13)"

          
      
          '判断需要填充多长，这里从第30行开始向后循环判断非空单元格
          i = 30
          While usedsheet.Cells(i + 1, 10) <> ""
              i = i + 1
          Wend
          cc = "u13:x" & i
      
          '使用range的方法，向下填充，还有其他功能，可以通过录制宏查看
          usedsheet.Range(cc).FillDown


          usedsheet.Range("w11").Formula = "=AVERAGE(w1412:w50012)"
          usedsheet.Range("x11").Formula = "=AVERAGE(x1412:x50012)"


        '把平均值输出到汇总表
        Worksheets("data").cells(1,1)= "sheet name"
        Worksheets("data").cells(1,2)="【绝对值】发电机转速（r/min）"
        Worksheets("data").cells(1,3)="【绝对值】齿轮箱输入转速(r/min)"
        Worksheets("data").cells(insertline+1,1)=usedsheet.name
        Worksheets("data").cells(insertline+1,2)=usedsheet.Range("w11").value
        Worksheets("data").cells(insertline+1,3)=usedsheet.Range("x11").value

      Else
  
      End If
  Next
' 显示适当列宽
  Worksheets("data").range("A:c").EntireColumn.AutoFit
'  Worksheets("data").Sort.SortFields.Clear
'  Worksheets("data").Sort.SortFields.Add2 Key:=Range("B1:B20"), _
'       SortOn:=xlSortOnValues, Order:=xlAscending, DataOption:=xlSortNormal
'   With Worksheets("data").Sort
'       .SetRange Range("A1:C20")
'       .Header = xlGuess
'       .MatchCase = False
'       .Orientation = xlTopToBottom
'       .SortMethod = xlPinYin
'       .Apply
'   End With
  End Sub

python-CAD绘制路基横断面

Tue, 24 May 2022 22:37:22 +0000

基于python的CAD路基横断面绘制

源代码

# _*_ coding:utf-8- _*_

from pyautocad import APoint, Autocad, aDouble
import win32com.client
import numpy as np
"""
设计高程：(2)1654.0m
填土边坡与路堑边坡坡率：(1)填方上部1:1.2(高度≤8m)，填方下部1:1.4(高度≤12m) ；挖方1:0.9(高度≤20m)
路基宽度：(2) 路面宽7.0米，硬路肩宽0.75米，土路肩宽0.75米
以下数据单位：m
"""
design_height = 1654.0  # 设计高程
slop_fill_up = 1.2      # 坡率填方上部(高度≤8m)
slop_fill_down = 1.4    # 坡率填方下部(高度≤12m)
slop_cut = 0.9          # 坡率挖方(高度≤20m)
surface_width = 7.0     # 路面宽
hard_shoulder = 0.75    # 硬路肩宽
dirt_shoulder = 0.75    # 土路肩宽

ctfs = 10               # 出图坐标放缩倍数
distance_sections = 60  # 修改为出图时横断面桩号间距
# 连接excel读取数据

xlapp = win32com.client.Dispatch("Excel.Application")
wb = xlapp.ActiveWorkbook
ws = xlapp.ActiveWorkbook.ActiveSheet


# points是原始地面点高程
ground_points = []
row = 2
while row <= 20:
    i = 6               # 左侧数据
    j = 8               # 右侧数据

    # 桩号点绘图坐标
    x0 = 0 + (distance_sections * (row - 2))//3
    y0 = ws.cells(row, 7).value
    p0 = APoint(x0, y0)
    pts = []
    pts.append(p0)


    #  拾取左侧点数据
    while i >= 1:
        try:
            y_ = ws.cells(row, i).value
            x_ = - ws.cells(row + 1, i).value
            pti = pts[0] + (x_, y_, 0)
            pts.insert(0, pti)
        except:
            pass
        i -= 1
    #  拾取右侧点数据
    while ws.cells(row, j).value:
        y_ = ws.cells(row, j).value
        x_ = ws.cells(row + 1, j).value
        ptj = pts[-1] + (x_, y_, 0)
        pts.append(ptj)
        j += 1
    ground_points.append(pts)

    for i in pts:
        print(i)
    print("")
    row = row + 3



# 连接cad，并输出结果

acad = Autocad(create_if_not_exists=True)
acad.prompt("connect with the Autocad successfully.")
print("connect with the Autocad successfully.")

dwg = acad.ActiveDocument
"原始地面线"
# 定义原始地面线 图层
ysdmx = acad.doc.layers.add("原始地面线")
ysdmx.color = 7
ysdmx.LineWeight = 25
dwg.ActiveLayer = ysdmx
# 绘制原始地面线，使用多段线（PLine）
for line_pts in ground_points:
    pt = aDouble([j for i in line_pts  for j in i * ctfs])
    acad.model.AddPolyLine(pt).LineWeight = 25

"路基线"
# 绘制路基线，所有横断面路基线按一样绘制，手动修改
# 生成路基线关键点坐标
orginal_middl_pts = []
pt0 = APoint(0, design_height)
orginal_middl_pts.append(pt0)
rightpts = [(surface_width/2, 0), (hard_shoulder, 0), (dirt_shoulder, 0), (0.6, - 0.6), (0.6, 0), (0.6, 0.6), (1, 0), (20 * slop_cut, 20)]
leftpts = [(surface_width/2, 0), (hard_shoulder, 0), (dirt_shoulder, 0), (8 * slop_fill_up, 8), (2, 0), (12 * slop_fill_down, 12)]


for i in leftpts:
    pti = orginal_middl_pts[0] - (i[0], i[1], 0)
    orginal_middl_pts.insert(0, pti)
for j in rightpts:
    ptj = orginal_middl_pts[-1] + (j[0], j[1], 0)
    orginal_middl_pts.append(ptj)
# 定义路基线 图层
ljx = acad.doc.layers.add("路基线")
ljx.color = 1
ljx.LineWeight = 100
dwg.ActiveLayer = ljx
# 绘制路基，使用多段线（PLine）
pt = aDouble([i for j in orginal_middl_pts  for i in j * ctfs])
ljx_drawn = acad.model.AddPolyLine(pt)
ljx_drawn.LineWeight = 100




"道路中心线"
# 得到坐标

# 定义图层属性
zxj = acad.doc.layers.add("道路中心线")
zxj.color = 2
zxj.LineWeight = 25
try:
    acad.ActiveDocument.Linetypes.Load("ACAD_ISO10W100", "acadiso.lin")
except:
    pass
zxj.Linetype = "ACAD_ISO10W100"
dwg.ActiveLayer = zxj
# 绘制中线
zxj_drewn = acad.model.addline(ctfs * (pt0 - (0, 2, 0)), ctfs * (pt0 + (0, 2, 0)))
zxj_drewn.LineWeight = 25
# [中心线、路基线]向右复制6次
for i in range(6):
    _ = ljx_drawn.copy()
    _.Move(APoint(0, 0), APoint(distance_sections * (i + 1), 0) * ctfs)
    del _
    _ = zxj_drewn.copy()
    _.Move(APoint(0, 0), APoint(distance_sections * (i + 1), 0) * ctfs)

"添加桩号标注"
# 得到坐标
# 定义图层
zh = acad.doc.layers.add("桩号标注")
dwg.ActiveLayer = zh
# 添加注释
for i in range(7):
    text_string = "K58+%d" % (70 + 10 * i)
    insertpt = APoint((distance_sections * i), design_height - 7) * ctfs
    height = 2.5 * ctfs
    textobj = acad.model.addtext(text_string, insertpt, height)
    textobj.Alignment = 7
    textobj.textalignmentpoint = insertpt

zerotier虚拟局域网/内网穿透

Tue, 24 May 2022 22:21:55 +0000

zerotier创建虚拟局域网

zerotier是一个可以创建基于p2p的虚拟局域网，第一次使用只需安装软件登录网页管理虚拟ip，以后每次启动都可以和局域网内的设备通讯，路由器也可以安装。

可以实现的操作

利用ip远程电脑，高画质，网速稳定（利用微软的默认程序）
虚拟局域网文件共享（类似于上述远程）
打印机

设置流程

注册账号

打开zerotier注册账号，可以用其他账号登录。

下载对应客户端

各种操作系统都有对应软件，下载安装即可。

创建子网络

创建
加入
二次确认
查看/设置ip

开始使用

后续如果需要设置网络，可以通过官网页面登录，也可以通过my.zerotier.com登录设置，少了一个页面。得到了局域网可以干什么，请参考上面写的用法。

python-CAD绘制等高线

Tue, 24 May 2022 22:07:53 +0000

基于python的CAD等高线绘制

这是我本科在做《路基工程课程设计》挡土墙设计的时候，老师只给了原始地面数据点，因为疫情在家时间充裕，我就摸索了python画CAD，最终实现了绘制不同等高距的等高线绘制。 python面向对象编程。

基础知识

我这个代码只适合地面趋势简单的，通过差值得到固定整数高程点，然后将相同高程点依次连接。

note：

当时是在windows先写的代码，因为要用到win32com.client这个库，所以我估计 Mac下应该用不了，但可以修改代码绘制dxf文件，不调用CAD。
后续我再补充 Excel文件、输出CAD的文件及效果图

源代码

# _*_ coding:utf-8- _*_

from pyautocad import APoint, Autocad, aDouble
import win32com.client
import numpy as np

dgj = 0.5       # 等高距
ctfs = 10       # 放大倍数，本身采取m为单位绘制，可能放大10倍后输出效果好些

"""
等高线计算
"""
#连接excel读取数据，这个地方可以做一个前处理，直接把excel数据作为数组或者numpy数组进行处理

xlapp = win32com.client.Dispatch("Excel.Application")
wb = xlapp.ActiveWorkbook
ws = xlapp.ActiveWorkbook.ActiveSheet
"""
K58+70向右走为前进方向，绘制在
"""

"""
First
定义一个类，输入 list，其元素为APoint
内部函数得到插值的数据点，打印点数据
"""
class nei_cha:
    def __init__(self, pnts = [APoint(0,1)], dgj = 0.5, gu_ding_cha_zhi = "X" ):
        # gu_ding_cz = "X" 意思是不进行固定插值操作
        self.pnts = pnts
        self.dgj = dgj
        self.gu_ding_cha_zhi = gu_ding_cha_zhi

    # 输入两个三维坐标，x/y是相等的，内插得到所需高程
    # 返回的是，该区间内得到的点组成的list
    def do_two_points(self, p1 = APoint(0), p2 = APoint(0)):
        if self.gu_ding_cha_zhi == "X":
            nmax = int(max(p1.z, p2.z) // self.dgj)
            nmin = int(min(p1.z, p2.z) // self.dgj)
            if nmax == nmin:        # 这两个点之间内插不出来
                return 0
        else:
            nmax = self.gu_ding_cha_zhi
            nmin = self.gu_ding_cha_zhi - 1
            if (p1.z - nmax) * (p2.z - nmin) > 0:        # 固定插值不在区间内
                return 0

        if p1.x == p2.x:        # x坐标一样，需要内插Z到y
            """
            求直线方程，y
            """
            k = (p1.y - p2.y)/(p1.z - p2.z)
            b = p1.y - k * p1.z

            if self.gu_ding_cha_zhi == "X":
                ni = np.array(range(nmin, nmax)).reshape(int(nmax - nmin), 1)
                zi = (ni + 1) * self.dgj
            else:
                ni = np.array([self.gu_ding_cha_zhi]).reshape(int(nmax - nmin), 1)
                zi = ni
            yi = np.around(k * zi + b, 3)
            #yi = np.interp(zi, [p1.z, p2.z], [p1.y, p2.y])
            xi = np.full((len(yi), 1), p1.x)
            pts_np = np.hstack((xi, yi, zi))
            return  pts_np

        else:                   # y坐标一样，需要内插Z得到x
            return 0
    def do_it(self):
        pnts_np = np.array([])
        for i in range(len(self.pnts) - 1):

            p_np_i = self.do_two_points(self.pnts[i], self.pnts[i+1])
            if type(p_np_i) != type(0):
                try:
                    pnts_np = np.vstack((pnts_np, p_np_i))
                except:
                    pnts_np = p_np_i
        if len(pnts_np) != len(np.array([])):
            return pnts_np
        else:
            return 0



# points是原始地面点高程
ground_points = []
row = 2

while row <= 20:
    i = 6               # 左侧数据
    j = 8               # 右侧数据

    # 桩号点绘图坐标
    x0 = ((row - 2) / 3) * 10 + 70
    y0 = 0
    z0 = ws.cells(row, 7).value
    p0 = APoint(x0, y0, z0)
    pts = []
    pts.append(p0)


    #  拾取左侧点数据
    while i >= 1:
        try:
            z_ = ws.cells(row, i).value
            y_ = ws.cells(row + 1, i).value
            pti = pts[-1] + (0, y_, z_)
            pts.append(pti)
        except:
            pass
        i -= 1
    #  拾取右侧点数据
    while ws.cells(row, j).value:
        try:
            z_ = ws.cells(row, j).value
            y_ = - ws.cells(row + 1, j).value
            ptj = pts[0] + (0, y_, z_)
            pts.insert(0, ptj)
        except:
            print("右侧点拾取失败了：cells(%d,%d)"% (row, j))
            pass
        j += 1
    ground_points.append(pts)

    for i in pts:
        print(i)
    print("")
    row = row + 3

# 利用上方点开始使用类进行，插值
all_np = []
for i in ground_points:
    cha_zhi_la = nei_cha(i, dgj)
    a = cha_zhi_la.do_it()
    try:
        all_np = np.vstack((all_np, a))
    except:
        all_np = a
pts_all = [APoint(list(i)) for i in all_np]
pts_all.sort(key = lambda x:x.z)
del all_np

h_min = pts_all[0].z
h_max = pts_all[-1].z


zzzi = h_min
sorted_pts = []
while zzzi <= h_max:
    _ = [i for i in pts_all if i.z == zzzi]
    _.sort(key = lambda x:x[0])
    sorted_pts.append(_)
    zzzi += dgj
del pts_all
for i in sorted_pts:
    elevation = i[0].z
    print("高程为 %s 的数据点"% elevation)
    for l_ in i:
        print(l_)
# 连接cad，并输出结果

del sorted_pts[0:11]
del sorted_pts[-10:-1]
acad = Autocad(create_if_not_exists=True)
acad.prompt("connect with the Autocad successfully.")
print("connect with the Autocad successfully.")

dwg = acad.ActiveDocument

# 定义等高线图层
dgx = acad.doc.layers.add("等高线")
dgx.LineWeight = 25
dwg.ActiveLayer = dgx
dgxbz = acad.doc.layers.add("等高线标注")


for m in range(len(sorted_pts)):
    i = sorted_pts[m]
    elevation = i[0].z
    print("高程为 %s 的数据点"% elevation)
    for l_ in i:
        print(l_)
    try:


        for j in i:
            j.z = 0
        pt = aDouble([j for _ in i for j in _ * ctfs])
        hl = acad.model.AddPolyLine(pt)
        hl.LineWeight = 25
        hl.Elevation = elevation
        hl.Thickness = 900
        if m % 5 == 0:
            hl.LineWeight = 100
            hl.color = 1
            text_string = elevation
            insertpt = i[0] * ctfs
            textobj = acad.model.addMtext(insertpt, 26, text_string)
            textobj.Height = 8
            textobj.BackgroundFill = True
            textobj.color = 1
            textobj.layer = "等高线标注"

    except:
        print('*' * 10, "高程为 %s 的数据点不足，跳过绘制此等高线"% elevation)
        pass

you-get命令行视频下载工具

Mon, 23 May 2022 23:57:36 +0000

you-get介绍

一、安装

更多方式请参考you-get GitHub主页

环境依赖

以下为必要配置:

Python 3.5 或更高版本
FFmpeg 1.0 或更高版本
(可选) RTMPDump

方法 1: 通过 pip 安装（Mac/Windows）

首先配置python3环境，然后通过pip命令直接安装you-get，如果电脑只有python3一个版本，则使用pip install you-get即可，如果有python2，则使用pip3命令，即pip3 install you-get

The official release of you-get is distributed on PyPI, and can be installed easily from a PyPI mirror via the pip package manager. Note that you must use the Python 3 version of pip:

pip3 install you-get

方法 2: Homebrew (Mac only)

You can install you-get easily via:

Excel vba修改数据图表格式

Mon, 23 May 2022 22:04:42 +0000

场景描述

修改坐标轴区间及交叉位置

代码

Sub change_pic()
'
' 修改数据图像格式，不需要的可以直接注释或修改参数
'    
    'y轴的显示区间，设置为自动，如果输入自定义，使用.MinimumScale进行赋值
    ActiveChart.Axes(xlValue).MinimumScaleIsAuto = True
    ActiveChart.Axes(xlValue).MaximumScaleIsAuto = True
    'x轴与y轴对齐的y坐标，这里读取了上一步的y轴最小值，即坐标轴放在下方
    ActiveChart.Axes(xlValue).CrossesAt = ActiveChart.Axes(xlValue).MinimumScale
End Sub

ffmpeg压缩图片尺寸

Mon, 23 May 2022 21:26:01 +0000

FFmpeg使用说明

简介

FFmpeg是一个跨平台可录制、转换、流式传输音视频工具。（自已从官网翻译的）

用户文档

图片处理

查看图像信息

ffprobe image_source 
#查看图像尺寸等信息

图片压缩

ffmpeg -i image_source -q compress_level out_source
#compress_level：压缩等级为正整数，数值越大压缩等级越高

约束比例压缩

ffmpeg -i image_source -vf scale=width:height out_source
#width为压缩后的图片宽度
#height为压缩后的图片高度
#当width或height其中有一个数值为-1时将保持原来图片尺寸比例压缩

图片裁剪

ffmpeg -i image_source -vf crop=width:height:from_x:from_y out_source
#width是要裁剪出的图片的宽度
#height是要裁剪出的图片的高度
#from_x是裁剪的起始X轴的位置
#from_y是裁剪的起始Y轴的位置
#如果不指定from_x和from_y表示从原始图片中心位置开始裁剪

示例

这个场景来源于我在课题组做推送，我进行视频分享时会录制gif动图，因为文案每次都是轮流写的，不同人有不同的录制方式，导致有时录屏转的文件格式很大（可能有15MB），但微信推送图片最大是10MB，所以必须进行压缩。

最简单的办法就是，直接转换录制的视频或者已经压缩后的超大gif文件，附件给了2个mov文件。

ffmpeg -i 1.mov 1-1.gif
ffmpeg -i 1.mov  -vf scale=500:-1 1-2.gif

附件下载

ffmpeg-test202207.zip

参考资料

FFmpeg—-图片处理

结构可靠度理论代码共享

Sat, 14 May 2022 22:45:14 +0000

本文将记录《结构可靠度理论》的一些代码及课堂笔记，采用Matlab编程

1 概率分布密度（pdf、cdf）

自己编写函数

clc
clear
close all
% unif
a=1;
b=3;
x=1:0.1:3;
for i=1:length(x)
    fx(i)=1/(b-a);
    Fx(i)=(x(i)-a)/(b-a);
end
h1=figure(1);
ax1=axes('Parent',h1,'Position',[0.1,0.4,0.35,0.35],'FontSize',10);
plot(x,fx,'r-')
xlabel('x','fontsize',10)
ylabel('pdf','FontSize',10)

ax2=axes('Parent',h1,'Position',[0.55,0.4,0.35,0.35],'FontSize',10);
plot(x,Fx,'r-')
xlabel('x','fontsize',10)
ylabel('cdf','FontSize',10)

%-------
mux=1;sigx=1;
x=-2:0.1:4;
syms t
for i=1:length(x)
    tx=x(i);
    a1=1/(sigx*(2*pi)^0.5);
    b1=exp(-0.5*((tx-mux)/sigx)^2);
    fx(i)=a1*b1;
    c1=int(exp(-0.5*((t-mux)/sigx)^2),t,-inf,tx);
    Fx(i)=a1*c1;
end

h1=figure(2);
ax1=axes('Parent',h1,'Position',[0.1,0.4,0.35,0.35],'FontSize',10);
plot(x,fx,'r-')
xlabel('x','fontsize',10)
ylabel('pdf','FontSize',10)

ax2=axes("Parent",h1,"Position",[0.55,0.4,0.35,0.35],"FontSize",10);
plot(x,Fx,'r-')
xlabel('x','fontsize',10)
ylabel('cdf','FontSize',10)

调用Matlab自带函数

clc
clear
close all
%-----
a=1;b=3;
x=-2:0.1:5;
fx=unifpdf(x,a,b);
Fx=unifcdf(x,a,b);
figure(1)
plot(x,fx,'-r')
hold on
plot(x,Fx,'--b')

%---------
mux=0;sigx=1;
x=-3:0.1:3;
fx=normpdf(x,mux,sigx);
Fx=normcdf(x,mux,sigx);

figure(2)
plot(x,fx,'r-')
hold on
plot(x,Fx,'r-')

Posts on 极客Vinkey

Cisco AnyConnect优雅脚本与openconnect进阶教程

概述

仓库结构

移动端

Cisco（桌面端，仅命令行）

Windows快速切换简体繁体拼音输入-Ctrl+Shift+F

VSCode+Latex Workshop环境配置

安装Tex

Windows

Ubuntu

我的科研绘图调色板

🌈 我的科研绘图调色板

双色

三色

端口转发

Windows

创建端口转发

查看端口使用

Ubuntu

如何将Cisco Anyconnect共享给局域网使用

所需软件

WSL配置

WSL中安装3Proxy并开启Socks代理

利用FRP将端口穿透至局域网

frps 服务端

frpc 客户端

记录一下脾胃虚弱中药方

症状表现

生活调理

中医调理

药方

痊愈

Mac查看端口占用和释放端口的教程

清除iPhone短信app右上角的小红点

来源

曲折

正解

问卷调查数据采用Excel处理生成数据看板

月总结-2023

3月很ok

4月是养精蓄锐

Matlab绘图代码共享-硕士阶段

拉压弯曲线图及柱状图2022/12

home assistant安装方法

docker安装

Mac下的docker

安装docker

部署镜像

访问hass

安装hacs商店

虚拟机安装

Mac虚拟机

windows虚拟机hyper

搭建自己的随机图API

生成图片链接

PHP生成随机链接

Mac下Nginx的安装与配置

安装

使用

Mac个人使用技巧 随时更

快速定位目标文件/文件夹

搬运官网 Mac 键盘快捷键

剪切、拷贝、粘贴和其他常用快捷键

Matlab 命令记录

cell

定义空cell

cell读取时{}()的区别

坐标轴科学记数法表示与轴范围的影响

Excel-vba创建表及粘贴公式结果

VBA创建新表汇总子sheet结果

源代码

V1.0

V2.0

python-CAD绘制路基横断面

基于python的CAD路基横断面绘制

源代码

zerotier虚拟局域网/内网穿透

zerotier创建虚拟局域网

可以实现的操作

Mac个人使用技巧随时更