for j=1:num_layers lamdpor(j) = lamds(j) .* (1 - phi(j)) + phi(j) .* lamdaf(j); dtt=dz*qg/lamdpor(j); %竖向网格温差 if(j==1) j00=2; j01=nz/4+1; elseif(j==2) j00=nz/4+2; j01=2*nz/4+1; elseif(j==3) j00=2*nz/4+2; j01=3*nz/4+1; else j00=3*nz/4+2; j01=nz+1; end for i=j00:j01 Tc=Tc+dtt; if(i<=nj) T0f1(i)=Tc; T0f2(i)=Tc; end end end优化代码

时间: 2024-01-15 07:02:36 浏览: 27
可以将外层循环和内层循环合并,使用向量化方式进行优化,并使用数组索引代替if-elseif-else语句。具体代码如下: ``` j0 = [2, nz/4+2, 2*nz/4+2, 3*nz/4+2, nz+1]; % 计算j00和j01的起始值 j00 = j0(1:end-1); j01 = j0(2:end) - 1; lamdpor = lamds .* (1 - phi) + phi .* lamdaf; % 计算lamdpor dzqg = dz * qg; % 计算dz*qg Tc = cumsum(dzqg ./ lamdpor) + Ta; % 计算温度,加上Ta for j = 1:num_layers T0f1(j00(j):min(j01(j),nj)) = Tc(j00(j):min(j01(j),nj)); % 将温度赋值给T0f1 T0f2(j00(j):min(j01(j),nj)) = Tc(j00(j):min(j01(j),nj)); % 将温度赋值给T0f2 end ``` 这样就可以将多重循环优化为一次矩阵运算,大大提高了代码执行效率,同时使用数组索引代替if-elseif-else语句,代码更加简洁易读。
相关问题

Tc=qg/ha+Ta; %地表温度 T0f1(1)=Tc; T0f2(1)=Tc; for j=1:num_layers lamdpor(j) = lamds(j) .* (1 - phi(j)) + phi(j) .* lamdaf(j); dtt=dz*qg/lamdpor(j); %竖向网格温差 if(j==1) j00=2; j01=nz/4+1; elseif(j==2) j00=nz/4+2; j01=2*nz/4+1; elseif(j==3) j00=2*nz/4+2; j01=3*nz/4+1; else j00=3*nz/4+2; j01=nz+1; end for i=j00:j01 Tc=Tc+dtt; if(i<=nj) T0f1(i)=Tc; T0f2(i)=Tc; end end end优化代码

可以将外层循环和内层循环合并,使用向量化方式进行优化。具体代码如下: ``` lamdpor = lamds .* (1 - phi) + phi .* lamdaf; % 计算lamdpor dzqg = dz * qg; % 计算dz*qg Tc = cumsum(dzqg ./ lamdpor) + Ta; % 计算温度,加上Ta T0f1(1:nj) = Tc(j00:j01); % 将温度赋值给T0f1 T0f2(1:nj) = Tc(j00:j01); % 将温度赋值给T0f2 ``` 这样就可以将多重循环优化为一次矩阵运算,大大提高了代码执行效率。注意要将地表温度Tc初始化为Ta,然后加上dz*qg/ha,最后再加上Ta。

Tc=qg/ha+Ta; %地表温度 T0f1(1)=Tc; T0f2(1)=Tc; for j=1:num_layers dtt=dzqg/lamdpor(j); %竖向网格温差 if(j==1) j00=2; j01=nz/4+1; elseif(j==2) j00=nz/4+2; j01=2nz/4+1; elseif(j==3) j00=2nz/4+2; j01=3nz/4+1; else j00=3*nz/4+2; j01=nz+1; end for i=j00:j01 Tc=Tc+dtt; if(i<=nj) T0f1(i)=Tc; T0f2(i)=Tc; end end end优化代码

这段代码可以进行如下优化: 1. 尽量减少不必要的变量赋值,如将dtt=dzqg/lamdpor(j)提到循环外面赋值。 2. 减少if-else语句的使用,如使用向量化操作代替if-else语句。 3. 将常数提取出来,避免重复计算。 下面是优化后的程序: ```matlab dtt_all = dz * qg ./ lamdpor; j00 = [2, nz/4+2, 2*nz/4+2, 3*nz/4+2]; j01 = [nz/4+1, 2*nz/4+1, 3*nz/4+1, nz+1]; J = 1:num_layers; T0f1 = ones(nj, 1) * (qg/ha+Ta); T0f2 = T0f1; for j = J jj = j00(j):j01(j); Tc = T0f1(1) + sum(dtt_all(1:j-1)); T0f1(jj) = Tc + dtt_all(j-j00(j)+1) * (1:nj>=j00(j) & 1:nj<=j01(j))'; T0f2(jj) = Tc + dtt_all(j-j00(j)+1) * (1:nj>=j00(j) & 1:nj<=j01(j))'; end ``` 这样做的好处是:通过减少变量赋值和if-else语句的使用,可以使程序更加简洁易懂;通过将常数提取出来,避免了重复计算;通过使用向量化操作代替if-else语句,可以使程序运行速度更快。

