N的平方次循环换算N+X次循环

本文介绍了一种高效处理地理数据重复问题的方法,通过使用哈希表结构优化了数据处理流程,显著减少了循环次数。

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

class GeoInfo
  attr :lat,:lng
end
arr = GeoInfo.find(:all)

这里有个GeoInfo(地理信息)类,分别有两个属性lat,lng(纬度,经度)。

我获得了数据库中所有的地理数据,由于需求这些地理数据会出现大量坐标相同的情况,现在我根据需求,需要把位置相同的数据放在同一个数组中,最后获取一个新的数组。

arr_re = []
arr.each_with_index do |geo_a,index_a|
  fg = true
  arr_fg = [geo_a]
  arr.each_with_index do |geo_b,index_b|
    if geo_a.lat == geo_b.lat && geo_a.lng == geo_b.lng && index_a != index_b
      arr_fg << geo_b
      fg = false if index_a > index_b
    end
  end
  arr_re << arr_fg if fg
end
以上,代码是按照最初的想法实现,需要进行arr.size的平方次循环判断,需要消耗大量的时间来循环。如果数据有100条,那将要循环10000次才能拿到想要的结果。

map_fg = {}
arr.each do |a|
  if map_fg["#{a.lat};#{a.lng}"]
    map_fg["#{a.lat};#{a.lng}"] << a
  else
    map_fg["#{a.lat};#{a.lng}"] = [a]
  end
end
arr_re = []
for k,v in map_fg do
  arr_re << v
end
以上,代码使用hash存储,使用更大的空间换取循环次数的减少。如果有一百条数据,其中有10条重复数据,则只需要100+(100-10)次循环。

优劣立现。

model FanCoefficientModule record FanPressureCoefficients "风机流量-压力性能曲线系数记录" Real a "三项系数"; Real b "二项系数"; Real c "一项系数"; Real d "常数项"; end FanPressureCoefficients; function getPressureCoefficients "根据转速获取风机流量-压力性能曲线系数" input Real n "风机转速 (rpm)"; output FanPressureCoefficients coeff "风机流量-压力性能曲线系数"; // 已知转速点 (升序排列) protected constant Real[:] knownSpeeds = {1280, 1440, 1500, 1600}; // 风机流量-压力性能曲线系数数据库 [a,b,c,d] constant Real[:,4] pressureCoeffs = { {-3.267975749307142e-10, 5.2912277771354885e-06, -0.03403187464621037, 741.9041721827268}, {-2.422194407553435e-10, 3.7698924924279502e-06, -0.02559742162953265, 931.5430549263344}, {-2.3878745827961e-10, 3.921270379695461e-06, -0.02640953385995719, 1002.1261618905066}, {-2.2348965671136672e-10, 3.886724444389961e-06, -0.027754898177200903, 1147.392301537067}}; // 1280rpm // 1440rpm // 1500rpm // 1600rpm Integer idx "索引"; Real frac "插值比例"; Real factor "相似比例"; algorithm // 检查是否与已知转速匹配 for i in 1:size(knownSpeeds,1) loop if abs(n - knownSpeeds[i]) < 1e-6 then coeff := FanPressureCoefficients( pressureCoeffs[i,1], pressureCoeffs[i,2], pressureCoeffs[i,3], pressureCoeffs[i,4]); return; end if; end for; // 转速在范围内插值 if n > knownSpeeds[1] and n < knownSpeeds[size(knownSpeeds,1)] then // 寻找插值区间 for i in 1:size(knownSpeeds,1)-1 loop if n > knownSpeeds[i] and n < knownSpeeds[i+1] then frac := (n - knownSpeeds[i])/(knownSpeeds[i+1] - knownSpeeds[i]); coeff.a := pressureCoeffs[i,1] + frac*(pressureCoeffs[i+1,1] - pressureCoeffs[i,1]); coeff.b := pressureCoeffs[i,2] + frac*(pressureCoeffs[i+1,2] - pressureCoeffs[i,2]); coeff.c := pressureCoeffs[i,3] + frac*(pressureCoeffs[i+1,3] - pressureCoeffs[i,3]); coeff.d := pressureCoeffs[i,4] + frac*(pressureCoeffs[i+1,4] - pressureCoeffs[i,4]); return; end if; end for; // 转速超出范围使用相似定律 else if n < knownSpeeds[1] then idx := 1; // 使用最低转速点 else idx := size(knownSpeeds,1); // 使用最高转速点 end if; // 压力与转速平方成正比 factor := (n/knownSpeeds[idx])^2; coeff.a := pressureCoeffs[idx,1] * factor; coeff.b := pressureCoeffs[idx,2] * factor; coeff.c := pressureCoeffs[idx,3] * factor; coeff.d := pressureCoeffs[idx,4] * factor; end if; end getPressureCoefficients; model FanCoefficientSolution "风机系数获取模块" // 参数定义 parameter Real n;// 风机转速 (rpm) // 输出变量 output FanPressureCoefficients coeff; // "风机流量-压力性能曲线系数"; equation // 调用系数计算函数 coeff = getPressureCoefficients(n); end FanCoefficientSolution; annotation (Icon(coordinateSystem(preserveAspectRatio=false)), Diagram( coordinateSystem(preserveAspectRatio=false)), uses(Modelica(version="4.0.0"))); end FanCoefficientModule; 帮我检查一下这段modelica代码的问题
07-18
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

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

余额充值