代码编织梦想

用sympy 求偏导

import numpy as np
import pandas as pd
import sympy as sy


theta, n, k = sy.symbols('theta, n, k')
z = theta+n+k
# print('z :{}'.format(z))
# result = z.subs({theta: 0, n: 1})  # 用数值分别对x、y进行替换
# print(result)
#
dx = sy.diff(z, theta)  # 对x求偏导
# print(dx)
# result = dx.subs({theta: 0, n: 2, k: 1})
# print(result)
#
dy = sy.diff(z, n)  # 对y求偏导
# print(dy)
# result = dy.subs({theta: 0, n: 2, k: 1})
# print(result)

dz = sy.diff(z, k)  # 对y求偏导
print(dz)
theta_set = [0, 0.7854, 1.5708, 2.3561, 3.1415]
n_set = [0, 0.5, 1, 2, 3, 8]
k_set = [1, 2, 3, 10]
result = np.zeros([len(k_set), len(n_set), len(theta_set)])
for k1 in k_set:
    for n1 in n_set:
        for theta1 in theta_set:
            # save = dz.subs({theta: theta1, n: n1, k: k1})
            # save = dx.subs({theta: theta1, n: n1, k: k1})
            save = dy.subs({theta: theta1, n: n1, k: k1})
            result[k_set.index(k1)][n_set.index(n1)][theta_set.index(theta1)] = save
for i in range(len(k_set)):
    df = pd.DataFrame(result[i])
    df.to_csv('result.csv', mode='a')

# result = dz.subs({theta: 3.14, n: 1, k: 1})
# print(result)

版权声明:本文为博主原创文章,遵循 CC 4.0 BY-SA 版权协议,转载请附上原文出处链接和本声明。 本文链接: https://blog.csdn.net/m0_38088084/article/details/109277746

Python学习(十):字典的魔法及常用函数,以及布尔值-爱代码爱编程

Python学习(十):字典的魔法 这是我关于Python学习的第十篇博客,加油! get、update、keys、values、items 字典 字典dict基本结构 info={‘k1’:‘v1’,‘k2’:‘v2’} 其中‘k1’:'v1’为一个键值对,字典中键值对用逗号分割字典可以嵌套,value可以是任何值info={‘k1’:True,

树莓派开发之——跑马灯-爱代码爱编程

raspberry实现跑马灯用例 使用跑马灯是学习树莓派的开始,下面的例子即为树莓派开发跑马灯的python代码用例。 import RPi.GPIO as GPIO import time WHITE = 5 YELLOW = 6 BLUE = 13 GREEN = 19 RED = 26 leds = [WHITE, YELLOW, BLUE,

爬虫-requests 使用百度翻译-爱代码爱编程

import requests import json if __name__ == '__main__': #1. 指定url post_url = 'https://fanyi.baidu.com/sug' #2. UA伪装 headers = { 'User Agent':'Mozilla/5.0 (W

python函数定义,局部变量、全局变量-爱代码爱编程

函数 使用函数的意义在于提高编程效率,避免大量的重复工作 函数分为: 内置函数:可以直接调用的函数 第三方,模块相关函数:即导入第三方包,直接调用第三方包函数 自定义函数:即根据自己的计算需求,编写属于自己的程序段 常用内置函数 (1)len() :计算字符串的长度 如` a =list((1,2,3,4,5)) print(len(a)) (2

PTA:Python字符串字母大小写转换-爱代码爱编程

第3章-14 字符串字母大小写转换 (15分) 本题要求编写程序,对一个以“#”结束的字符串,将其小写字母全部转换成大写字母,把大写字母全部转换成小写字母,其他字符不变输出。 输入格式: 本题要求编写程序,对一个以“#”结束的字符串,将其小写字母全部转换成大写字母,把大写字母全部转换成小写字母,其他字符不变输出。 输出格式: 在一行中输出大小写转

【Practical】Python-Numpy 计算二分类问题准确率与召回率-爱代码爱编程

文章目录 问题描述.实现思路.实际代码. 问题描述. 【情境】监督学习中的二分类问题(正反例以1-0标记),已经训练好了模型并且获得了预测值向量 Y p