相关推荐

这段代码中加一个test loss功能 class LSTM(nn.Module): def __init__(self, input_size, hidden_size, num_layers, output_size, batch_size, device): super().__init__() self.device = device self.input_size = input_size self.hidden_size = hidden_size self.num_layers = num_layers self.output_size = output_size self.num_directions = 1 # 单向LSTM self.batch_size = batch_size self.lstm = nn.LSTM(self.input_size, self.hidden_size, self.num_layers, batch_first=True) self.linear = nn.Linear(65536, self.output_size) def forward(self, input_seq): h_0 = torch.randn(self.num_directions * self.num_layers, self.batch_size, self.hidden_size).to(self.device) c_0 = torch.randn(self.num_directions * self.num_layers, self.batch_size, self.hidden_size).to(self.device) output, _ = self.lstm(input_seq, (h_0, c_0)) pred = self.linear(output.contiguous().view(self.batch_size, -1)) return pred if __name__ == '__main__': # 加载已保存的模型参数 saved_model_path = '/content/drive/MyDrive/危急值/model/dangerous.pth' device = 'cuda:0' lstm_model = LSTM(input_size=1, hidden_size=64, num_layers=1, output_size=3, batch_size=256, device='cuda:0').to(device) state_dict = torch.load(saved_model_path) lstm_model.load_state_dict(state_dict) dataset = ECGDataset(X_train_df.to_numpy()) dataloader = DataLoader(dataset, batch_size=256, shuffle=True, num_workers=0, drop_last=True) loss_fn = nn.CrossEntropyLoss() optimizer = optim.SGD(lstm_model.parameters(), lr=1e-4) for epoch in range(200000): print(f'epoch:{epoch}') lstm_model.train() epoch_bar = tqdm(dataloader) for x, y in epoch_bar: optimizer.zero_grad() x_out = lstm_model(x.to(device).type(torch.cuda.FloatTensor)) loss = loss_fn(x_out, y.long().to(device)) loss.backward() epoch_bar.set_description(f'loss:{loss.item():.4f}') optimizer.step() if epoch % 100 == 0 or epoch == epoch - 1: torch.save(lstm_model.state_dict(), "/content/drive/MyDrive/危急值/model/dangerous.pth") print("权重成功保存一次")

最新推荐

recommend-type

Linux 平台基于 Qt5 的网速浮窗.zip

Linux 平台基于 Qt5 的网速浮窗
recommend-type

手机游戏峡谷沼泽农田关卡地图Ai+EPS+PSD源文件.zip

游戏开发资源,游戏UI,游戏GUI,游戏图标,PSD格式,XD格式,PNG下载,源文件,可编辑下载,游戏购物充值界面,宝石,图标,PS格式,AI格式等,游戏APP
recommend-type

上市公司-企业资本结构动态调整数据及代码(2001-2022年).txt

数据存放网盘,txt文件内包含下载链接及提取码,永久有效。 样例数据及详细介绍参见文章:https://blog.csdn.net/li514006030/article/details/138324717
recommend-type

RTL8188FU-Linux-v5.7.4.2-36687.20200602.tar(20765).gz

REALTEK 8188FTV 8188eus 8188etv linux驱动程序稳定版本, 支持AP,STA 以及AP+STA 共存模式。 稳定支持linux4.0以上内核。
recommend-type

管理建模和仿真的文件

管理Boualem Benatallah引用此版本:布阿利姆·贝纳塔拉。管理建模和仿真。约瑟夫-傅立叶大学-格勒诺布尔第一大学,1996年。法语。NNT:电话:00345357HAL ID:电话:00345357https://theses.hal.science/tel-003453572008年12月9日提交HAL是一个多学科的开放存取档案馆,用于存放和传播科学研究论文,无论它们是否被公开。论文可以来自法国或国外的教学和研究机构,也可以来自公共或私人研究中心。L’archive ouverte pluridisciplinaire
recommend-type

:YOLOv1目标检测算法:实时目标检测的先驱,开启计算机视觉新篇章

