2021-09-29
2021/9/29 23:40:46
本文主要是介绍2021-09-29,对大家解决编程问题具有一定的参考价值,需要的程序猿们随着小编来一起学习吧!
from collections import defaultdict states = ("Rainy", "Sunny") observations = ("Walk", "Shop", "Clean") start_probability = {"Rainy": 0.6, "Sunny": 0.4} transition_probability = { "Rainy": {"Rainy": 0.7, "Sunny": 0.3}, "Sunny": {"Rainy": 0.4, "Sunny": 0.6}, } emission_probability = { "Rainy": {"Walk": 0.1, "Shop": 0.4, "Clean": 0.5}, "Sunny": {"Walk": 0.6, "Shop": 0.3, "Clean": 0.1}, } def compute(obs, states, start_p, trans_p, emit_p): v = [{} for _ in range(len(obs))] path = defaultdict(list) for state in states: v[0][state] = start_p[state] * emit_p[state][obs[0]] path[state].append(state) for t in range(1, len(obs)): for y1 in states: max_prob = -1 for y0 in v[t - 1]: nprob = v[t - 1][y0] * trans_p[y0][y1] * emit_p[y1][obs[t]] if nprob > max_prob: max_prob = nprob max_state = y0 v[t][y1] = max_prob newpath = [] for state1 in path[max_state]: newpath.append(state1) newpath.append(y1) path[y1] = newpath prob = -1 for y1 in states: if v[len(obs) - 1][y1] > prob: prob = v[len(obs) - 1][y1] state = y1 return path[state] if __name__ == "__main__": max_path = compute(observations, states, start_probability, transition_probability, emission_probability) print(max_path)
这篇关于2021-09-29的文章就介绍到这儿,希望我们推荐的文章对大家有所帮助,也希望大家多多支持为之网!
- 2024-06-26结对编程到底难不难?答案在这里
- 2024-06-19《2023版Java工程师》课程升级公告
- 2024-06-15matplotlib作图不显示3D图,怎么办?
- 2024-06-1503-Loki 日志监控
- 2024-06-1504-让LLM理解知识 -Prompt
- 2024-06-05做软件测试需要懂代码吗?
- 2024-06-0514-ShardingSphere的分布式主键实现
- 2024-06-03为什么以及如何要进行架构设计权衡?
- 2024-05-31全网首发第二弹!软考2024年5月《软件设计师》真题+解析+答案!(11-20题)
- 2024-05-31全网首发!软考2024年5月《软件设计师》真题+解析+答案!(21-30题)