![:YOLOv1目标检测算法:实时目标检测的先驱,开启计算机视觉新篇章](https://img-blog.csdnimg.cn/img_convert/69b98e1a619b1bb3c59cf98f4e397cd2.png) # 1. 目标检测算法概述 目标检测算法是一种计算机视觉技术,用于识别和定位图像或视频中的对象。它在各种应用中至关重要,例如自动驾驶、视频监控和医疗诊断。 目标检测算法通常分为两类:两阶段算法和单阶段算法。两阶段算法,如 R-CNN 和 Fast R-CNN,首先生成候选区域,然后对每个区域进行分类和边界框回归。单阶段算法,如 YOLO 和 SSD,一次性执行检
recommend-type

设计算法实现将单链表中数据逆置后输出。用C语言代码

如下所示: ```c #include <stdio.h> #include <stdlib.h> // 定义单链表节点结构体 struct node { int data; struct node *next; }; // 定义单链表逆置函数 struct node* reverse(struct node *head) { struct node *prev = NULL; struct node *curr = head; struct node *next; while (curr != NULL) { next
recommend-type

c++校园超市商品信息管理系统课程设计说明书(含源代码) (2).pdf

校园超市商品信息管理系统课程设计旨在帮助学生深入理解程序设计的基础知识,同时锻炼他们的实际操作能力。通过设计和实现一个校园超市商品信息管理系统,学生掌握了如何利用计算机科学与技术知识解决实际问题的能力。在课程设计过程中,学生需要对超市商品和销售员的关系进行有效管理,使系统功能更全面、实用,从而提高用户体验和便利性。 学生在课程设计过程中展现了积极的学习态度和纪律,没有缺勤情况,演示过程流畅且作品具有很强的使用价值。设计报告完整详细,展现了对问题的深入思考和解决能力。在答辩环节中,学生能够自信地回答问题,展示出扎实的专业知识和逻辑思维能力。教师对学生的表现予以肯定,认为学生在课程设计中表现出色,值得称赞。 整个课程设计过程包括平时成绩、报告成绩和演示与答辩成绩三个部分,其中平时表现占比20%,报告成绩占比40%,演示与答辩成绩占比40%。通过这三个部分的综合评定,最终为学生总成绩提供参考。总评分以百分制计算,全面评估学生在课程设计中的各项表现,最终为学生提供综合评价和反馈意见。 通过校园超市商品信息管理系统课程设计,学生不仅提升了对程序设计基础知识的理解与应用能力,同时也增强了团队协作和沟通能力。这一过程旨在培养学生综合运用技术解决问题的能力,为其未来的专业发展打下坚实基础。学生在进行校园超市商品信息管理系统课程设计过程中,不仅获得了理论知识的提升,同时也锻炼了实践能力和创新思维,为其未来的职业发展奠定了坚实基础。 校园超市商品信息管理系统课程设计的目的在于促进学生对程序设计基础知识的深入理解与掌握,同时培养学生解决实际问题的能力。通过对系统功能和用户需求的全面考量,学生设计了一个实用、高效的校园超市商品信息管理系统,为用户提供了更便捷、更高效的管理和使用体验。 综上所述,校园超市商品信息管理系统课程设计是一项旨在提升学生综合能力和实践技能的重要教学活动。通过此次设计,学生不仅深化了对程序设计基础知识的理解,还培养了解决实际问题的能力和团队合作精神。这一过程将为学生未来的专业发展提供坚实基础,使其在实际工作中能够胜任更多挑战。
recommend-type

"互动学习:行动中的多样性与论文攻读经历"

多样性她- 事实上SCI NCES你的时间表ECOLEDO C Tora SC和NCESPOUR l’Ingén学习互动,互动学习以行动为中心的强化学习学会互动,互动学习,以行动为中心的强化学习计算机科学博士论文于2021年9月28日在Villeneuve d'Asq公开支持马修·瑟林评审团主席法布里斯·勒菲弗尔阿维尼翁大学教授论文指导奥利维尔·皮耶昆谷歌研究教授:智囊团论文联合主任菲利普·普雷教授,大学。里尔/CRISTAL/因里亚报告员奥利维耶·西格德索邦大学报告员卢多维奇·德诺耶教授,Facebook /索邦大学审查员越南圣迈IMT Atlantic高级讲师邀请弗洛里安·斯特鲁布博士,Deepmind对于那些及时看到自己错误的人...3谢谢你首先,我要感谢我的两位博士生导师Olivier和Philippe。奥利维尔,"站在巨人的肩膀上"这句话对你来说完全有意义了。从科学上讲,你知道在这篇论文的(许多)错误中,你是我可以依
recommend-type

:YOLO目标检测算法的挑战与机遇:数据质量、计算资源与算法优化,探索未来发展方向

![:YOLO目标检测算法的挑战与机遇:数据质量、计算资源与算法优化,探索未来发展方向](https://img-blog.csdnimg.cn/7e3d12895feb4651b9748135c91e0f1a.png?x-oss-process=image/watermark,type_ZHJvaWRzYW5zZmFsbGJhY2s,shadow_50,text_Q1NETiBA5rKJ6YaJ77yM5LqO6aOO5Lit,size_20,color_FFFFFF,t_70,g_se,x_16) # 1. YOLO目标检测算法简介 YOLO(You Only Look Once)是一种